Exception for HairBSDF, Mp, upper

\[\left(\left(\left(\left(\left(-1 \leq cosTheta_i \land cosTheta_i \leq 1\right) \land \left(-1 \leq cosTheta_O \land cosTheta_O \leq 1\right)\right) \land \left(-1 \leq sinTheta_i \land sinTheta_i \leq 1\right)\right) \land \left(-1 \leq sinTheta_O \land sinTheta_O \leq 1\right)\right) \land 0.1 < v\right) \land v \leq 1.5707964\]

Math

\[\frac{e^{-\frac{sinTheta_i \cdot sinTheta_O}{v}} \cdot \frac{cosTheta_i \cdot cosTheta_O}{v}}{\left(\sinh \left(\frac{1}{v}\right) \cdot 2\right) \cdot v} \]

FPCore

(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
 :precision binary64
 (/
  (* (exp (- (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v))
  (* (* (sinh (/ 1.0 v)) 2.0) v)))

C

double code(double cosTheta_i, double cosTheta_O, double sinTheta_i, double sinTheta_O, double v) {
	return (exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinh((1.0 / v)) * 2.0) * v);
}

Fortran

real(8) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
    real(8), intent (in) :: costheta_i
    real(8), intent (in) :: costheta_o
    real(8), intent (in) :: sintheta_i
    real(8), intent (in) :: sintheta_o
    real(8), intent (in) :: v
    code = (exp(-((sintheta_i * sintheta_o) / v)) * ((costheta_i * costheta_o) / v)) / ((sinh((1.0d0 / v)) * 2.0d0) * v)
end function

Java

public static double code(double cosTheta_i, double cosTheta_O, double sinTheta_i, double sinTheta_O, double v) {
	return (Math.exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((Math.sinh((1.0 / v)) * 2.0) * v);
}

Python

def code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v):
	return (math.exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((math.sinh((1.0 / v)) * 2.0) * v)

Julia

function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
	return Float64(Float64(exp(Float64(-Float64(Float64(sinTheta_i * sinTheta_O) / v))) * Float64(Float64(cosTheta_i * cosTheta_O) / v)) / Float64(Float64(sinh(Float64(1.0 / v)) * 2.0) * v))
end

MATLAB

function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
	tmp = (exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinh((1.0 / v)) * 2.0) * v);
end

Wolfram

code[cosTheta$95$i_, cosTheta$95$O_, sinTheta$95$i_, sinTheta$95$O_, v_] := N[(N[(N[Exp[(-N[(N[(sinTheta$95$i * sinTheta$95$O), $MachinePrecision] / v), $MachinePrecision])], $MachinePrecision] * N[(N[(cosTheta$95$i * cosTheta$95$O), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Sinh[N[(1.0 / v), $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision] * v), $MachinePrecision]), $MachinePrecision]