crossproduct Function

private elemental function crossproduct(lhs, rhs) result(cross)

Compute the cross product.

$$\vec V=\left({y_1 z_2 - z_1 y_2}\right)\vec i + \left({z_1 x_2 - x_1 z_2}\right)\vec j + \left({x_1 y_2 - y_1 x_2}\right)\vec k$$ where $x_i$, $y_i$ and $z_i$ $i=1,2$ are the components of the vectors.

 type(vector) :: pt(0:2)
 pt(1) = 2 * ex
 pt(2) = ex
 pt(0) = pt(1).cross.pt(2)
 print "(3(F3.1,1X))", abs(pt(0)%x), abs(pt(0)%y), abs(pt(0)%z)

Arguments

Return Value type(vector)

Cross product vector.