Ice Sheet System Model  4.18
Code documentation
Public Member Functions | Data Fields | Private Attributes
GaussTetra Class Reference

#include <GaussTetra.h>

Inheritance diagram for GaussTetra:
Gauss

Public Member Functions

 GaussTetra ()
 
 GaussTetra (int order)
 
 GaussTetra (int index1, int index2, int index3, int order)
 
 ~GaussTetra ()
 
int begin (void)
 
void Echo (void)
 
int end (void)
 
int Enum (void)
 
void GaussPoint (int ig)
 
void GaussNode (int finitelement, int iv)
 
void GaussVertex (int iv)
 
void SynchronizeGaussBase (Gauss *gauss)
 
- Public Member Functions inherited from Gauss
virtual ~Gauss ()
 

Data Fields

IssmDouble coord1
 
IssmDouble coord2
 
IssmDouble coord3
 
IssmDouble coord4
 
- Data Fields inherited from Gauss
IssmDouble weight
 

Private Attributes

int numgauss
 
IssmDouble * weights
 
IssmDouble * coords1
 
IssmDouble * coords2
 
IssmDouble * coords3
 
IssmDouble * coords4
 

Detailed Description

Definition at line 12 of file GaussTetra.h.

Constructor & Destructor Documentation

◆ GaussTetra() [1/3]

GaussTetra::GaussTetra ( )

Definition at line 15 of file GaussTetra.cpp.

15  {/*{{{*/
16 
17  numgauss=-1;
18 
19  weights=NULL;
20  coords1=NULL;
21  coords2=NULL;
22  coords3=NULL;
23  coords4=NULL;
24 
25  weight=UNDEF;
26  coord1=UNDEF;
27  coord2=UNDEF;
28  coord3=UNDEF;
29  coord4=UNDEF;
30 }

◆ GaussTetra() [2/3]

GaussTetra::GaussTetra ( int  order)

Definition at line 32 of file GaussTetra.cpp.

32  {/*{{{*/
33 
34  /*Get gauss points*/
35  GaussLegendreTetra(&numgauss,&coords1,&coords2,&coords3,&coords4,&weights,order);
36 
37  /*Rescale weights if necessary*/
38  IssmDouble sumweights = 0.;
39  for(int i=0;i<numgauss;i++) sumweights += this->weights[i];
40  if(sumweights==1.){
41  for(int i=0;i<numgauss;i++) this->weights[i] = this->weights[i]/6.; /*rescale volume to 1/6*/
42  }
43 
44  /*Check final weights in debugging mode*/
45  #ifdef _ISSM_DEBUG_
46  sumweights = 0.; for(int i=0;i<numgauss;i++) sumweights += this->weights[i];
47  _assert_(sumweights>1./6.-1e-10);
48  _assert_(sumweights<1./6.+1e-10);
49  #endif
50 
51  /*Initialize static fields as undefinite*/
52  weight=UNDEF;
53  coord1=UNDEF;
54  coord2=UNDEF;
55  coord3=UNDEF;
56 }

◆ GaussTetra() [3/3]

GaussTetra::GaussTetra ( int  index1,
int  index2,
int  index3,
int  order 
)

Definition at line 58 of file GaussTetra.cpp.

58  {/*{{{*/
59 
60  /*Basal Tria*/
61  if(index1==0 && index2==1 && index3==2){
62  GaussLegendreTria(&numgauss,&coords1,&coords2,&coords3,&weights,order);
63  coords4=xNew<IssmDouble>(numgauss);
64  for(int i=0;i<numgauss;i++) coords4[i]=0.;
65  }
66  else if(index1==0 && index2==3 && index3==1){
67  GaussLegendreTria(&numgauss,&coords1,&coords2,&coords4,&weights,order);
68  coords3=xNew<IssmDouble>(numgauss);
69  for(int i=0;i<numgauss;i++) coords3[i]=0.;
70  }
71  else if(index1==1 && index2==3 && index3==2){
72  GaussLegendreTria(&numgauss,&coords2,&coords3,&coords4,&weights,order);
73  coords1=xNew<IssmDouble>(numgauss);
74  for(int i=0;i<numgauss;i++) coords1[i]=0.;
75  }
76  else if(index1==0 && index2==2 && index3==3){
77  GaussLegendreTria(&numgauss,&coords1,&coords3,&coords4,&weights,order);
78  coords2=xNew<IssmDouble>(numgauss);
79  for(int i=0;i<numgauss;i++) coords2[i]=0.;
80  }
81  else{
82  _error_(index1 <<" "<<index2 <<" "<<index3 <<" Not supported yet");
83  }
84 }

