Changeset 24964
- Timestamp:
- 06/03/20 17:54:21 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r24963 r24964 5597 5597 #ifdef _HAVE_RESTRICT_ 5598 5598 IssmDouble* __restrict__ G=NULL; 5599 IssmDouble* __restrict__GU=NULL; 5600 IssmDouble* __restrict__GN=NULL; 5601 IssmDouble* __restrict__GE=NULL; 5602 IssmDouble* __restrict__G_elastic_precomputed=NULL; 5603 IssmDouble* __restrict__G_rigid_precomputed=NULL; 5604 IssmDouble* __restrict__U_elastic_precomputed=NULL; 5605 IssmDouble* __restrict__H_elastic_precomputed=NULL; 5599 IssmDouble* __restrict__ GU=NULL; 5600 IssmDouble* __restrict__ GN=NULL; 5601 IssmDouble* __restrict__ GE=NULL; 5602 IssmDouble* __restrict__ G_elastic_precomputed=NULL; 5603 IssmDouble* __restrict__ G_rigid_precomputed=NULL; 5604 IssmDouble* __restrict__ U_elastic_precomputed=NULL; 5605 IssmDouble* __restrict__ H_elastic_precomputed=NULL; 5606 IssmDouble* __restrict__ indices=NULL; 5606 5607 #else 5607 5608 IssmDouble* G=NULL; … … 5613 5614 IssmDouble* U_elastic_precomputed=NULL; 5614 5615 IssmDouble* H_elastic_precomputed=NULL; 5616 IssmDouble* indices=NULL; 5615 5617 #endif 5616 5618 5617 5619 /*elastic green function:*/ 5618 IssmDouble* indices=NULL;5619 5620 int index; 5620 5621 int M; … … 5636 5637 /*recover precomputed green function kernels:*/ 5637 5638 DoubleVecParam* parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelriseGRigidEnum)); _assert_(parameter); 5638 parameter->GetParameterValueByPointer( &G_rigid_precomputed,&M);5639 parameter->GetParameterValueByPointer((IssmDouble**)&G_rigid_precomputed,&M); 5639 5640 5640 5641 parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelriseGElasticEnum)); _assert_(parameter); 5641 parameter->GetParameterValueByPointer( &G_elastic_precomputed,&M);5642 parameter->GetParameterValueByPointer((IssmDouble**)&G_elastic_precomputed,&M); 5642 5643 5643 5644 parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelriseHElasticEnum)); _assert_(parameter); 5644 parameter->GetParameterValueByPointer( &H_elastic_precomputed,&M);5645 parameter->GetParameterValueByPointer((IssmDouble**)&H_elastic_precomputed,&M); 5645 5646 5646 5647 parameter = static_cast<DoubleVecParam*>(this->parameters->FindParamObject(SealevelriseUElasticEnum)); _assert_(parameter); 5647 parameter->GetParameterValueByPointer( &U_elastic_precomputed,&M);5648 parameter->GetParameterValueByPointer((IssmDouble**)&U_elastic_precomputed,&M); 5648 5649 5649 5650 /*allocate indices:*/
Note:
See TracChangeset
for help on using the changeset viewer.