Source: det_33
Keywords: source: det_33 source

| Prev | Next |  |  |  |  |  | det_33.hpp | Headings

Source: det_33
`# ifndef CPPAD_DET_33_HPP `
` # define CPPAD_DET_33_HPP

    # include <cppad/utility/near_equal.hpp>
    namespace CppAD {
    template <class Vector>
    	bool det_33(const Vector &x, const Vector &d)
    	{	bool ok = true;

    		// use expansion by minors to compute the determinant by hand
    		double check = 0.;
    		check += x[0] * ( x[4] * x[8] - x[5] * x[7] );
    		check -= x[1] * ( x[3] * x[8] - x[5] * x[6] );
    		check += x[2] * ( x[3] * x[7] - x[4] * x[6] );

    		ok &= CppAD::NearEqual(check, d[0], 1e-10, 1e-10);

    		return ok;
    	}
    }

# endif`

Input File: omh/det_33_hpp.omh