Prev Next det_33.hpp Headings

Source: det_33
# ifndef CPPAD_DET_33_INCLUDED
# define CPPAD_DET_33_INCLUDED
 
# include <cppad/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