Example #1
0
bool testDet(Mat_DP* A, Vec_CPLX_DP* eigens, double tol) {
    std::complex<double> prodOfEigens = 0;
    for (int i = 0; i < A->nrows(); i++) {
        prodOfEigens *= (*eigens)[i];
    }
    double det = detmat((*A));
    return (abs(det - prodOfEigens) > tol);
}
Example #2
0
void invmat(double *mat, double *imat)

/* ======== Inverse of a 3x3 matrix ========= */
/*	DO NOT MAKE THE OUTPUT THE SAME AS ONE OF THE INPUTS!! */

{
double det;
int count;

det = detmat(mat);	/* Determinant */
adjmat(mat, imat);	/* Adjoint */

for (count=0; count<9; count++)
	*imat++ /= det;		
}