Changeset a52369 for src/Fragmentation/Exporters/SphericalPointDistribution.cpp

Timestamp:

May 25, 2016, 7:13:59 AM (9 years ago)

Author:

Frederik Heber <heber@…>

Children:

9cf90e

Parents:

46e561

git-author:

Frederik Heber <heber@…> (05/18/16 08:23:05)

git-committer:

Frederik Heber <heber@…> (05/25/16 07:13:59)

Message:

FIX: Need to check whether old and newcenter are not linear dependent.

File:

• src/Fragmentation/Exporters/SphericalPointDistribution.cpp (modified) (1 diff)

src/Fragmentation/Exporters/SphericalPointDistribution.cpp

@@ -704 +704 @@
     // calculate rotation axis and angle
     Rotation.first = oldCenter;
-    Rotation.first.VectorProduct(newCenter);
-    Rotation.first.Normalize();
+    if (oldCenter.IsParallelTo(newCenter, 1e-6))
+      Rotation.first.GetOneNormalVector(oldCenter);
+    else {
+      Rotation.first.VectorProduct(newCenter);
+      Rotation.first.Normalize();
+    }
     // construct reference vector to determine direction of rotation
     const double sign = determineSignOfRotation(newCenter, oldCenter, Rotation.first);