Commit 967bd17d by Francois Gygi

Define cell volume as fabs(a0*(a1^a2))

git-svn-id: http://qboxcode.org/svn/qb/trunk@913 cba15fb0-1239-40c8-b417-11db7ca47a34
parent 30d854df
......@@ -37,7 +37,7 @@ void UnitCell::set(const D3vector& a0, const D3vector& a1, const D3vector& a2)
amat_[8] = a2.z;
// volume = det(A)
volume_ = a0 * ( a1 ^ a2 );
volume_ = fabs(a0 * ( a1 ^ a2 ));
if ( volume_ > 0.0 )
{
// Compute rows of A-1 (columns of A^-T)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment