Changeset 10430

Timestamp:

11/02/11 11:47:08 (13 years ago)

Author:

Mathieu Morlighem

Message:

Better test for rotation: only change CS on the boudary

Location:

issm/trunk/test

Files:

• Archives/Archive1601.mat (modified) ( previous)
• NightlyRun/test1601.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test1601.m

@@ -27 +27 @@
 disp(['Error between Cartesian and rotated CS: ' num2str(max(abs(vel0-vel1))/(max(abs(vel0))+eps)) ]);
 
+%Now, put CS back to normal except on the side where the spc are applied
+pos=find(x==0 | x==1000000);
+md.diagnostic.referential(:)=NaN;
+md.diagnostic.referential(pos,1:3)=repmat([cos(theta),sin(theta),0],size(pos,1),1);
+md.diagnostic.referential(pos,4:6)=repmat([0,0,1],size(pos,1),1);
+md=solve(md,DiagnosticSolutionEnum);
+vel2=md.results.DiagnosticSolution.Vel;
+
+plotmodel(md,'data',vel0,'data',vel2,'data',vel2-vel0,'title','Cartesian CS','title','Rotated CS','title','difference')
+disp(['Error between Cartesian and rotated CS: ' num2str(max(abs(vel0-vel2))/(max(abs(vel0))+eps)) ]);
+
 %Fields and tolerances to track changes
-field_names     ={'vel1'};
-field_tolerances={1e-11};
+field_names     ={'vel1','vel2'};
+field_tolerances={1e-11,1e-11};
 field_values={...
         vel1, ...
+        vel2, ...
         };