◆ ~GaussTetra()

GaussTetra::~GaussTetra ( )

Definition at line 86 of file GaussTetra.cpp.

86  {/*{{{*/
87  xDelete<IssmDouble>(weights);
88  xDelete<IssmDouble>(coords1);
89  xDelete<IssmDouble>(coords2);
90  xDelete<IssmDouble>(coords3);
91  xDelete<IssmDouble>(coords4);
92 }

Member Function Documentation

◆ begin()

int GaussTetra::begin ( void  )
virtual

Implements Gauss.

Definition at line 96 of file GaussTetra.cpp.

96  {/*{{{*/
97 
98  /*Check that this has been initialized*/
99  _assert_(numgauss>0);
100  _assert_(weights);
101  _assert_(coords1);
102  _assert_(coords2);
103  _assert_(coords3);
104  _assert_(coords4);
105 
106  /*return first gauss index*/
107  return 0;
108 }

◆ Echo()

void GaussTetra::Echo ( void  )
virtual

Implements Gauss.

Definition at line 110 of file GaussTetra.cpp.

110  {/*{{{*/
111 
112  _printf_("GaussTetra:\n");
113  _printf_(" numgauss: " << numgauss << "\n");
114 
115  if (weights){
116  _printf_(" weights = [");
117  for(int i=0;i<numgauss;i++) _printf_(" " << weights[i] << "\n");
118  _printf_("]\n");
119  }
120  else _printf_("weights = NULL\n");
121  if (coords1){
122  _printf_(" coords1 = [");
123  for(int i=0;i<numgauss;i++) _printf_(" " << coords1[i] << "\n");
124  _printf_("]\n");
125  }
126  else _printf_("coords1 = NULL\n");
127  if (coords2){
128  _printf_(" coords2 = [");
129  for(int i=0;i<numgauss;i++) _printf_(" " << coords2[i] << "\n");
130  _printf_("]\n");
131  }
132  else _printf_("coords2 = NULL\n");
133  if (coords3){
134  _printf_(" coords3 = [");
135  for(int i=0;i<numgauss;i++) _printf_(" " << coords3[i] << "\n");
136  _printf_("]\n");
137  }
138  else _printf_("coords3 = NULL\n");
139  if (coords4){
140  _printf_(" coords4 = [");
141  for(int i=0;i<numgauss;i++) _printf_(" " << coords4[i] << "\n");
142  _printf_("]\n");
143  }
144  else _printf_("coords4 = NULL\n");
145 
146  _printf_(" weight = " << weight << "\n");
147  _printf_(" coord1 = " << coord1 << "\n");
148  _printf_(" coord2 = " << coord2 << "\n");
149  _printf_(" coord3 = " << coord3 << "\n");
150  _printf_(" coord4 = " << coord4 << "\n");
151 
152 }

◆ end()

int GaussTetra::end ( void  )
virtual

Implements Gauss.

Definition at line 154 of file GaussTetra.cpp.

154  {/*{{{*/
155 
156  /*Check that this has been initialized*/
157  _assert_(numgauss>0);
158  _assert_(weights);
159  _assert_(coords1);
160  _assert_(coords2);
161  _assert_(coords3);
162  _assert_(coords4);
163 
164  /*return last gauss index +1*/
165  return numgauss;
166 }

◆ Enum()

int GaussTetra::Enum ( void  )
virtual

Implements Gauss.

Definition at line 168 of file GaussTetra.cpp.

168  {/*{{{*/
169  return GaussTetraEnum;
170 }

◆ GaussPoint()

