Cfsurfacelogvel Class Reference

#include <Cfsurfacelogvel.h>

Inheritance diagram for Cfsurfacelogvel:

Public Member Functions
	Cfsurfacelogvel ()
	Cfsurfacelogvel (char *in_name, int in_definitionenum, IssmDouble in_datatime, bool timepassedflag)
	~Cfsurfacelogvel ()
Object *	copy ()
void	DeepEcho (void)
void	Echo (void)
int	Id (void)
void	Marshall (char **pmarshalled_data, int *pmarshalled_data_size, int marshall_direction)
int	ObjectEnum (void)
int	DefinitionEnum ()
char *	Name ()
IssmDouble	Response (FemModel *femmodel)
IssmDouble	Cfsurfacelogvel_Calculation (Element *element, int definitionenum)
Public Member Functions inherited from Object
virtual	~Object ()
Public Member Functions inherited from Definition
virtual	~Definition ()

Data Fields
int	definitionenum
char *	name
IssmDouble	datatime
bool	timepassedflag
IssmDouble	last_time
int	lock
IssmDouble	misfit

Detailed Description

Definition at line 15 of file Cfsurfacelogvel.h.

Constructor & Destructor Documentation

Cfsurfacelogvel() [1/2]

Cfsurfacelogvel::Cfsurfacelogvel

(

)

Definition at line 27 of file Cfsurfacelogvel.cpp.

                                {/*{{{*/
 
   this->definitionenum = -1;
   this->name = NULL;
   this->misfit=0;
   this->lock=0;
   this->datatime=0.;
   this->timepassedflag = false;
   this->last_time=0.;
 
}

Cfsurfacelogvel() [2/2]

Cfsurfacelogvel::Cfsurfacelogvel	(	char *	in_name,
		int	in_definitionenum,
		IssmDouble	in_datatime,
		bool	timepassedflag
	)

Definition at line 39 of file Cfsurfacelogvel.cpp.

                                                                                                                    {/*{{{*/
 
   this->definitionenum=in_definitionenum;
 
   this->name     = xNew<char>(strlen(in_name)+1);
   xMemCpy<char>(this->name,in_name,strlen(in_name)+1);
 
   this->datatime=in_datatime;
   this->timepassedflag=in_timepassedflag;
 
   this->misfit=0;
   this->lock=0;
}

~Cfsurfacelogvel()

Cfsurfacelogvel::~Cfsurfacelogvel

(

)

Definition at line 53 of file Cfsurfacelogvel.cpp.

                                 {/*{{{*/
   if(this->name)xDelete(this->name);
   this->misfit=0;
   this->lock=0;
}

Member Function Documentation

copy()

Object * Cfsurfacelogvel::copy ( void  )

virtual

Implements Object.

Definition at line 60 of file Cfsurfacelogvel.cpp.

                              {/*{{{*/
   Cfsurfacelogvel* mf = new Cfsurfacelogvel(this->name,this->definitionenum,this->datatime,this->timepassedflag);
   mf->misfit=this->misfit;
   mf->lock=this->lock;
   return (Object*) mf;
}

DeepEcho()

void Cfsurfacelogvel::DeepEcho ( void  )

virtual

Implements Object.

Definition at line 67 of file Cfsurfacelogvel.cpp.

                                  {/*{{{*/
   this->Echo();
}

Echo()

void Cfsurfacelogvel::Echo ( void  )

virtual

Implements Object.

Definition at line 71 of file Cfsurfacelogvel.cpp.

                              {/*{{{*/
   _printf_(" Cfsurfacelogvel: " << name << " " << this->definitionenum << "\n");
   _printf_("    datatime: " << datatime << "\n");
   _printf_("    timepassedflag: "<<timepassedflag<<"\n");
}

Id()

int Cfsurfacelogvel::Id ( void  )

virtual

Implements Object.

Definition at line 77 of file Cfsurfacelogvel.cpp.

                           {/*{{{*/
   return -1;
}

Marshall()

void Cfsurfacelogvel::Marshall ( char **  pmarshalled_data, int *  pmarshalled_data_size, int  marshall_direction  )

virtual

Implements Object.

Definition at line 81 of file Cfsurfacelogvel.cpp.

                                                                                                        {/*{{{*/
   _error_("not implemented yet!"); 
} 

ObjectEnum()

int Cfsurfacelogvel::ObjectEnum ( void  )

virtual

Implements Object.

Definition at line 85 of file Cfsurfacelogvel.cpp.

                                   {/*{{{*/
   return CfsurfacelogvelEnum;
}

DefinitionEnum()

int Cfsurfacelogvel::DefinitionEnum ( )

virtual

Implements Definition.

Definition at line 90 of file Cfsurfacelogvel.cpp.

                                   {/*{{{*/
   return this->definitionenum;
}

Name()

char * Cfsurfacelogvel::Name ( )

virtual

Implements Definition.

Definition at line 94 of file Cfsurfacelogvel.cpp.

                           {/*{{{*/
   char* name2=xNew<char>(strlen(this->name)+1);
   xMemCpy(name2,this->name,strlen(this->name)+1);
 
   return name2;
}

Response()

IssmDouble Cfsurfacelogvel::Response ( FemModel *  femmodel )

virtual

Implements Definition.

