balancevelocity_core.cpp

Go to the documentation of this file.

 
#include "./cores.h"
#include "../toolkits/toolkits.h"
#include "../classes/classes.h"
#include "../shared/shared.h"
#include "../modules/modules.h"
#include "../solutionsequences/solutionsequences.h"
 
void balancevelocity_core(FemModel* femmodel){
 
   /*parameters: */
   bool        save_results;
   IssmDouble  l = 8.;
 
   /*recover parameters: */
   femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
 
   if(VerboseSolution()) _printf0_("computing smooth driving stress:\n");
   femmodel->parameters->SetParam(l,SmoothThicknessMultiplierEnum);
   femmodel->SetCurrentConfiguration(SmoothAnalysisEnum);
   femmodel->parameters->SetParam(DrivingStressXEnum,InputToSmoothEnum);
   solutionsequence_linear(femmodel);
   femmodel->parameters->SetParam(DrivingStressYEnum,InputToSmoothEnum);
   solutionsequence_linear(femmodel);
 
   if(VerboseSolution()) _printf0_("   call computational core:\n");
   femmodel->SetCurrentConfiguration(BalancevelocityAnalysisEnum);
   solutionsequence_linear(femmodel);
 
   if(save_results){
      if(VerboseSolution()) _printf0_("   saving results\n");
      int outputs[3] = {DrivingStressXEnum,DrivingStressYEnum,VelEnum};
      femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3);
   }
 
}