source: issm/oecreview/Archive/13393-13976/ISSM-13507-13508.diff

../trunk-jpl/test/NightlyRun/ad.m

@@ -7 +7 @@
 md.autodiff.isautodiff=true;
 md.verbose.autodiff=true;
 md.autodiff.independents={...
-        independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices,'fov_forward_indices',(1:md.mesh.numberofvertices)')
+        independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices);
+        %independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices,'fov_forward_indices',(1:md.mesh.numberofvertices)')
         %independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices,'fos_forward_index',1)
         };
 
 md.autodiff.dependents={...
-        dependent('name','IceVolume','type','scalar')...
-        dependent('name','MaxVel','type','scalar')...
-        dependent('name','MassFlux','type','scalar','exp','../Exp/MassFlux1.exp','mesh',md.mesh,'index',1)...
-        dependent('name','MassFlux','type','scalar','exp','../Exp/MassFlux2.exp','mesh',md.mesh,'index',2)...
+        dependent('name','MaxVel','type','scalar','fos_reverse_index',1)...
+        %dependent('name','IceVolume','type','scalar')...
+        %dependent('name','MassFlux','type','scalar','exp','../Exp/MassFlux1.exp','mesh',md.mesh,'index',1)...
+        %dependent('name','MassFlux','type','scalar','exp','../Exp/MassFlux2.exp','mesh',md.mesh,'index',2)...
         };
 
 %md.autodiff.driver='fos_forward';
-md.autodiff.driver='fov_forward';
+md.autodiff.driver='fos_reverse';
 
 md=solve(md,TransientSolutionEnum);
 

../trunk-jpl/test/NightlyRun/test3019.m

@@ +1 @@
+md=triangle(model,'../Exp/Square.exp',100000);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setflowequation(md,'macayeal','all');
+md.cluster=generic('name',oshostname(),'np',3);
+
+md.autodiff.isautodiff=true;
+md.verbose.autodiff=true;
+
+%first run scalar reverse mode: 
+md.autodiff.independents={independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices)};
+md.autodiff.dependents={dependent('name','MaxVel','type','scalar','fos_reverse_index',1)};
+md.autodiff.driver='fos_reverse';
+
+md=solve(md,TransientSolutionEnum);
+
+%recover jacobian: 
+jac_reverse=md.results.TransientSolution(1).AutodiffJacobian;
+
+
+%now run vectorial forward mode
+md.autodiff.independents={independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices,'fov_forward_indices',(1:md.mesh.numberofvertices)')};
+md.autodiff.dependents={dependent('name','MaxVel','type','scalar')};
+md.autodiff.driver='fov_forward';
+
+md=solve(md,TransientSolutionEnum);
+
+%recover jacobian: 
+jac_forward=md.results.TransientSolution(1).AutodiffJacobian;
+
+
+%Fields and tolerances to track changes
+field_names     ={'Jac Forward','Jac Reverse','Jac Forward - Reverse'};
+field_tolerances={1e-13,1e-13,1e-13};
+field_values={jac_forward,jac_reverse,jac_forward-jac_reverse};