[2333] | 1 | /*!\file Numpar.c
|
---|
| 2 | * \brief: implementation of the Numpar object
|
---|
| 3 | */
|
---|
| 4 |
|
---|
| 5 | #ifdef HAVE_CONFIG_H
|
---|
| 6 | #include "config.h"
|
---|
| 7 | #else
|
---|
| 8 | #error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
|
---|
| 9 | #endif
|
---|
| 10 |
|
---|
| 11 | #include "stdio.h"
|
---|
| 12 | #include "./Numpar.h"
|
---|
| 13 | #include <string.h>
|
---|
| 14 | #include "../EnumDefinitions/EnumDefinitions.h"
|
---|
| 15 | #include "../shared/shared.h"
|
---|
| 16 | #include "../DataSet/DataSet.h"
|
---|
| 17 | #include "../include/typedefs.h"
|
---|
[3332] | 18 | #include "../include/macros.h"
|
---|
[2333] | 19 |
|
---|
[2907] | 20 | /*Object constructors and destructor*/
|
---|
| 21 | /*FUNCTION Numpar::constructor {{{1*/
|
---|
[2333] | 22 | Numpar::Numpar(){
|
---|
| 23 | return;
|
---|
| 24 | }
|
---|
[2907] | 25 | /*}}}*/
|
---|
| 26 | /*FUNCTION Numpar::creation {{{1*/
|
---|
[2333] | 27 | Numpar::Numpar(int numpar_id){
|
---|
| 28 | id=numpar_id;
|
---|
| 29 |
|
---|
| 30 | meanvel=UNDEF;
|
---|
| 31 | epsvel=UNDEF;
|
---|
| 32 | artdiff=UNDEF;
|
---|
| 33 | viscosity_overshoot=UNDEF;
|
---|
| 34 | stokesreconditioning=UNDEF;
|
---|
[2347] | 35 | control_type=NULL;
|
---|
[2354] | 36 | cm_noisedmp=UNDEF;
|
---|
| 37 | cm_mindmp_value=UNDEF;
|
---|
| 38 | cm_mindmp_slope=UNDEF;
|
---|
| 39 | cm_maxdmp_value=UNDEF;
|
---|
| 40 | cm_maxdmp_slope=UNDEF;
|
---|
[2333] | 41 |
|
---|
| 42 | return;
|
---|
| 43 | }
|
---|
[2907] | 44 | /*}}}*/
|
---|
| 45 | /*FUNCTION Numpar::destructor {{{1*/
|
---|
[2333] | 46 | Numpar::~Numpar(){
|
---|
[3595] | 47 |
|
---|
| 48 | /*Free the only pointer*/
|
---|
[3607] | 49 | //xfree((void**)&control_type);
|
---|
[3595] | 50 |
|
---|
[2333] | 51 | return;
|
---|
| 52 | }
|
---|
[2907] | 53 | /*}}}*/
|
---|
[2333] | 54 |
|
---|
[2907] | 55 | /*Object marshall*/
|
---|
| 56 | /*FUNCTION Numpar::Marshall {{{1*/
|
---|
[2333] | 57 | void Numpar::Marshall(char** pmarshalled_dataset){
|
---|
| 58 |
|
---|
| 59 | char* marshalled_dataset=NULL;
|
---|
| 60 | int enum_type=0;
|
---|
| 61 |
|
---|
| 62 | /*recover marshalled_dataset: */
|
---|
| 63 | marshalled_dataset=*pmarshalled_dataset;
|
---|
| 64 |
|
---|
| 65 | /*get enum type of Numpar: */
|
---|
[3567] | 66 | enum_type=NumparEnum;
|
---|
[2333] | 67 |
|
---|
| 68 | /*marshall enum: */
|
---|
| 69 | memcpy(marshalled_dataset,&enum_type,sizeof(enum_type));marshalled_dataset+=sizeof(enum_type);
|
---|
| 70 |
|
---|
| 71 | /*marshall Numpar data: */
|
---|
| 72 | memcpy(marshalled_dataset,&id,sizeof(id));marshalled_dataset+=sizeof(id);
|
---|
| 73 | memcpy(marshalled_dataset,&meanvel,sizeof(meanvel));marshalled_dataset+=sizeof(meanvel);
|
---|
| 74 | memcpy(marshalled_dataset,&epsvel,sizeof(epsvel));marshalled_dataset+=sizeof(epsvel);
|
---|
| 75 | memcpy(marshalled_dataset,&artdiff,sizeof(artdiff));marshalled_dataset+=sizeof(artdiff);
|
---|
| 76 | memcpy(marshalled_dataset,&viscosity_overshoot,sizeof(viscosity_overshoot));marshalled_dataset+=sizeof(viscosity_overshoot);
|
---|
| 77 | memcpy(marshalled_dataset,&stokesreconditioning,sizeof(stokesreconditioning));marshalled_dataset+=sizeof(stokesreconditioning);
|
---|
[2347] | 78 | memcpy(marshalled_dataset,&control_type,sizeof(control_type));marshalled_dataset+=sizeof(control_type);
|
---|
[2354] | 79 | memcpy(marshalled_dataset,&cm_noisedmp,sizeof(cm_noisedmp));marshalled_dataset+=sizeof(cm_noisedmp);
|
---|
| 80 | memcpy(marshalled_dataset,&cm_mindmp_value,sizeof(cm_mindmp_value));marshalled_dataset+=sizeof(cm_mindmp_value);
|
---|
| 81 | memcpy(marshalled_dataset,&cm_mindmp_slope,sizeof(cm_mindmp_slope));marshalled_dataset+=sizeof(cm_mindmp_slope);
|
---|
| 82 | memcpy(marshalled_dataset,&cm_maxdmp_value,sizeof(cm_maxdmp_value));marshalled_dataset+=sizeof(cm_maxdmp_value);
|
---|
| 83 | memcpy(marshalled_dataset,&cm_maxdmp_slope,sizeof(cm_maxdmp_slope));marshalled_dataset+=sizeof(cm_maxdmp_slope);
|
---|
[2333] | 84 |
|
---|
| 85 | *pmarshalled_dataset=marshalled_dataset;
|
---|
| 86 | return;
|
---|
| 87 | }
|
---|
[2907] | 88 | /*}}}*/
|
---|
| 89 | /*FUNCTION Numpar::MarshallSize{{{1*/
|
---|
[2333] | 90 | int Numpar::MarshallSize(){
|
---|
| 91 | return sizeof(id)
|
---|
| 92 | +sizeof(meanvel)
|
---|
| 93 | +sizeof(epsvel)
|
---|
| 94 | +sizeof(artdiff)
|
---|
| 95 | +sizeof(viscosity_overshoot)
|
---|
| 96 | +sizeof(stokesreconditioning)
|
---|
[2347] | 97 | +sizeof(control_type)
|
---|
[2354] | 98 | +sizeof(cm_noisedmp)
|
---|
| 99 | +sizeof(cm_mindmp_value)
|
---|
| 100 | +sizeof(cm_mindmp_slope)
|
---|
| 101 | +sizeof(cm_maxdmp_value)
|
---|
| 102 | +sizeof(cm_maxdmp_slope)
|
---|
[2333] | 103 | +sizeof(int); //sizeof(int) for enum type
|
---|
| 104 | }
|
---|
[2907] | 105 | /*}}}*/
|
---|
| 106 | /*FUNCTION Numpar::Demarshall {{{1*/
|
---|
[2333] | 107 | void Numpar::Demarshall(char** pmarshalled_dataset){
|
---|
| 108 |
|
---|
| 109 | char* marshalled_dataset=NULL;
|
---|
| 110 | int i;
|
---|
| 111 |
|
---|
| 112 | /*recover marshalled_dataset: */
|
---|
| 113 | marshalled_dataset=*pmarshalled_dataset;
|
---|
| 114 |
|
---|
| 115 | /*this time, no need to get enum type, the pointer directly points to the beginning of the
|
---|
| 116 | *object data (thanks to DataSet::Demarshall):*/
|
---|
| 117 |
|
---|
| 118 | memcpy(&id,marshalled_dataset,sizeof(id));marshalled_dataset+=sizeof(id);
|
---|
| 119 | memcpy(&meanvel,marshalled_dataset,sizeof(meanvel));marshalled_dataset+=sizeof(meanvel);
|
---|
| 120 | memcpy(&epsvel,marshalled_dataset,sizeof(epsvel));marshalled_dataset+=sizeof(epsvel);
|
---|
| 121 | memcpy(&artdiff,marshalled_dataset,sizeof(artdiff));marshalled_dataset+=sizeof(artdiff);
|
---|
| 122 | memcpy(&viscosity_overshoot,marshalled_dataset,sizeof(viscosity_overshoot));marshalled_dataset+=sizeof(viscosity_overshoot);
|
---|
| 123 | memcpy(&stokesreconditioning,marshalled_dataset,sizeof(stokesreconditioning));marshalled_dataset+=sizeof(stokesreconditioning);
|
---|
[2347] | 124 | memcpy(&control_type,marshalled_dataset,sizeof(control_type));marshalled_dataset+=sizeof(control_type);
|
---|
[2354] | 125 | memcpy(&cm_noisedmp,marshalled_dataset,sizeof(cm_noisedmp));marshalled_dataset+=sizeof(cm_noisedmp);
|
---|
| 126 | memcpy(&cm_mindmp_value,marshalled_dataset,sizeof(cm_mindmp_value));marshalled_dataset+=sizeof(cm_mindmp_value);
|
---|
| 127 | memcpy(&cm_mindmp_slope,marshalled_dataset,sizeof(cm_mindmp_slope));marshalled_dataset+=sizeof(cm_mindmp_slope);
|
---|
| 128 | memcpy(&cm_maxdmp_value,marshalled_dataset,sizeof(cm_maxdmp_value));marshalled_dataset+=sizeof(cm_maxdmp_value);
|
---|
| 129 | memcpy(&cm_maxdmp_slope,marshalled_dataset,sizeof(cm_maxdmp_slope));marshalled_dataset+=sizeof(cm_maxdmp_slope);
|
---|
[2333] | 130 |
|
---|
| 131 | /*return: */
|
---|
| 132 | *pmarshalled_dataset=marshalled_dataset;
|
---|
| 133 | return;
|
---|
| 134 | }
|
---|
[2907] | 135 | /*}}}*/
|
---|
[2333] | 136 |
|
---|
[2907] | 137 | /*Object functions*/
|
---|
| 138 | /*FUNCTION Numpar::Configure {{{1*/
|
---|
[2333] | 139 | void Numpar::Configure(void* pparametersin){
|
---|
| 140 |
|
---|
| 141 | DataSet* parameters=NULL;
|
---|
| 142 |
|
---|
| 143 | /*Recover virtual pointer:*/
|
---|
| 144 | parameters=(DataSet*)pparametersin;
|
---|
| 145 |
|
---|
| 146 | /*Go through parameters dataset, and find the Param object corresponding to our fields,
|
---|
| 147 | * and update the fields: */
|
---|
[2354] | 148 | parameters->FindParam(&meanvel,"meanvel");
|
---|
| 149 | parameters->FindParam(&epsvel,"epsvel");
|
---|
[3332] | 150 | if(!parameters->FindParam(&artdiff,"artdiff"))ISSMERROR(" error message: could not update artdiff field");
|
---|
| 151 | if(!parameters->FindParam(&viscosity_overshoot,"viscosity_overshoot"))ISSMERROR(" error message: could not update viscosity_overshoot field");
|
---|
| 152 | if(!parameters->FindParam(&stokesreconditioning,"stokesreconditioning"))ISSMERROR(" error message: could not update stokesreconditioning field");
|
---|
[3053] | 153 | xfree((void**)&control_type);
|
---|
| 154 | parameters->FindParam(&control_type,"control_type");
|
---|
[2354] | 155 | parameters->FindParam(&cm_noisedmp,"cm_noisedmp");
|
---|
| 156 | parameters->FindParam(&cm_mindmp_value,"cm_mindmp_value");
|
---|
| 157 | parameters->FindParam(&cm_mindmp_slope,"cm_mindmp_slope");
|
---|
| 158 | parameters->FindParam(&cm_maxdmp_value,"cm_maxdmp_value");
|
---|
| 159 | parameters->FindParam(&cm_maxdmp_slope,"cm_maxdmp_slope");
|
---|
[2333] | 160 |
|
---|
| 161 | return;
|
---|
| 162 | }
|
---|
[2907] | 163 | /*}}}*/
|
---|
| 164 | /*FUNCTION Numpar::copy {{{1*/
|
---|
| 165 | Object* Numpar::copy() {
|
---|
| 166 |
|
---|
| 167 | return new Numpar(*this);
|
---|
[2333] | 168 |
|
---|
[2907] | 169 | }
|
---|
| 170 | /*}}}*/
|
---|
| 171 | /*FUNCTION Numpar::DeepEcho {{{1*/
|
---|
| 172 | void Numpar::DeepEcho(void){
|
---|
[2333] | 173 |
|
---|
[2907] | 174 | printf("Numpar:\n");
|
---|
| 175 | printf(" id: %i\n",id);
|
---|
| 176 | printf(" meanvel: %g\n",meanvel);
|
---|
| 177 | printf(" epsvel: %g\n",epsvel);
|
---|
| 178 | printf(" artdiff: %i\n",artdiff);
|
---|
| 179 | printf(" viscosity_overshoot: %g\n",viscosity_overshoot);
|
---|
| 180 | printf(" stokesreconditioning: %g\n",stokesreconditioning);
|
---|
| 181 | printf(" control_type: %s\n",control_type);
|
---|
| 182 | printf(" cm_noisedmp: %g\n",cm_noisedmp);
|
---|
| 183 | printf(" cm_mindmp_value: %g\n",cm_mindmp_value);
|
---|
| 184 | printf(" cm_mindmp_slope: %g\n",cm_mindmp_slope);
|
---|
| 185 | printf(" cm_maxdmp_value: %g\n",cm_maxdmp_value);
|
---|
| 186 | printf(" cm_maxdmp_slope: %g\n",cm_maxdmp_slope);
|
---|
| 187 | }
|
---|
| 188 | /*}}}*/
|
---|
| 189 | /*FUNCTION Numpar::Echo{{{1*/
|
---|
| 190 | void Numpar::Echo(void){
|
---|
| 191 |
|
---|
| 192 | printf("Numpar:\n");
|
---|
| 193 | printf(" id: %i\n",id);
|
---|
| 194 | printf(" meanvel: %g\n",meanvel);
|
---|
| 195 | printf(" epsvel: %g\n",epsvel);
|
---|
| 196 | printf(" artdiff: %i\n",artdiff);
|
---|
| 197 | printf(" viscosity_overshoot: %g\n",viscosity_overshoot);
|
---|
| 198 | printf(" stokesreconditioning: %g\n",stokesreconditioning);
|
---|
| 199 | printf(" control_type: %s\n",control_type);
|
---|
| 200 | printf(" cm_noisedmp: %g\n",cm_noisedmp);
|
---|
| 201 | printf(" cm_mindmp_value: %g\n",cm_mindmp_value);
|
---|
| 202 | printf(" cm_mindmp_slope: %g\n",cm_mindmp_slope);
|
---|
| 203 | printf(" cm_maxdmp_value: %g\n",cm_maxdmp_value);
|
---|
| 204 | printf(" cm_maxdmp_slope: %g\n",cm_maxdmp_slope);
|
---|
| 205 | }
|
---|
| 206 | /*}}}*/
|
---|
| 207 | /*FUNCTION Numpar::Enum {{{1*/
|
---|
| 208 | int Numpar::Enum(void){
|
---|
| 209 |
|
---|
[3567] | 210 | return NumparEnum;
|
---|
[2907] | 211 |
|
---|
| 212 | }
|
---|
| 213 | /*}}}*/
|
---|
| 214 | /*FUNCTION Numpar::GetId {{{1*/
|
---|
| 215 | int Numpar::GetId(void){ return id; }
|
---|
| 216 | /*}}}*/
|
---|
| 217 | /*FUNCTION Numpar::GetName {{{1*/
|
---|
| 218 | char* Numpar::GetName(void){
|
---|
| 219 | return "beam";
|
---|
| 220 | }
|
---|
| 221 | /*}}}*/
|
---|
| 222 | /*FUNCTION Numpar::MyRank {{{1*/
|
---|
| 223 | int Numpar::MyRank(void){
|
---|
| 224 | extern int my_rank;
|
---|
| 225 | return my_rank;
|
---|
| 226 | }
|
---|
| 227 | /*}}}*/
|
---|
| 228 | /*FUNCTION Numpar::UpdateFromInputs {{{1*/
|
---|
[2333] | 229 | void Numpar::UpdateFromInputs(void* vinputs){
|
---|
| 230 |
|
---|
| 231 | ParameterInputs* inputs=NULL;
|
---|
| 232 |
|
---|
| 233 | /*recover pointers: */
|
---|
| 234 | inputs=(ParameterInputs*)vinputs;
|
---|
| 235 |
|
---|
| 236 | /*Update internal data if inputs holds new values: */
|
---|
| 237 | inputs->Recover("meanvel",&meanvel);
|
---|
| 238 | inputs->Recover("epsvel",&epsvel);
|
---|
| 239 | inputs->Recover("artdiff",&artdiff);
|
---|
| 240 | inputs->Recover("viscosity_overshoot",&viscosity_overshoot);
|
---|
| 241 | inputs->Recover("stokesreconditioning",&stokesreconditioning);
|
---|
[2347] | 242 | inputs->Recover("control_type",&control_type);
|
---|
[2354] | 243 | inputs->Recover("cm_noisedmp",&cm_noisedmp);
|
---|
| 244 | inputs->Recover("cm_mindmp_value",&cm_mindmp_value);
|
---|
| 245 | inputs->Recover("cm_mindmp_slope",&cm_mindmp_slope);
|
---|
| 246 | inputs->Recover("cm_maxdmp_value",&cm_maxdmp_value);
|
---|
| 247 | inputs->Recover("cm_maxdmp_slope",&cm_maxdmp_slope);
|
---|
[2333] | 248 |
|
---|
| 249 | }
|
---|
[2907] | 250 | /*}}}*/
|
---|