Changeset 18521 for issm/trunk-jpl/src/c/classes/kriging/Observations.cpp
- Timestamp:
- 09/15/14 16:17:28 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/kriging/Observations.cpp
r18064 r18521 89 89 this->quadtree->ClosestObs(&index,x[i],y[i]); 90 90 if(index>=0){ 91 observation= dynamic_cast<Observation*>(this->GetObjectByOffset(index));91 observation=xDynamicCast<Observation*>(this->GetObjectByOffset(index)); 92 92 if(pow(observation->x-x[i],2)+pow(observation->y-y[i],2) < minspacing) continue; 93 93 } … … 131 131 this->quadtree->ClosestObs(&index,x_interp,y_interp); 132 132 if(index>=0){ 133 observation= dynamic_cast<Observation*>(this->GetObjectByOffset(index));133 observation=xDynamicCast<Observation*>(this->GetObjectByOffset(index)); 134 134 hmin = sqrt((observation->x-x_interp)*(observation->x-x_interp) + (observation->y-y_interp)*(observation->y-y_interp)); 135 135 if(hmin<radius) radius=hmin; … … 139 139 this->quadtree->RangeSearch(&indices,&nobs,x_interp,y_interp,radius); 140 140 for (i=0;i<nobs;i++){ 141 observation= dynamic_cast<Observation*>(this->GetObjectByOffset(indices[i]));141 observation=xDynamicCast<Observation*>(this->GetObjectByOffset(indices[i])); 142 142 h2 = (observation->x-x_interp)*(observation->x-x_interp) + (observation->y-y_interp)*(observation->y-y_interp); 143 143 if(i==0){ … … 155 155 /*Assign output pointer*/ 156 156 if(nobs || hmin==radius){ 157 observation= dynamic_cast<Observation*>(this->GetObjectByOffset(index));157 observation=xDynamicCast<Observation*>(this->GetObjectByOffset(index)); 158 158 *px = observation->x; 159 159 *py = observation->y; … … 210 210 nobs = 0; 211 211 for(i=0;i<tempnobs;i++){ 212 observation= dynamic_cast<Observation*>(this->GetObjectByOffset(tempindices[i]));212 observation=xDynamicCast<Observation*>(this->GetObjectByOffset(tempindices[i])); 213 213 h2 = (observation->x-x_interp)*(observation->x-x_interp) + (observation->y-y_interp)*(observation->y-y_interp); 214 214 … … 252 252 /*Loop over all observations and fill in x, y and obs*/ 253 253 for(i=0;i<nobs;i++){ 254 observation= dynamic_cast<Observation*>(this->GetObjectByOffset(indices[i]));254 observation=xDynamicCast<Observation*>(this->GetObjectByOffset(indices[i])); 255 255 observation->WriteXYObs(&x[i],&y[i],&obs[i]); 256 256 } … … 280 280 obs = xNew<IssmPDouble>(nobs); 281 281 for(int i=0;i<this->Size();i++){ 282 observation= dynamic_cast<Observation*>(this->GetObjectByOffset(i));282 observation=xDynamicCast<Observation*>(this->GetObjectByOffset(i)); 283 283 observation->WriteXYObs(&x[i],&y[i],&obs[i]); 284 284 } … … 532 532 533 533 for(i=0;i<this->Size();i++){ 534 observation1= dynamic_cast<Observation*>(this->GetObjectByOffset(i));534 observation1=xDynamicCast<Observation*>(this->GetObjectByOffset(i)); 535 535 536 536 for(j=i+1;j<this->Size();j++){ 537 observation2= dynamic_cast<Observation*>(this->GetObjectByOffset(j));537 observation2=xDynamicCast<Observation*>(this->GetObjectByOffset(j)); 538 538 539 539 distance=sqrt(pow(observation1->x - observation2->x,2) + pow(observation1->y - observation2->y,2));
Note:
See TracChangeset
for help on using the changeset viewer.