pub fn compute_z(&self, r_ij: V3) -> (f64, V3) {
        let z_hat = V3([0.0, 0.0, 1.0]);
        let Output {
            value, grad_rij,
            grad_ni: _,
            grad_nj: _,
        } = self.compute(r_ij, z_hat, z_hat);

        (value, grad_rij)
    }
/// Computes crespi assuming normals are all +Z.
Computes crespi assuming normals are all +Z.
Computes crespi assuming normals are all + Z .
pub fn compute_z(&self, r_ij: V3) -> (f64, V3) {
        let z_hat = V3([0.0, 0.0, 1.0]);
        let Output {
            value, grad_rij,
            grad_ni: _,
            grad_nj: _,
        } = self.compute(r_ij, z_hat, z_hat);
        (value, grad_rij)
    }
pub fn compute_z ( & self , r_ij : V3 ) -> ( f64 , V3 ) { let z_hat = V3 ( [ 0.0 , 0.0 , 1.0 ] ) ; let Output { value , grad_rij , grad_ni : _ , grad_nj : _ , } = self . compute ( r_ij , z_hat , z_hat ) ; ( value , grad_rij ) }
Computes crespi assuming normals are all +Z.
Computes crespi assuming normals are all + Z .
&self
r_ij
V3
&self
r_ij
V3