Commit 057ddb96 by Francois Gygi

fixed sign of torsion angle


git-svn-id: http://qboxcode.org/svn/qb/trunk@613 cba15fb0-1239-40c8-b417-11db7ca47a34
parent e9f5cfb7
...@@ -3,7 +3,7 @@ ...@@ -3,7 +3,7 @@
// TorsionCmd.h: // TorsionCmd.h:
// //
//////////////////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////////////////
// $Id: TorsionCmd.h,v 1.3 2008-02-12 05:39:18 fgygi Exp $ // $Id: TorsionCmd.h,v 1.4 2008-05-23 03:23:36 fgygi Exp $
#ifndef TORSIONCMD_H #ifndef TORSIONCMD_H
#define TORSIONCMD_H #define TORSIONCMD_H
...@@ -103,7 +103,7 @@ class TorsionCmd : public Cmd ...@@ -103,7 +103,7 @@ class TorsionCmd : public Cmd
D3vector e123 = normalized(e12^e32); D3vector e123 = normalized(e12^e32);
D3vector e234 = normalized(e23^e43); D3vector e234 = normalized(e23^e43);
double cc = max(min(e123*e234,1.0),-1.0); double cc = max(min(e123*e234,1.0),-1.0);
double ss = max(min((e123^e234)*e32,1.0),-1.0); double ss = -max(min((e123^e234)*e32,1.0),-1.0);
a = (180.0/M_PI) * atan2(ss,cc); a = (180.0/M_PI) * atan2(ss,cc);
} }
......
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