--- branches/mmeineke/mdtools/md_code/DirectionalAtom.cpp 2002/07/09 18:40:59 10 +++ trunk/mdtools/md_code/DirectionalAtom.cpp 2002/07/16 21:10:49 30 @@ -39,6 +39,22 @@ void DirectionalAtom::setQ( double the_q[4] ){ Azz = q0Sqr - q1Sqr -q2Sqr +q3Sqr; } +void DirectionalAtom::getA( double the_A[3][3] ){ + + A[0][0] = Axx; + A[0][1] = Axy; + A[0][2] = Axz; + + A[1][0] = Ayx; + A[1][1] = Ayy; + A[1][2] = Ayz; + + A[2][0] = Azx; + A[2][1] = Azy; + A[2][2] = Azz; +} + + void DirectionalAtom::getU( double the_u[3] ){ the_u[0] = sux;