void GaussTetra::GaussPoint ( int  ig)
virtual

Implements Gauss.

Definition at line 172 of file GaussTetra.cpp.

172  {/*{{{*/
173 
174  /*Check input in debugging mode*/
175  _assert_(ig>=0 && ig< numgauss);
176 
177  /*update static arrays*/
178  weight=weights[ig];
179  coord1=coords1[ig];
180  coord2=coords2[ig];
181  coord3=coords3[ig];
182  coord4=coords4[ig];
183 
184 }

◆ GaussNode()

void GaussTetra::GaussNode ( int  finitelement,
int  iv 
)
virtual

Implements Gauss.

Definition at line 186 of file GaussTetra.cpp.

186  {/*{{{*/
187 
188  /*in debugging mode: check that the default constructor has been called*/
189  _assert_(numgauss==-1);
190 
191  /*update static arrays*/
192  switch(finiteelement){
193  case P1Enum: case P1DGEnum:
194  switch(iv){
195  case 0: coord1=1.; coord2=0.; coord3=0.; coord4=0.; break;
196  case 1: coord1=0.; coord2=1.; coord3=0.; coord4=0.; break;
197  case 2: coord1=0.; coord2=0.; coord3=1.; coord4=0.; break;
198  case 3: coord1=0.; coord2=0.; coord3=0.; coord4=1.; break;
199  default: _error_("node index should be in [0 3]");
200  }
201  break;
202  case P1bubbleEnum: case P1bubblecondensedEnum:
203  switch(iv){
204  case 0: coord1=1.; coord2=0.; coord3=0.; coord4=0.; break;
205  case 1: coord1=0.; coord2=1.; coord3=0.; coord4=0.; break;
206  case 2: coord1=0.; coord2=0.; coord3=1.; coord4=0.; break;
207  case 3: coord1=0.; coord2=0.; coord3=0.; coord4=1.; break;
208  case 4: coord1=1./4.; coord2=1./4.; coord3=1./4.; coord4=1./4.; break;
209  default: _error_("node index should be in [0 4]");
210  }
211  break;
212  case P2Enum:
213  switch(iv){
214  case 0: coord1=1.; coord2=0.; coord3=0.; coord4=0.; break;
215  case 1: coord1=0.; coord2=1.; coord3=0.; coord4=0.; break;
216  case 2: coord1=0.; coord2=0.; coord3=1.; coord4=0.; break;
217  case 3: coord1=0.; coord2=0.; coord3=0.; coord4=1.; break;
218 
219  case 4: coord1=0.; coord2=.5; coord3=.5; coord4=0.; break;
220  case 5: coord1=.5; coord2=0.; coord3=.5; coord4=0.; break;
221  case 6: coord1=.5; coord2=.5; coord3=0.; coord4=0.; break;
222  case 7: coord1=.5; coord2=0.; coord3=0.; coord4=.5; break;
223  case 8: coord1=0.; coord2=.5; coord3=0.; coord4=.5; break;
224  case 9: coord1=0.; coord2=0.; coord3=.5; coord4=.5; break;
225  default: _error_("node index should be in [0 9]");
226  }
227  break;
228  default: _error_("Finite element "<<EnumToStringx(finiteelement)<<" not supported");
229  }
230 
231 }

◆ GaussVertex()

void GaussTetra::GaussVertex ( int  iv)
virtual

Implements Gauss.

Definition at line 233 of file GaussTetra.cpp.

233  {/*{{{*/
234 
235  /*in debugging mode: check that the default constructor has been called*/
236  _assert_(numgauss==-1);
237 
238  /*update static arrays*/
239  switch(iv){
240  case 0: coord1=1.; coord2=0.; coord3=0.; coord4=0.; break;
241  case 1: coord1=0.; coord2=1.; coord3=0.; coord4=0.; break;
242  case 2: coord1=0.; coord2=0.; coord3=1.; coord4=0.; break;
243  case 3: coord1=0.; coord2=0.; coord3=0.; coord4=1.; break;
244  default: _error_("vertex index should be in [0 3]");
245 
246  }
247 
248 }

◆ SynchronizeGaussBase()

