Lines Matching refs:f2
686 double* f1, *f2;
722 f2 = v + 8*9;
724 // f1, f2 is a basis => lambda*f1 + mu*f2 is an arbitrary f. matrix.
726 // so f ~ lambda*f1 + (1 - lambda)*f2.
727 // use the additional constraint det(f) = det(lambda*f1 + (1-lambda)*f2) to find lambda.
731 f1[i] -= f2[i];
733 t0 = f2[4]*f2[8] - f2[5]*f2[7];
734 t1 = f2[3]*f2[8] - f2[5]*f2[6];
735 t2 = f2[3]*f2[7] - f2[4]*f2[6];
737 c[3] = f2[0]*t0 - f2[1]*t1 + f2[2]*t2;
740 f1[3]*(f2[1]*f2[8] - f2[2]*f2[7]) +
741 f1[4]*(f2[0]*f2[8] - f2[2]*f2[6]) -
742 f1[5]*(f2[0]*f2[7] - f2[1]*f2[6]) +
743 f1[6]*(f2[1]*f2[5] - f2[2]*f2[4]) -
744 f1[7]*(f2[0]*f2[5] - f2[2]*f2[3]) +
745 f1[8]*(f2[0]*f2[4] - f2[1]*f2[3]);
751 c[1] = f2[0]*t0 - f2[1]*t1 + f2[2]*t2 -
752 f2[3]*(f1[1]*f1[8] - f1[2]*f1[7]) +
753 f2[4]*(f1[0]*f1[8] - f1[2]*f1[6]) -
754 f2[5]*(f1[0]*f1[7] - f1[1]*f1[6]) +
755 f2[6]*(f1[1]*f1[5] - f1[2]*f1[4]) -
756 f2[7]*(f1[0]*f1[5] - f1[2]*f1[3]) +
757 f2[8]*(f1[0]*f1[4] - f1[1]*f1[3]);
771 double s = f1[8]*r[k] + f2[8];
784 fmatrix[i] = f1[i]*lambda + f2[i]*mu;
