Changeset 10398

Timestamp:

11/01/11 10:42:52 (13 years ago)

Author:

Mathieu Morlighem

Message:

Fixed referencial consistency

File:

• issm/trunk/src/m/classes/diagnostic.m (modified) (1 diff)

issm/trunk/src/m/classes/diagnostic.m

@@ -113 +113 @@
                         if any(sum(isnan(md.diagnostic.referential),2)==0),
                                 pos=find(sum(isnan(md.diagnostic.referential),2)==0);
-                                if any(dot(md.diagnostic.referential(pos,1:3),md.diagnostic.referential(pos,4:6))),
-                                        dot(md.diagnostic.referential(pos,1:3),md.diagnostic.referential(pos,4:6))
-                                        checkmessage(['model ' md.miscellaneous.name ' has problem with rotated spc. Vectors in diagnostic.referential (colums 1 to 3 and 4 to 6) must be orthogonal']);
-                                end
-                        end
-                        %CHECK THAT ROTATION IS IN THE (X,Y) PLANE FOR 2D MODELS
-                        if md.mesh.dimension==2,
-                                pos=find(sum(isnan(md.diagnostic.referential),2)==0  & md.flowequation.vertex_equation==2);
-                                if any(md.diagnostic.referential(pos,3:5)~=0);
-                                        checkmessage(['model ' md.miscellaneous.name ' has problem with rotated spc. The rotation should be in the (x,y) plane for 2D diagnostic models (nodeonmacayeal)']);
+                                if any(abs(dot(md.diagnostic.referential(pos,1:3)',md.diagnostic.referential(pos,4:6)'))>eps),
+                                        checkmessage(['model ' md.miscellaneous.name ' has problem with referential. Vectors in diagnostic.referential (colums 1 to 3 and 4 to 6) must be orthogonal']);
                                 end
                         end