naushad-rahman / ellipsoids

Automatically exported from code.google.com/p/ellipsoids
0 stars 0 forks source link

Shape matrix symmetry #3

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago
I am using the ellipsoidal constructor and found the shape matrix test to 
be sensitive to numerical precision, resulting in the "ELLIPSOID: shape 
matrix must be symmetric." due to a 10e-19 difference in the non-diagonal 
elements.

Suggested fix:
Replace the second condition in (line 74 in \@ellipsoid\ellipsoid.m)

%if (m ~= n) | (min(min((Q == Q'))) == 0)

by

norm(Q-Q')>= ellOptions.abs_tol

Original issue reported on code.google.com by pvazteix...@gmail.com on 18 Jun 2009 at 3:23

Attachments: