povW.h

Go to the documentation of this file.

/* This software and supporting documentation are distributed by
 *     Institut Federatif de Recherche 49
 *     CEA/NeuroSpin, Batiment 145,
 *     91191 Gif-sur-Yvette cedex
 *     France
 *
 * This software is governed by the CeCILL-B license under
 * French law and abiding by the rules of distribution of free software.
 * You can  use, modify and/or redistribute the software under the
 * terms of the CeCILL-B license as circulated by CEA, CNRS
 * and INRIA at the following URL "http://www.cecill.info".
 *
 * As a counterpart to the access to the source code and  rights to copy,
 * modify and redistribute granted by the license, users are provided only
 * with a limited warranty  and the software's author,  the holder of the
 * economic rights,  and the successive licensors  have only  limited
 * liability.
 *
 * In this respect, the user's attention is drawn to the risks associated
 * with loading,  using,  modifying and/or developing or reproducing the
 * software by the user in light of its specific status of free software,
 * that may mean  that it is complicated to manipulate,  and  that  also
 * therefore means  that it is reserved for developers  and  experienced
 * professionals having in-depth computer knowledge. Users are therefore
 * encouraged to load and test the software's suitability as regards their
 * requirements in conditions enabling the security of their systems and/or
 * data to be ensured and,  more generally, to use and operate it in the
 * same conditions as regards security.
 *
 * The fact that you are presently reading this means that you have had
 * knowledge of the CeCILL-B license and that you accept its terms.
 */

#ifndef AIMS_IO_POVW_H
#define AIMS_IO_POVW_H

#include <aims/config/aimsdata_config.h>
#include <aims/io/defaultItemW.h>
#include <aims/io/meshheader.h>
#include <aims/mesh/surface.h>
#include <aims/io/datatypecode.h>
#include <cartobase/exception/file.h>


namespace aims
{

  template<int D, class T=Void>
  class AIMSDATA_API PovWriter
  {
  public:
    PovWriter( const std::string & name ) 
      : _name( name ) { }
    ~PovWriter() {}

    inline void write( const AimsTimeSurface<D,T> & thing );

    inline std::string removeExtension(const std::string& name) const;

  private:
    std::string _name;
  };


  template <int D, class T>
  inline 
  AIMSDATA_API PovWriter<D,T> &
  operator << ( PovWriter<D,T> & writer, const AimsTimeSurface<D,T> & thing )
  {
    writer.write( thing );
    return( writer );
  }


  template <int D, class T> 
  std::string PovWriter<D,T>::removeExtension( const std::string& name ) const
  {
    std::string res = name;
    std::string ext="";
    if( res.length() > 4 )
      ext = res.substr( int(res.length() - 4), 4 );
    if( ext == ".pov" )
      res = res.substr( 0, res.length() - 4 );
    return res;
  }


  template <int D, class T>
  void PovWriter<D,T>::write( const AimsTimeSurface<D,T> & thing )
  {
    std::string fname = removeExtension( _name ) + ".pov";
    std::ofstream       os( fname.c_str() );
    if( !os )
      throw carto::file_error( fname );

    //  write data
    os << "mesh2 {" << std::endl;

    typename AimsTimeSurface<D,T>::const_iterator is = thing.begin();

    // Only fisrt time step is written
    if ( is != thing.end() ) {
      const std::vector<Point3df> &vert = (*is).second.vertex();
      const std::vector<Point3df> &norm = (*is).second.normal();
      const std::vector<AimsVector<uint,D> > &poly 
        = (*is).second.polygon();


      // vertices
      os << "  vertex_vectors {" << std::endl 
         << vert.size() << "," << std::endl << "    ";
      for( std::vector<Point3df>::const_iterator iv = vert.begin();
           iv != vert.end(); ++iv ) {
        os << "<" << (*iv)[0] << "," << (*iv)[1] << "," << (*iv)[2] << ">,";
      }
      os << "  }" << std::endl;
      
      // normals
      if ( ! norm.empty() ) {
        os << "  normal_vectors {" << std::endl 
           << norm.size() << "," << std::endl << "    ";
        for( std::vector<Point3df>::const_iterator in = norm.begin();
             in != norm.end(); ++in ) {
          os << "<" << (*in)[0] << "," << (*in)[1] << "," << (*in)[2] << ">,";
        }
        os << "  }" << std::endl;
      }
      
      // polygons
      os << "  face_indices {" << std::endl 
         << poly.size() << "," << std::endl << "    ";
      for( typename std::vector<AimsVector<uint,D> >::const_iterator ip = 
             poly.begin(); ip != poly.end(); ++ip ) {
        os << "<" << (*ip)[0] << "," << (*ip)[1] << "," << (*ip)[2] << ">,";
      }
      os << "  }" << std::endl;
      
      os << "  pigment { Gray80 }\n}" << std::endl;
    }
  }
}


#endif