Changeset 16202
- Timestamp:
- 09/20/13 09:54:36 (12 years ago)
- Location:
- issm/trunk-jpl/test/NightlyRun
- Files:
-
- 4 deleted
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test3119.m
r16194 r16202 1 %test reverse scalar vs forward vectorial driversin ADOLC, using the test3009 setup, equivalent to test109 setup.1 %test reverse scalar driver in ADOLC, using the test3009 setup, equivalent to test109 setup. 2 2 md=triangle(model(),'../Exp/Square.exp',100000.); 3 3 md=setmask(md,'all',''); … … 20 20 jac_reverse=md.results.TransientSolution(1).AutodiffJacobian; 21 21 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 22 %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};23 field_names ={'Jac Reverse'}; 24 field_tolerances={1e-8}; 25 field_values={jac_reverse,}; -
issm/trunk-jpl/test/NightlyRun/test3119.py
r16194 r16202 32 32 jac_reverse=md.results.TransientSolution[0].AutodiffJacobian 33 33 34 #now run vectorial forward mode35 md.autodiff.independents=[independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices,'fov_forward_indices',numpy.arange(0,md.mesh.numberofvertices))]36 md.autodiff.dependents=[dependent('name','MaxVel','type','scalar')]37 md.autodiff.driver='fov_forward'38 39 md=solve(md,TransientSolutionEnum())40 41 #recover jacobian:42 jac_forward=md.results.TransientSolution[0].AutodiffJacobian43 44 34 #Fields and tolerances to track changes 45 field_names =['Jac Forward','Jac Reverse','Jac Forward -Reverse']46 field_tolerances=[1e-13 ,1e-13,1e-13]47 field_values=[jac_ forward,jac_reverse,jac_forward-jac_reverse]35 field_names =['Jac Reverse'] 36 field_tolerances=[1e-13] 37 field_values=[jac_reverse]
Note:
See TracChangeset
for help on using the changeset viewer.