Changeset 4335


Ignore:
Timestamp:
06/30/10 10:08:42 (15 years ago)
Author:
seroussi
Message:

save results only if needed for serial solutions

Location:
issm/trunk/src/m/solutions
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/m/solutions/balancedthickness2_core.m

    r4222 r4335  
    88        verbose=femmodel.parameters.Verbose;
    99        dim=femmodel.parameters.Dim;
     10        solution_type=femmodel.parameters.SolutionType;
    1011
    1112        %Activate formulation
    … …  
    2526        %femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
    2627
    27         displaystring(verbose,'\n%s',['saving results...']);
    28         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     28        if solution_type==Balancedthickness2SolutionEnum,
     29                displaystring(verbose,'\n%s',['saving results...']);
     30                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     31        end
    2932       
    3033end %end function

  • issm/trunk/src/m/solutions/balancedthickness_core.m

    r4222 r4335  
    88        verbose=femmodel.parameters.Verbose;
    99        dim=femmodel.parameters.Dim;
     10        solution_type=femmodel.parameters.SolutionType;
    1011
    1112        %Activate formulation
    … …  
    2526        femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
    2627
    27         displaystring(verbose,'\n%s',['saving results...']);
    28         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     28        if solution_type==BalancethicknessSolutionEnum,
     29                displaystring(verbose,'\n%s',['saving results...']);
     30                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     31        end
    2932       
    3033end %end function

  • issm/trunk/src/m/solutions/balancedvelocities_core.m

    r4222 r4335  
    88        verbose=femmodel.parameters.Verbose;
    99        dim=femmodel.parameters.Dim;
     10        solution_type=femmodel.parameters.SolutionType;
    1011
    1112        %Activate formulation
    … …  
    2627        femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
    2728
    28         displaystring(verbose,'\n%s',['saving results...']);
    29         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
    30         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
     29        if solution_type==BalancedvelocitiesSolutionEnum,
     30                displaystring(verbose,'\n%s',['saving results...']);
     31                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
     32                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
     33        end
    3134
    3235end %end function

  • issm/trunk/src/m/solutions/bedslope_core.m

    r4302 r4335  
    99        verbose=femmodel.parameters.Verbose;
    1010        dim=femmodel.parameters.Dim;
     11        solution_type=femmodel.parameters.SolutionType;
    1112
    1213        displaystring(verbose,'\n%s',['computing bed slope...']);
    … …  
    2526        }
    2627       
    27         displaystring(verbose,'\n%s',['saving results...']);
    28         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeXEnum);
    29         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeYEnum);
     28        if solution_type==BesSlopeSolutionEnum,
     29                displaystring(verbose,'\n%s',['saving results...']);
     30                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeXEnum);
     31                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeYEnum);
     32        end

  • issm/trunk/src/m/solutions/diagnostic_core.m

    r4254 r4335  
    1818        if(isstokes) stokesreconditioning=femmodel.parameters.stokesreconditioning; end
    1919        qmu_analysis=femmodel.parameters.QmuAnalysis;
     20        solution_type=femmodel.parameters.SolutionType;
    2021
    2122        %for qmu analysis, be sure the velocity input we are starting from  is the one in the parameters:
    … …  
    7071        end
    7172
    72         displaystring(verbose,'\n%s',['saving results...']);
    73         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
    74         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
    75         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);
    76         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
    77         if(dim==3) femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum); end
     73        if solution_type==DiagnosticSolutionEnum,
     74                displaystring(verbose,'\n%s',['saving results...']);
     75                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
     76                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
     77                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);
     78                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
     79                if(dim==3) femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum); end
     80        end

  • issm/trunk/src/m/solutions/prognostic2_core.m

    r4222 r4335  
    77        %recover parameters common to all solutions
    88        verbose=femmodel.parameters.Verbose;
     9        solution_type=femmodel.parameters.SolutionType;
    910
    1011        %Activate formulation
    … …  
    2122%       femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
    2223
    23         displaystring(verbose,'\n%s',['saving results...']);
    24         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     24        if solution_type==Prognostic2SolutionEnum,
     25                displaystring(verbose,'\n%s',['saving results...']);
     26                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     27        end
    2528       
    2629end %end function

  • issm/trunk/src/m/solutions/prognostic_core.m

    r4222 r4335  
    77        %recover parameters common to all solutions
    88        verbose=femmodel.parameters.Verbose;
     9        solution_type=femmodel.parameters.SolutionType;
    910
    1011        %Activate formulation
    … …  
    2122        femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
    2223
    23         displaystring(verbose,'\n%s',['saving results...']);
    24         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     24        if solution_type==PrognosticSolutionEnum,
     25                displaystring(verbose,'\n%s',['saving results...']);
     26                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
     27        end
    2528
    2629end %end function

  • issm/trunk/src/m/solutions/steadystate_core.m

    r4222 r4335  
    99        verbose=femmodel.parameters.Verbose;
    1010        dim=femmodel.parameters.Dim;
     11        solution_type=femmodel.parameters.SolutionType;
    1112
    1213        %Initialize counter
    … …  
    4344        end
    4445
    45         displaystring(verbose,'\n%s',['saving results...']);
    46         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
    47         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
    48         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);
    49         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
    50         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum);
     46        if solution_type==SteadystateSolutionEnum,
     47                displaystring(verbose,'\n%s',['saving results...']);
     48                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
     49                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
     50                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);
     51                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
     52                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum);
    5153        end
     54end
    5255
    5356end %end of function

  • issm/trunk/src/m/solutions/surfaceslope_core.m

    r4302 r4335  
    99        verbose=femmodel.parameters.Verbose;
    1010        dim=femmodel.parameters.Dim;
     11        solution_type=femmodel.parameters.SolutionType;
    1112
    1213        displaystring(verbose,'\n%s',['computing surface slope...']);
    … …  
    2526        }
    2627       
    27         displaystring(verbose,'\n%s',['saving results...']);
    28         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeXEnum);
    29         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeYEnum);
     28        if solution_type==SurfaceSlopeSolutionEnum,
     29                displaystring(verbose,'\n%s',['saving results...']);
     30                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeXEnum);
     31                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeYEnum);
     32        end
    3033
    3134end %end function

  • issm/trunk/src/m/solutions/thermal_core.m

    r4222 r4335  
    1010        ndt=femmodel.parameters.Ndt;
    1111        dt=femmodel.parameters.Dt;
     12        solution_type=femmodel.parameters.SolutionType;
    1213
    1314        %Compute number of timesteps
    … …  
    2728                femmodel=thermal_core_step(femmodel,i,time);
    2829
    29                 displaystring(verbose,'\n%s',['saving results...']);
    30                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,i,time);
    31                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingRateEnum,i,time);
     30                if solution_type==ThermalSolutionEnum,
     31                        displaystring(verbose,'\n%s',['saving results...']);
     32                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,i,time);
     33                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingRateEnum,i,time);
     34                end
    3235        end
    3336

  • issm/trunk/src/m/solutions/transient2d_core.m

    r4325 r4335  
    99        ndt=femmodel.parameters.Ndt;
    1010        dt=femmodel.parameters.Dt;
     11        solution_type=femmodel.parameters.SolutionType;
    1112
    1213        %Initialize
    … …  
    2930                [femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=UpdateGeometry(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
    3031
    31                 displaystring(verbose,'\n%s',['saving results...']);
    32                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
    33                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
    34                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
    35                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
    36                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
    37                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
     32                if solution_type==Transient2DSolutionEnum,
     33                        displaystring(verbose,'\n%s',['saving results...']);
     34                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
     35                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
     36                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
     37                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
     38                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
     39                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
     40                end
    3841        end
    39 
    4042end %end of function

  • issm/trunk/src/m/solutions/transient3d_core.m

    r4222 r4335  
    99        ndt=femmodel.parameters.Ndt;
    1010        dt=femmodel.parameters.Dt;
     11        solution_type=femmodel.parameters.SolutionType;
    1112
    1213        %Initialize
    … …  
    3839                [femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=UpdateVertexPositions(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
    3940
    40                 displaystring(verbose,'\n%s',['saving results...']);
    41                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
    42                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
    43                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,step,time);
    44                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
    45                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
    46                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
    47                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
    48                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,step,time);
    49                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingEnum,step,time);
     41                if solution_type==Transient3DSolutionEnum,
     42                        displaystring(verbose,'\n%s',['saving results...']);
     43                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
     44                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
     45                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,step,time);
     46                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
     47                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
     48                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
     49                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
     50                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,step,time);
     51                        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingEnum,step,time);
     52                end
    5053        end
    5154end %end of function
Note: See TracChangeset for help on using the changeset viewer.