 |
Ice Sheet System Model
4.18
Code documentation
|
Go to the documentation of this file.
7 #include "../../shared/shared.h"
8 #include "../../toolkits/toolkits.h"
9 #include "../../classes/Inputs2/DatasetInput2.h"
34 int domaintype,numcomponents;
54 default:
_error_(
"not supported yet");
68 Input2* vyobs_input = NULL;
76 for(
int ig=gauss->
begin();ig<gauss->end();ig++){
100 if(numcomponents==2){
101 scalex=pow(meanvel/(vxobs+epsvel),2);
if(vxobs==0)scalex=0;
102 scaley=pow(meanvel/(vyobs+epsvel),2);
if(vyobs==0)scaley=0;
103 misfit=0.5*(scalex*pow((vx-vxobs),2)+scaley*pow((vy-vyobs),2));
106 scalex=pow(meanvel/(vxobs+epsvel),2);
if(vxobs==0)scalex=0;
107 misfit=0.5*(scalex*pow((vx-vxobs),2));
111 Jelem+=misfit*weight*Jdet*gauss->
weight;
116 xDelete<IssmDouble>(xyz_list);
Declaration of Vertices class.
Declaration of Nodes class.
void FindParam(bool *pvalue, int paramenum)
@ InversionCostFunctionsCoefficientsEnum
Declaration of Parameters class.
static ISSM_MPI_Comm GetComm(void)
Declaration of Elements class.
@ SurfaceRelVelMisfitEnum
virtual Input2 * GetInput2(int inputenum)=0
void DeleteMaterials(void)
virtual DatasetInput2 * GetDatasetInput2(int inputenum)
header file for inverse methods misfit computation
Declaration of Materials class.
virtual Gauss * NewGauss(void)=0
void SurfaceRelVelMisfitx(IssmDouble *pJ, Elements *elements, Nodes *nodes, Vertices *vertices, Loads *loads, Materials *materials, Parameters *parameters)
void GetVerticesCoordinates(IssmDouble **xyz_list)
int ISSM_MPI_Bcast(void *buffer, int count, ISSM_MPI_Datatype datatype, int root, ISSM_MPI_Comm comm)
IssmDouble SurfaceRelVelMisfit(Element *element)
Declaration of Loads class.
#define _error_(StreamArgs)
virtual Element * SpawnTopElement(void)=0
virtual int begin(void)=0
Object * GetObjectByOffset(int offset)
virtual void JacobianDeterminant(IssmDouble *Jdet, IssmDouble *xyz_list, Gauss *gauss)=0
virtual void GaussPoint(int ig)=0
int ISSM_MPI_Reduce(void *sendbuf, void *recvbuf, int count, ISSM_MPI_Datatype datatype, ISSM_MPI_Op op, int root, ISSM_MPI_Comm comm)