triv_lib/geomat4d.c:180 void TrivVectCross3Vecs(TrivVType A, TrivVType B, TrivVType C, TrivVType Res)
with(linalg);
readlib(C);
d := det( matrix( [ [I, J, K, L],
[A[0], A[1], A[2], A[3]],
[B[0], B[1], B[2], B[3]],
[C[0], C[1], C[2], C[3]] ] ) );
coeff( d, I );
coeff( d, J );
coeff( d, K );
coeff( d, L );
A, B, C: The three vectors to compute their cross product.
Res: Where the output goes into.