Changeset 1612 for XIOS/dev/dev_olga/src/transformation/domain_algorithm_interpolate.cpp

Timestamp:

11/23/18 14:48:14 (6 years ago)

Author:

oabramkina

Message:

Dev: adding exception handling.

To activate it, compilation flag -DXIOS_EXCEPTION should be added.

File:

• XIOS/dev/dev_olga/src/transformation/domain_algorithm_interpolate.cpp (modified) (22 diffs)

XIOS/dev/dev_olga/src/transformation/domain_algorithm_interpolate.cpp

@@ -31 +31 @@
                                                                      std::map<int, int>& elementPositionInGridDst2AxisPosition,
                                                                      std::map<int, int>& elementPositionInGridDst2DomainPosition)
+TRY
 {
   std::vector<CDomain*> domainListDestP = gridDst->getDomains();
@@ -41 +42 @@
   return (new CDomainAlgorithmInterpolate(domainListDestP[domainDstIndex], domainListSrcP[domainSrcIndex], interpolateDomain));
 }
+CATCH
 
 bool CDomainAlgorithmInterpolate::registerTrans()
+TRY
 {
   CGridTransformationFactory<CDomain>::registerTransformation(TRANS_INTERPOLATE_DOMAIN, create);
 }
+CATCH
 
 CDomainAlgorithmInterpolate::CDomainAlgorithmInterpolate(CDomain* domainDestination, CDomain* domainSource, CInterpolateDomain* interpDomain)
 : CDomainAlgorithmTransformation(domainDestination, domainSource), interpDomain_(interpDomain), writeToFile_(false), readFromFile_(false)
+TRY
 {
   CContext* context = CContext::getCurrent();
@@ -93 +98 @@
     
 }
+CATCH
 
 /*!
@@ -98 +104 @@
 */
 void CDomainAlgorithmInterpolate::computeRemap()
+TRY
 {
   using namespace sphereRemap;
@@ -399 +406 @@
 
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::processPole(std::map<int,std::vector<std::pair<int,double> > >& interMapValuePole,
                                               int nbGlobalPointOnPole)
+TRY
 {
   CContext* context = CContext::getCurrent();
@@ -468 +477 @@
 
 }
+CATCH
 
 /*!
@@ -473 +483 @@
 */
 void CDomainAlgorithmInterpolate::computeIndexSourceMapping_(const std::vector<CArray<double,1>* >& dataAuxInputs)
+TRY
 {
   if (readFromFile_)  
@@ -481 +492 @@
   }
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::writeRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue)
+TRY
 {  
   writeInterpolationInfo(fileToReadWrite_, interpMapValue);
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::readRemapInfo()
+TRY
 {  
   std::map<int,std::vector<std::pair<int,double> > > interpMapValue;
@@ -494 +509 @@
   exchangeRemapInfo(interpMapValue);
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::convertRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue)
+TRY
 {
   CContext* context = CContext::getCurrent();
@@ -520 +537 @@
   }      
 }
+CATCH
 
 /*!
@@ -525 +543 @@
 */
 void CDomainAlgorithmInterpolate::exchangeRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue)
+TRY
 {
   CContext* context = CContext::getCurrent();
@@ -704 +723 @@
   delete [] recvBuff;
 }
+CATCH
  
 /*! Redefined some functions of CONetCDF4 to make use of them */
@@ -710 +730 @@
 int CDomainAlgorithmInterpolate::WriteNetCdf::addDimensionWrite(const StdString& name, 
                                                                 const StdSize size)
+TRY
 {
   return CONetCDF4::addDimension(name, size);  
 }
+CATCH
 
 int CDomainAlgorithmInterpolate::WriteNetCdf::addVariableWrite(const StdString& name, nc_type type,
                                                                const std::vector<StdString>& dim)
+TRY
 {
   return CONetCDF4::addVariable(name, type, dim);
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::WriteNetCdf::endDefinition()
+TRY
 {
   CONetCDF4::definition_end();
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::WriteNetCdf::writeDataIndex(const CArray<int,1>& data, const StdString& name,
@@ -729 +755 @@
                                                               const std::vector<StdSize>* start,
                                                               const std::vector<StdSize>* count)
+TRY
 {
   CONetCDF4::writeData<int,1>(data, name, collective, record, start, count);
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::WriteNetCdf::writeDataIndex(const CArray<double,1>& data, const StdString& name,
@@ -737 +765 @@
                                                               const std::vector<StdSize>* start,
                                                               const std::vector<StdSize>* count)
+TRY
 {
   CONetCDF4::writeData<double,1>(data, name, collective, record, start, count);
 }
+CATCH
 
 /*
@@ -748 +778 @@
 void CDomainAlgorithmInterpolate::writeInterpolationInfo(std::string& filename,
                                                          std::map<int,std::vector<std::pair<int,double> > >& interpMapValue)
+TRY
 {
   CContext* context = CContext::getCurrent();
@@ -827 +858 @@
   netCdfWriter.closeFile();
 }
+CATCH
 
 /*!
@@ -836 +868 @@
 void CDomainAlgorithmInterpolate::readInterpolationInfo(std::string& filename,
                                                         std::map<int,std::vector<std::pair<int,double> > >& interpMapValue)
+TRY
 {
   int ncid ;
@@ -894 +927 @@
       interpMapValue[dstIndex[ind]-indexOffset].push_back(make_pair(srcIndex[ind]-indexOffset,weight[ind]));
  }
+CATCH
 
 void CDomainAlgorithmInterpolate::apply(const std::vector<std::pair<int,double> >& localIndex,
@@ -900 +934 @@
                                             std::vector<bool>& flagInitial,
                                             bool ignoreMissingValue, bool firstPass  )
+TRY
 {
   int nbLocalIndex = localIndex.size();   
@@ -941 +976 @@
   }
 }
+CATCH
 
 void CDomainAlgorithmInterpolate::updateData(CArray<double,1>& dataOut)
+TRY
 {
   if (detectMissingValue)
@@ -961 +998 @@
   }
 }
-
-}
+CATCH
+
+}