xp(2.*a);
    return (a - 1.) / (a + 1.);
  }
  vec4 tanh_(vec4 v) {
    v = clamp(v, -10., 10.);
    v = exp(2.*v);
    return (v - 1.) / (v + 1.);
  }