Changeset 18521 for issm/trunk-jpl/src/c/classes/kriging/Observations.cpp

Timestamp:

09/15/14 16:17:28 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: removed ALL dynamic casts, now change the template shared/Numerics/recast.h to change back to dynamic_casts

File:

• issm/trunk-jpl/src/c/classes/kriging/Observations.cpp (modified) (8 diffs)

issm/trunk-jpl/src/c/classes/kriging/Observations.cpp

@@ -89 +89 @@
                 this->quadtree->ClosestObs(&index,x[i],y[i]);
                 if(index>=0){
-                        observation=dynamic_cast<Observation*>(this->GetObjectByOffset(index));
+                        observation=xDynamicCast<Observation*>(this->GetObjectByOffset(index));
                         if(pow(observation->x-x[i],2)+pow(observation->y-y[i],2) < minspacing) continue;
                 }
@@ -131 +131 @@
         this->quadtree->ClosestObs(&index,x_interp,y_interp);
         if(index>=0){
-                observation=dynamic_cast<Observation*>(this->GetObjectByOffset(index));
+                observation=xDynamicCast<Observation*>(this->GetObjectByOffset(index));
                 hmin = sqrt((observation->x-x_interp)*(observation->x-x_interp) + (observation->y-y_interp)*(observation->y-y_interp));
                 if(hmin<radius) radius=hmin;
@@ -139 +139 @@
         this->quadtree->RangeSearch(&indices,&nobs,x_interp,y_interp,radius);
         for (i=0;i<nobs;i++){
-                observation=dynamic_cast<Observation*>(this->GetObjectByOffset(indices[i]));
+                observation=xDynamicCast<Observation*>(this->GetObjectByOffset(indices[i]));
                 h2 = (observation->x-x_interp)*(observation->x-x_interp) + (observation->y-y_interp)*(observation->y-y_interp);
                 if(i==0){
@@ -155 +155 @@
         /*Assign output pointer*/
         if(nobs || hmin==radius){
-                observation=dynamic_cast<Observation*>(this->GetObjectByOffset(index));
+                observation=xDynamicCast<Observation*>(this->GetObjectByOffset(index));
                 *px   = observation->x;
                 *py   = observation->y;
@@ -210 +210 @@
         nobs = 0;
         for(i=0;i<tempnobs;i++){
-                observation=dynamic_cast<Observation*>(this->GetObjectByOffset(tempindices[i]));
+                observation=xDynamicCast<Observation*>(this->GetObjectByOffset(tempindices[i]));
                 h2 = (observation->x-x_interp)*(observation->x-x_interp) + (observation->y-y_interp)*(observation->y-y_interp);
 
@@ -252 +252 @@
                 /*Loop over all observations and fill in x, y and obs*/
                 for(i=0;i<nobs;i++){
-                        observation=dynamic_cast<Observation*>(this->GetObjectByOffset(indices[i]));
+                        observation=xDynamicCast<Observation*>(this->GetObjectByOffset(indices[i]));
                         observation->WriteXYObs(&x[i],&y[i],&obs[i]);
                 }
@@ -280 +280 @@
                 obs = xNew<IssmPDouble>(nobs);
                 for(int i=0;i<this->Size();i++){
-                        observation=dynamic_cast<Observation*>(this->GetObjectByOffset(i));
+                        observation=xDynamicCast<Observation*>(this->GetObjectByOffset(i));
                         observation->WriteXYObs(&x[i],&y[i],&obs[i]);
                 }
@@ -532 +532 @@
 
         for(i=0;i<this->Size();i++){
-                observation1=dynamic_cast<Observation*>(this->GetObjectByOffset(i));
+                observation1=xDynamicCast<Observation*>(this->GetObjectByOffset(i));
 
                 for(j=i+1;j<this->Size();j++){
-                        observation2=dynamic_cast<Observation*>(this->GetObjectByOffset(j));
+                        observation2=xDynamicCast<Observation*>(this->GetObjectByOffset(j));
 
                         distance=sqrt(pow(observation1->x - observation2->x,2) + pow(observation1->y - observation2->y,2));