/[escript]/trunk/escript/src/DataExpanded.cpp
ViewVC logotype

Diff of /trunk/escript/src/DataExpanded.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 965 by gross, Tue Feb 13 05:38:45 2007 UTC revision 966 by gross, Tue Feb 13 05:45:18 2007 UTC
# Line 617  DataExpanded::eigenvalues_and_eigenvecto Line 617  DataExpanded::eigenvalues_and_eigenvecto
617  void  void
618  DataExpanded::dump(const std::string fileName) const  DataExpanded::dump(const std::string fileName) const
619  {  {
    #ifdef PASO_MPI  
    throw DataException("Error - DataExpanded:: dump is not implemented for MPI yet.")  
    #endif  
    const NcDim* ncdims[2+DataArrayView::maxRank];  
    NcVar* var;  
    int rank = getPointDataView().getRank();  
    int type=  getFunctionSpace().getTypeCode();  
    int ndims =0;  
    long dims[2+DataArrayView::maxRank];  
    DataArrayView::ShapeType shape = getPointDataView().getShape();  
620    
621     // netCDF error handler     throw DataException("Error - DataExpanded:: dump is not implemented yet.");
    NcError err(NcError::verbose_nonfatal);  
    // Create the file.  
    NcFile dataFile(fileName.c_str(), NcFile::Replace);  
    // check if writing was successful  
    if (!dataFile.is_valid())  
         throw DataException("Error - DataExpanded:: opening of netCDF file for output failed.");  
    if (!dataFile.add_att("type","constant") )  
         throw DataException("Error - DataExpanded:: appending data type to netCDF file failed.");  
    if (!dataFile.add_att("rank",rank) )  
         throw DataException("Error - DataExpanded:: appending rank attribute to netCDF file failed.");  
    if (!dataFile.add_att("function_space_type",type))  
         throw DataException("Error - DataExpanded:: appending function space attribute to netCDF file failed.");  
    ndims=rank+2;  
    if ( rank >1 ) {  
        dims[0]=shape[0];  
        if (! (ncdims[0] = dataFile.add_dim("d0",shape[0])) )  
             throw DataException("Error - DataExpanded:: appending ncdimsion 0 to netCDF file failed.");  
    }  
    if ( rank >1 ) {  
        dims[1]=shape[1];  
        if (! (ncdims[1] = dataFile.add_dim("d1",shape[1])) )  
             throw DataException("Error - DataExpanded:: appending ncdimsion 1 to netCDF file failed.");  
    }  
    if ( rank >2 ) {  
        dims[2]=shape[2];  
        if (! (ncdims[2] = dataFile.add_dim("d2", shape[2])) )  
             throw DataException("Error - DataExpanded:: appending ncdimsion 2 to netCDF file failed.");  
    }  
    if ( rank >3 ) {  
        dims[3]=shape[3];  
        if (! (ncdims[3] = dataFile.add_dim("d3", shape[3])) )  
             throw DataException("Error - DataExpanded:: appending ncdimsion 3 to netCDF file failed.");  
    }  
    dims[rank]= getNumDataPointsPerSample();  
    if (! (ncdims[rank] = dataFile.add_dim("num_data_points_per_sample", dims[rank])) )  
             throw DataException("Error - DataExpanded:: appending num_data_points_per_sample to netCDF file failed.");  
    dims[rank+1]= getNumSamples();  
    if (! (ncdims[rank+1] = dataFile.add_dim("num_sample", dims[rank+1])) )  
             throw DataException("Error - DataExpanded:: appending num_sample to netCDF file failed.");  
    if (! ( var = dataFile.add_var("data", ncDouble, ndims, ncdims)) )  
         throw DataException("Error - DataExpanded:: appending variable to netCDF file failed.");  
    if (! (var->put(&m_data[0],dims)) )  
         throw DataException("Error - DataExpanded:: copy data to netCDF buffer failed.");  
622  }  }
623    
624    
625    
626  }  // end of namespace  }  // end of namespace
Legend:
Removed from v.965  
changed lines
  Added in v.966
  ViewVC Help
Powered by ViewVC 1.1.26