void GaussTetra::SynchronizeGaussBase ( Gauss *  gauss)
virtual

Implements Gauss.

Definition at line 250 of file GaussTetra.cpp.

250  {/*{{{*/
251 
252  _error_("not supported");
253 }

Field Documentation

◆ numgauss

int GaussTetra::numgauss
private

Definition at line 15 of file GaussTetra.h.

◆ weights

IssmDouble* GaussTetra::weights
private

Definition at line 16 of file GaussTetra.h.

◆ coords1

IssmDouble* GaussTetra::coords1
private

Definition at line 17 of file GaussTetra.h.

◆ coords2

IssmDouble* GaussTetra::coords2
private

Definition at line 18 of file GaussTetra.h.

◆ coords3

IssmDouble* GaussTetra::coords3
private

Definition at line 19 of file GaussTetra.h.

◆ coords4

IssmDouble* GaussTetra::coords4
private

Definition at line 20 of file GaussTetra.h.

◆ coord1

IssmDouble GaussTetra::coord1

Definition at line 23 of file GaussTetra.h.

◆ coord2

IssmDouble GaussTetra::coord2

Definition at line 24 of file GaussTetra.h.

◆ coord3

IssmDouble GaussTetra::coord3

Definition at line 25 of file GaussTetra.h.

◆ coord4

IssmDouble GaussTetra::coord4

Definition at line 26 of file GaussTetra.h.

The documentation for this class was generated from the following files:
GaussTetra::weights
IssmDouble * weights
Definition: GaussTetra.h:16
_assert_
#define _assert_(ignore)
Definition: exceptions.h:37
IssmDouble
double IssmDouble
Definition: types.h:37
_printf_
#define _printf_(StreamArgs)
Definition: Print.h:22
GaussTetra::coords3
IssmDouble * coords3
Definition: GaussTetra.h:19
GaussTetra::coord4
IssmDouble coord4
Definition: GaussTetra.h:26
P1DGEnum
@ P1DGEnum
Definition: EnumDefinitions.h:1215
P1bubblecondensedEnum
@ P1bubblecondensedEnum
Definition: EnumDefinitions.h:1219
GaussTetraEnum
@ GaussTetraEnum
Definition: EnumDefinitions.h:1080
GaussTetra::coords4
IssmDouble * coords4
Definition: GaussTetra.h:20
P1Enum
@ P1Enum
Definition: EnumDefinitions.h:662
GaussTetra::coord1
IssmDouble coord1
Definition: GaussTetra.h:23
GaussTetra::coord3
IssmDouble coord3
Definition: GaussTetra.h:25
EnumToStringx
const char * EnumToStringx(int enum_in)
Definition: EnumToStringx.cpp:15
GaussTetra::coords2
IssmDouble * coords2
Definition: GaussTetra.h:18
UNDEF
#define UNDEF
Definition: constants.h:8
P1bubbleEnum
@ P1bubbleEnum
Definition: EnumDefinitions.h:1218
GaussLegendreTetra
void GaussLegendreTetra(int *pngaus, IssmDouble **pl1, IssmDouble **pl2, IssmDouble **pl3, IssmDouble **pl4, IssmDouble **pwgt, int iord)
Definition: GaussPoints.cpp:1212
GaussTetra::numgauss
int numgauss
Definition: GaussTetra.h:15
_error_
#define _error_(StreamArgs)
Definition: exceptions.h:49
GaussLegendreTria
void GaussLegendreTria(int *pngaus, IssmDouble **pl1, IssmDouble **pl2, IssmDouble **pl3, IssmDouble **pwgt, int iord)
Definition: GaussPoints.cpp:95
P2Enum
@ P2Enum
Definition: EnumDefinitions.h:1223
Gauss::weight
IssmDouble weight
Definition: Gauss.h:11
GaussTetra::coords1
IssmDouble * coords1
Definition: GaussTetra.h:17
GaussTetra::coord2
IssmDouble coord2
Definition: GaussTetra.h:24
Generated on Thu Jul 2 2020 08:09:27 for Ice Sheet System Model by   doxygen 1.8.19