Definition at line 101 of file Cfsurfacelogvel.cpp.

                                                      {/*{{{*/
 
    /*diverse: */
    IssmDouble time;
 
    /*recover time parameters: */
    femmodel->parameters->FindParam(&time,TimeEnum);
    if(time < last_time) timepassedflag = false;
    last_time = time;
 
       int i;
       IssmDouble J=0.;
       IssmDouble J_sum=0.;
 
    if(datatime<=time && !timepassedflag){
       for(i=0;i<femmodel->elements->Size();i++){
          Element* element=(Element*)femmodel->elements->GetObjectByOffset(i);
          J+=this->Cfsurfacelogvel_Calculation(element,definitionenum);
       }
 
       ISSM_MPI_Allreduce ( (void*)&J,(void*)&J_sum,1,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,IssmComm::GetComm());
       ISSM_MPI_Bcast(&J_sum,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm());
       J=J_sum;
 
       timepassedflag = true;
       return J;
      }
    else return J;
 }

Cfsurfacelogvel_Calculation()

IssmDouble Cfsurfacelogvel::Cfsurfacelogvel_Calculation	(	Element *	element,
		int	definitionenum
	)

Definition at line 131 of file Cfsurfacelogvel.cpp.

                                                                                           {/*{{{*/
 
   int        domaintype,numcomponents;
   IssmDouble Jelem=0.;
   IssmDouble epsvel=2.220446049250313e-16;
   IssmDouble meanvel=3.170979198376458e-05; /*1000 m/yr*/
   IssmDouble velocity_mag,obs_velocity_mag;
   IssmDouble misfit,Jdet;
   IssmDouble vx,vy,vxobs,vyobs,weight;
   IssmDouble* xyz_list = NULL;
 
   /*Get basal element*/
   if(!element->IsOnSurface()) return 0.;
 
   /*If on water, return 0: */
   if(!element->IsIceInElement()) return 0.;
 
   /*Get problem dimension*/
   element->FindParam(&domaintype,DomainTypeEnum);
   switch(domaintype){
      case Domain2DverticalEnum:   numcomponents   = 1; break;
      case Domain3DEnum:           numcomponents   = 2; break;
      case Domain2DhorizontalEnum: numcomponents   = 2; break;
      default: _error_("not supported yet");
   }
 
   /*Spawn surface element*/
   Element* topelement = element->SpawnTopElement();
 
   /* Get node coordinates*/
   topelement->GetVerticesCoordinates(&xyz_list);
 
   /*Get model values*/
   Input2 *vx_input = topelement->GetInput2(VxEnum); _assert_(vx_input);
   Input2 *vy_input = NULL;
   if(numcomponents==2){
         vy_input = topelement->GetInput2(VyEnum); _assert_(vy_input);
   }
 
   /*Retrieve all inputs we will be needing: */
   DatasetInput2 *datasetinput = topelement->GetDatasetInput2(definitionenum); _assert_(datasetinput);
 
   /* Start  looping on the number of gaussian points: */
   Gauss* gauss=topelement->NewGauss(2);
   for(int ig=gauss->begin();ig<gauss->end();ig++){
 
      gauss->GaussPoint(ig);
 
      /* Get Jacobian determinant: */
      topelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
 
      /*Get all parameters at gaussian point*/
      datasetinput->GetInputValue(&weight,gauss,WeightsSurfaceObservationEnum);
      vx_input->GetInputValue(&vx,gauss);
      datasetinput->GetInputValue(&vxobs,gauss,VxObsEnum);
      if(numcomponents==2){
         vy_input->GetInputValue(&vy,gauss);
         datasetinput->GetInputValue(&vyobs,gauss,VyObsEnum);
      }
 
      /*Compute SurfaceLogVelMisfit:
       *        *                 [        vel + eps     ] 2
       *               * J = 4 \bar{v}^2 | log ( -----------  ) |
       *                      *                 [       vel   + eps    ]
       *                             *                            obs
       *                                    */
      if(numcomponents==1){
         velocity_mag    =fabs(vx)+epsvel;
         obs_velocity_mag=fabs(vxobs)+epsvel;
      }
      else{
         velocity_mag    =sqrt(vx*vx+vy*vy)+epsvel;
         obs_velocity_mag=sqrt(vxobs*vxobs+vyobs*vyobs)+epsvel;
      }
 
      misfit=4*pow(meanvel,2)*pow(log(velocity_mag/obs_velocity_mag),2);
 
      /*Add to cost function*/
      Jelem+=misfit*weight*Jdet*gauss->weight;
 
   }
 
   /*clean up and Return: */
   if(domaintype!=Domain2DhorizontalEnum){topelement->DeleteMaterials(); delete topelement;};
   xDelete<IssmDouble>(xyz_list);
   delete gauss;
   return Jelem;
}/*}}}*/

Field Documentation

definitionenum

int Cfsurfacelogvel::definitionenum

Definition at line 19 of file Cfsurfacelogvel.h.

name

char* Cfsurfacelogvel::name

Definition at line 20 of file Cfsurfacelogvel.h.

datatime

IssmDouble Cfsurfacelogvel::datatime

Definition at line 21 of file Cfsurfacelogvel.h.

timepassedflag

bool Cfsurfacelogvel::timepassedflag

Definition at line 22 of file Cfsurfacelogvel.h.

last_time

IssmDouble Cfsurfacelogvel::last_time

Definition at line 23 of file Cfsurfacelogvel.h.

lock

int Cfsurfacelogvel::lock

Definition at line 25 of file Cfsurfacelogvel.h.

misfit

IssmDouble Cfsurfacelogvel::misfit

Definition at line 26 of file Cfsurfacelogvel.h.

---

The documentation for this class was generated from the following files:

• src/c/classes/Cfsurfacelogvel.h
• src/c/classes/Cfsurfacelogvel.cpp