Changeset 13900


Ignore:
Timestamp:
11/07/12 16:04:40 (12 years ago)
Author:
jschierm
Message:

NEW: Working (up to ADOLC) test3???.py (plus cosmetic Matlab changes).

Location:
issm/trunk-jpl/test/NightlyRun
Files:
3 added
13 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/test/NightlyRun/test3001.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',50000)
     10md=triangle(model(),'../Exp/Square.exp',50000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1514md.cluster=generic('name',oshostname(),'np',3)
    1615md.diagnostic.requested_outputs=StressTensorEnum()
    17 md.autodiff.isautodiff=true
     16md.autodiff.isautodiff=True
    1817md=solve(md,DiagnosticSolutionEnum())
    1918
    20 
    21 # Fields and tolerances to track changes
    22 
     19#Fields and tolerances to track changes
    2320field_names     =['Vx','Vy','Vel','Pressure',\
    2421        'StressTensorxx','StressTensoryy','StressTensorxy']

  • issm/trunk-jpl/test/NightlyRun/test3002.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',180000)
     10md=triangle(model(),'../Exp/Square.exp',180000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1514md=setflowequation(md,'macayeal','all')
    1615md.cluster=generic('name',oshostname(),'np',3)
    17 md.autodiff.isautodiff=true
     16md.autodiff.isautodiff=True
    1817md=solve(md,DiagnosticSolutionEnum())
    1918
    20 
    21 # Fields and tolerances to track changes
    22 
     19#Fields and tolerances to track changes
    2320field_names     =['Vx','Vy','Vz','Vel','Pressure']
    2421field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13]

  • issm/trunk-jpl/test/NightlyRun/test3003.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',180000)
     10md=triangle(model(),'../Exp/Square.exp',180000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1615md.cluster=generic('name',oshostname(),'np',3)
    1716md.diagnostic.requested_outputs=StressTensorEnum()
    18 md.autodiff.isautodiff=true
     17md.autodiff.isautodiff=True
    1918md=solve(md,DiagnosticSolutionEnum())
    2019
    21 
    22 # Fields and tolerances to track changes
    23 
     20#Fields and tolerances to track changes
    2421field_names     =['Vx','Vy','Vz','Vel','Pressure',\
    2522        'StressTensorxx','StressTensoryy','StressTensorzz','StressTensorxy','StressTensorxz','StressTensoryz']

  • issm/trunk-jpl/test/NightlyRun/test3004.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',180000)
     10md=triangle(model(),'../Exp/Square.exp',180000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1514md=setflowequation(md,'stokes','all')
    1615md.cluster=generic('name',oshostname(),'np',3)
    17 md.autodiff.isautodiff=true
     16md.autodiff.isautodiff=True
    1817md=solve(md,DiagnosticSolutionEnum())
    1918
    20 
    21 # Fields and tolerances to track changes
    22 
     19#Fields and tolerances to track changes
    2320field_names     =['Vx','Vy','Vz','Vel','Pressure']
    2421field_tolerances=[1e-08,1e-08,1e-07,1e-08,1e-08]

  • issm/trunk-jpl/test/NightlyRun/test3005.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',150000)
     10md=triangle(model(),'../Exp/Square.exp',150000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    1413md=setflowequation(md,'macayeal','all')
    1514md.cluster=generic('name',oshostname(),'np',3)
    16 md.autodiff.isautodiff=true
     15md.autodiff.isautodiff=True
    1716md=solve(md,PrognosticSolutionEnum())
    1817
    19 
    20 # Fields and tolerances to track changes
    21 
     18#Fields and tolerances to track changes
    2219field_names     =['Thickness']
    2320field_tolerances=[1e-13]

  • issm/trunk-jpl/test/NightlyRun/test3006.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from meshconvert import *
    … …  
    85from parameterize import *
    96from setflowequation import *
     7from EnumDefinitions import *
    108from solve import *
     9from MatlabFuncs import *
    1110
    12 md=triangle(model(),'../Exp/Square.exp',150000)
     11md=triangle(model(),'../Exp/Square.exp',150000.)
    1312md=meshconvert(md)
    1413md=setmask(md,'all','')
    … …  
    1817md.prognostic.stabilization=3
    1918md.prognostic.spcthickness=md.geometry.thickness
    20 md.autodiff.isautodiff=true
     19md.autodiff.isautodiff=True
    2120md=solve(md,PrognosticSolutionEnum())
    2221
    23 
    24 # Fields and tolerances to track changes
    25 
     22#Fields and tolerances to track changes
    2623field_names     =['Thickness']
    2724field_tolerances=[1e-13]

  • issm/trunk-jpl/test/NightlyRun/test3007.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',150000)
     10md=triangle(model(),'../Exp/Square.exp',150000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1514md.extrude(5,3.)
    1615md.cluster=generic('name',oshostname(),'np',3)
    17 md.autodiff.isautodiff=true
     16md.autodiff.isautodiff=True
    1817md=solve(md,PrognosticSolutionEnum())
    1918
    20 
    21 # Fields and tolerances to track changes
    22 
     19#Fields and tolerances to track changes
    2320field_names     =['Thickness']
    2421field_tolerances=[1e-13]

  • issm/trunk-jpl/test/NightlyRun/test3008.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',180000)
     10md=triangle(model(),'../Exp/Square.exp',180000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1615md.timestepping.time_step=0
    1716md.cluster=generic('name',oshostname(),'np',3)
    18 md.autodiff.isautodiff=true
     17md.autodiff.isautodiff=True
    1918md=solve(md,ThermalSolutionEnum())
    2019
    21 
    22 # Fields and tolerances to track changes
    23 
     20#Fields and tolerances to track changes
    2421field_names     =['Temperature','BasalforcingsMeltingRate']
    2522field_tolerances=[1e-13,1e-13]

  • issm/trunk-jpl/test/NightlyRun/test3009.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',180000)
     10md=triangle(model(),'../Exp/Square.exp',180000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1918md.transient.isthermal=1
    2019md.transient.isgroundingline=0
    21 md.autodiff.isautodiff=true
     20md.autodiff.isautodiff=True
    2221md=solve(md,TransientSolutionEnum())
    2322
    24 
    25 # Fields and tolerances to track changes
    26 
     23#Fields and tolerances to track changes
    2724field_names     =['Temperature','BasalforcingsMeltingRate']
    2825field_tolerances=[1e-13,1e-13]

  • issm/trunk-jpl/test/NightlyRun/test3010.py

    r13879 r13900  
    1 from MatlabFuncs import *
    21from model import *
    3 from EnumDefinitions import *
    4 import numpy
    52from triangle import *
    63from setmask import *
    74from parameterize import *
    85from setflowequation import *
     6from EnumDefinitions import *
    97from solve import *
     8from MatlabFuncs import *
    109
    11 md=triangle(model(),'../Exp/Square.exp',150000)
     10md=triangle(model(),'../Exp/Square.exp',150000.)
    1211md=setmask(md,'all','')
    1312md=parameterize(md,'../Par/SquareShelfConstrained.py')
    … …  
    1615md.transient.requested_outputs=IceVolumeEnum()
    1716
    18 
    19 md.autodiff.isautodiff=true
     17md.autodiff.isautodiff=True
    2018md=solve(md,TransientSolutionEnum())
    2119
    22 
    23 # Fields and tolerances to track changes
    24 
     20#Fields and tolerances to track changes
    2521field_names     =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Volume1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Volume2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Volume3']
    2622field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\

  • issm/trunk-jpl/test/NightlyRun/test3015.m

    r13670 r13900  
    2222md.autodiff.driver='fos_forward';
    2323
    24 %parameters for the step-wise devivative
    25 delta=.001;
     24%parameters for the step-wise derivative
     25delta=0.001;
    2626h1=md.geometry.thickness(index);
    27 h0=h1*(1-delta);
    28 h2=h1*(1+delta);
     27h0=h1*(1.-delta);
     28h2=h1*(1.+delta);
    2929deltaH=(h2-h0);
    3030

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

    r13670 r13900  
    1919jac_reverse=md.results.TransientSolution(1).AutodiffJacobian;
    2020
    21 
    2221%now run vectorial forward mode
    2322md.autodiff.independents={independent('name','Thickness','type','vertex','nods',md.mesh.numberofvertices,'fov_forward_indices',(1:md.mesh.numberofvertices)')};
    … …  
    3029jac_forward=md.results.TransientSolution(1).AutodiffJacobian;
    3130
    32 
    3331%Fields and tolerances to track changes
    3432field_names     ={'Jac Forward','Jac Reverse','Jac Forward - Reverse'};

  • issm/trunk-jpl/test/NightlyRun/test3020.m

    r13670 r13900  
    88md=setflowequation(md,'macayeal','all');
    99md.cluster=generic('name',oshostname(),'np',3);
    10 md.transient.requested_outputs=[IceVolumeEnum();MaxVelEnum];
     10md.transient.requested_outputs=[IceVolumeEnum();MaxVelEnum()];
    1111md.verbose=verbose('autodiff',true);
    12 md.diagnostic.restol=.000001;
     12md.diagnostic.restol=0.000001;
    1313
    1414%setup autodiff parameters
    … …  
    2424md.autodiff.driver='fos_forward';
    2525
    26 %parameters for the step-wise devivative
    27 delta=.00001;
     26%parameters for the step-wise derivative
     27delta=0.00001;
    2828h1=md.geometry.thickness(index);
    29 h0=h1*(1-delta);
    30 h2=h1*(1+delta);
     29h0=h1*(1.-delta);
     30h2=h1*(1.+delta);
    3131deltaH=(h2-h0);
    3232
    … …  
    6363dMaxVdh_an=(MaxV2-MaxV0)/deltaH;
    6464
    65 
    6665%evaluate derivative using ADOLC
    6766md=md2;
    … …  
    8079disp(sprintf('dMaxV/dh: analytical:  %16.16g\n       using adolc:  %16.16g\n',dMaxVdh_an,dMaxVdh_ad));
    8180
    82 
    8381%Fields and tolerances to track changes
    8482field_names     ={'dV/dh-dV/dh0','dMaxV/dh-dMaxV/dh0'};
Note: See TracChangeset for help on using the changeset viewer.