update acos issues

jdtuck · Jan 29, 2021 · 7ef4b36 · 7ef4b36
1 parent 3d69d38
commit 7ef4b36
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 6 deletions.
diff --git a/Find_Rotation_and_Seed_unique.m b/Find_Rotation_and_Seed_unique.m
@@ -69,7 +69,13 @@
     else
         q2new  = q2n;
     end
-    Ec = acos(InnerProd_Q(q1,q2new));
+
+    tmp = InnerProd_Q(q1,q2new);
+    if tmp > 1
+        tmp = 1;
+    end
+
+    Ec = acos(tmp);
     if Ec < minE
         Rbest = R;
         q2best = q2new;

diff --git a/fdacurve.m b/fdacurve.m
@@ -25,10 +25,10 @@
         cent            % center
         scale           % scale
         E               % energy
-        len             % scale of curve
-        len_q           % scale of SRVF
-        mean_scale      % mean scale
-        mean_scale_q    % mean scale
+        len             % length of curves
+        len_q           % length of SRVFs
+        mean_scale      % mean length
+        mean_scale_q    % mean length SRVF
     end
 
     methods
@@ -610,7 +610,7 @@ function plot_pca(obj, n)
 
                     p = q_to_curve(q2n,tmp_scale);
                     if obj.scale
-                        mv = 0.2*mean(obj.len);
+                        mv = 0.2*obj.mean_scale;
                     else
                         mv = 0.2;
                     end