Changeset 16195 for issm/trunk-jpl/test/NightlyRun/test3019.m
- Timestamp:
- 09/19/13 21:50:25 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test3019.m
r16192 r16195 4 4 md=parameterize(md,'../Par/SquareShelfConstrained.par'); 5 5 md=setflowequation(md,'SSA','all'); 6 md.cluster=generic('name',oshostname(),'np',1); 7 md.toolkits.DefaultAnalysis=issmgslsolver(); 6 md.cluster=generic('name',oshostname(),'np',3); 7 %md.debug.valgrind=true; 8 %md.verbose=verbose('11111111') 9 md.toolkits.DefaultAnalysis=issmmumpssolver(); 10 %md.toolkits.DefaultAnalysis.solver_type='gsl'; 8 11 9 12 md.autodiff.isautodiff=true; … … 20 23 jac_reverse=md.results.TransientSolution(1).AutodiffJacobian; 21 24 22 %now run vectorial forward mode23 md.autodiff.independents={independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices,'fov_forward_indices',(1:md.mesh.numberofvertices)')};24 md.autodiff.dependents={dependent('name','MaxVel','type','scalar')};25 md.autodiff.driver='fov_forward';26 27 md=solve(md,TransientSolutionEnum());28 29 %recover jacobian:30 jac_forward=md.results.TransientSolution(1).AutodiffJacobian;31 32 25 %Fields and tolerances to track changes 33 field_names ={'Jac Forward','Jac Reverse','Jac Forward -Reverse'};34 field_tolerances={1e-8 ,1e-8,1e-6};35 field_values={jac_ forward,jac_reverse,jac_forward-jac_reverse};26 field_names ={'Jac Reverse'}; 27 field_tolerances={1e-8}; 28 field_values={jac_reverse};
Note:
See TracChangeset
for help on using the changeset viewer.