Changeset 13518


Ignore:
Timestamp:
10/03/12 09:36:29 (12 years ago)
Author:
jschierm
Message:

NEW: Working test212.py.

Location:
issm/trunk-jpl/test/NightlyRun
Files:
1 edited
1 moved

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test212.m

    r12791 r13518  
    2525field_values={...
    2626        (md.results.DiagnosticSolution.Gradient1),...
    27         md.results.DiagnosticSolution.J,...
     27        (md.results.DiagnosticSolution.J),...
    2828        (md.results.DiagnosticSolution.MaterialsRheologyBbar),...
    2929        (md.results.DiagnosticSolution.Pressure),...
  • issm/trunk-jpl/test/NightlyRun/test212.py

    r13497 r13518  
    3131md.inversion.iscontrol=1
    3232md.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)
     33md.inversion.min_parameters=10**6*ones((md.mesh.numberofvertices,1))
     34md.inversion.max_parameters=2*10**9*ones((md.mesh.numberofvertices,1))
    3535md.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)
     36md.inversion.cost_functions=101*ones((md.inversion.nsteps,1))
     37md.inversion.cost_functions_coefficients=ones((md.mesh.numberofvertices,1))
     38md.inversion.gradient_scaling=10**7*ones((md.inversion.nsteps,1))
     39md.inversion.maxiter_per_step=2*ones((md.inversion.nsteps,1))
     40md.inversion.step_threshold=0.3*ones((md.inversion.nsteps,1))
    4141md.inversion.vx_obs=md.initialization.vx
    4242md.inversion.vy_obs=md.initialization.vy
     
    4949# Fields and tolerances to track changes
    5050
    51 field_names     =['Gradient' 'Misfits' 'MaterialsRheologyBbar' 'Pressure' 'Vel' 'Vx' 'Vy']
     51field_names     =['Gradient','Misfits','MaterialsRheologyBbar','Pressure','Vel','Vx','Vy']
    5252field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
    5353field_values=[\
    5454md.results['DiagnosticSolution'][1]['Gradient1'],\
    55 md.results['DiagnosticSolution.J'],\
     55md.results['DiagnosticSolution'][1]['J'],\
    5656md.results['DiagnosticSolution'][1]['MaterialsRheologyBbar'],\
    5757md.results['DiagnosticSolution'][1]['Pressure'],\
Note: See TracChangeset for help on using the changeset viewer.