source:
issm/oecreview/Archive/13393-13976/ISSM-13517-13518.diff
Last change on this file was 13980, checked in by , 12 years ago | |
---|---|
File size: 5.4 KB |
-
../trunk-jpl/test/NightlyRun/test212.py
1 """ 2 == == == == == == == == == == == == == == == == == == == 3 Auto generated python script for ISSM: test212.m 4 Created on 2012-09-27 via translateToPy.py Ver 1.0 by mikep 5 == == == == == == == == == == == == == == == == == == == 6 7 Matlab script conversion into python 8 translateToPy.py Author: Michael Pellegrin 9 translateToPy.py Date: 09/24/12 10 == == == == == == == == == == == == == == == == == == == 11 """ 12 13 from MatlabFuncs import * 14 from model import * 15 from EnumDefinitions import * 16 from numpy import * 17 from triangle import * 18 from setmask import * 19 from parameterize import * 20 from setflowequation import * 21 from solve import * 22 23 md=triangle(model(),'../Exp/Square.exp',200000) 24 md=setmask(md,'all','') 25 md=parameterize(md,'../Par/SquareShelf.py') 26 md=setflowequation(md,'macayeal','all') 27 28 29 # control parameters 30 31 md.inversion.iscontrol=1 32 md.inversion.control_parameters=['MaterialsRheologyBbar'] 33 md.inversion.min_parameters=10**6*ones((md.mesh.numberofvertices,1)) 34 md.inversion.max_parameters=2*10**9*ones((md.mesh.numberofvertices,1)) 35 md.inversion.nsteps=2 36 md.inversion.cost_functions=101*ones((md.inversion.nsteps,1)) 37 md.inversion.cost_functions_coefficients=ones((md.mesh.numberofvertices,1)) 38 md.inversion.gradient_scaling=10**7*ones((md.inversion.nsteps,1)) 39 md.inversion.maxiter_per_step=2*ones((md.inversion.nsteps,1)) 40 md.inversion.step_threshold=0.3*ones((md.inversion.nsteps,1)) 41 md.inversion.vx_obs=md.initialization.vx 42 md.inversion.vy_obs=md.initialization.vy 43 44 45 md.cluster=generic('name',oshostname(),'np',3) 46 md=solve(md,DiagnosticSolutionEnum()) 47 48 49 # Fields and tolerances to track changes 50 51 field_names =['Gradient','Misfits','MaterialsRheologyBbar','Pressure','Vel','Vx','Vy'] 52 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13] 53 field_values=[\ 54 md.results['DiagnosticSolution'][1]['Gradient1'],\ 55 md.results['DiagnosticSolution'][1]['J'],\ 56 md.results['DiagnosticSolution'][1]['MaterialsRheologyBbar'],\ 57 md.results['DiagnosticSolution'][1]['Pressure'],\ 58 md.results['DiagnosticSolution'][1]['Vel'],\ 59 md.results['DiagnosticSolution'][1]['Vx'],\ 60 md.results['DiagnosticSolution'][1]['Vy'],\ 61 ] -
../trunk-jpl/test/NightlyRun/test212.m
24 24 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 25 25 field_values={... 26 26 (md.results.DiagnosticSolution.Gradient1),... 27 md.results.DiagnosticSolution.J,...27 (md.results.DiagnosticSolution.J),... 28 28 (md.results.DiagnosticSolution.MaterialsRheologyBbar),... 29 29 (md.results.DiagnosticSolution.Pressure),... 30 30 (md.results.DiagnosticSolution.Vel),... -
../trunk-jpl/test/NightlyRun/InNeedOfDebugging/test212.py
1 """2 == == == == == == == == == == == == == == == == == == ==3 Auto generated python script for ISSM: test212.m4 Created on 2012-09-27 via translateToPy.py Ver 1.0 by mikep5 == == == == == == == == == == == == == == == == == == ==6 7 Matlab script conversion into python8 translateToPy.py Author: Michael Pellegrin9 translateToPy.py Date: 09/24/1210 == == == == == == == == == == == == == == == == == == ==11 """12 13 from MatlabFuncs import *14 from model import *15 from EnumDefinitions import *16 from numpy import *17 from triangle import *18 from setmask import *19 from parameterize import *20 from setflowequation import *21 from solve import *22 23 md=triangle(model(),'../Exp/Square.exp',200000)24 md=setmask(md,'all','')25 md=parameterize(md,'../Par/SquareShelf.py')26 md=setflowequation(md,'macayeal','all')27 28 29 # control parameters30 31 md.inversion.iscontrol=132 md.inversion.control_parameters=['MaterialsRheologyBbar']33 md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1)34 md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1)35 md.inversion.nsteps=236 md.inversion.cost_functions=101*ones(md.inversion.nsteps,1)37 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1)38 md.inversion.gradient_scaling=10^7*ones(md.inversion.nsteps,1)39 md.inversion.maxiter_per_step=2*ones(md.inversion.nsteps,1)40 md.inversion.step_threshold=0.3*ones(md.inversion.nsteps,1)41 md.inversion.vx_obs=md.initialization.vx42 md.inversion.vy_obs=md.initialization.vy43 44 45 md.cluster=generic('name',oshostname(),'np',3)46 md=solve(md,DiagnosticSolutionEnum())47 48 49 # Fields and tolerances to track changes50 51 field_names =['Gradient' 'Misfits' 'MaterialsRheologyBbar' 'Pressure' 'Vel' 'Vx' 'Vy']52 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]53 field_values=[\54 md.results['DiagnosticSolution'][1]['Gradient1'],\55 md.results['DiagnosticSolution.J'],\56 md.results['DiagnosticSolution'][1]['MaterialsRheologyBbar'],\57 md.results['DiagnosticSolution'][1]['Pressure'],\58 md.results['DiagnosticSolution'][1]['Vel'],\59 md.results['DiagnosticSolution'][1]['Vx'],\60 md.results['DiagnosticSolution'][1]['Vy'],\61 ]
Note:
See TracBrowser
for help on using the repository browser.