basics/outputTecplot.hh

Go to the documentation of this file.

/* tecplot output function
 */

#ifndef OutputTecplot_hh
#define OutputTecplot_hh

#include <iostream>
#include <complex>
#include <vector>
#include "outputOperator.hh"
#include "vectorsMatrices.hh"
#include "toolbox/array.hh"

namespace concepts {

  // ********************************************************* OutputTecplot **

  template<typename F>
  class OutputTecplot : public OutputOperator {
  public:
    OutputTecplot(const F& val) : val_(val) {}
    void header(std::ostream& os, std::string name, uint cnt = 0);
  protected:
    std::ostream& info(std::ostream& os) const {
      return os << val_;
    }
  private:
    const F& val_;
  };

  template<typename F>
  void OutputTecplot<F>::header(std::ostream& os, std::string name, uint cnt) {
    os << name << cnt;
  }

  // *************************************** OutputTecplot<std::complex<F> > **

  template<typename F>
  class OutputTecplot<std::complex<F> > : public OutputOperator {
  public:
    OutputTecplot(const std::complex<F>& val) : val_(val) {}
    void header(std::ostream& os, std::string name, uint cnt = 0);
  protected:
    std::ostream& info(std::ostream& os) const {
      return os << OutputTecplot<F>(std::real(val_)) << " "
    << OutputTecplot<F>(std::imag(val_));
    }
  private:
    const std::complex<F>& val_;
  };

  template<>
  template<typename F>
  void OutputTecplot<std::complex<F> >::header
  (std::ostream& os, std::string name, uint cnt) {
    os << name << cnt << "r, " << name << cnt << "i";
  }

  // ********************************************** OutputTecplot<Array<F> > **

  template<typename F>
  class OutputTecplot<Array<F> > : public OutputOperator {
  public:
    OutputTecplot(const Array<F>& val) : val_(val) {}
    void header(std::ostream& os, std::string name, uint cnt = 0);
  protected:
    std::ostream& info(std::ostream& os) const {
      for(uint i = 0; i < val_.size(); ++i)
  os << OutputTecplot<F>(val_[i]) << " ";
      return os;
    }
  private:
    const Array<F>& val_;
  };

  template<>
  template<typename F>
  void OutputTecplot<Array<F> >::header
  (std::ostream& os, std::string name, uint cnt) {
    OutputTecplot<F> o(val_[0]);
    for(uint i = 0; i < val_.size(); ) {
      o.header(os, name, cnt+i);
      if (++i < val_.size()) os << ", ";
    }
  }

  // ****************************************** OutputTecplot<Point<F,dim> > **

  template<typename F, int dim>
  class OutputTecplot<Point<F,dim> > : public OutputOperator {
  public:
    OutputTecplot(const Point<F,dim>& val) : val_(val) {}
    void header(std::ostream& os, std::string name, uint cnt = 0);
  protected:
    std::ostream& info(std::ostream& os) const {
      for(uint i = 0; i < dim; ++i)
  os << OutputTecplot<F>(val_[i]) << " ";
      return os;
    }
  private:
    const Point<F,dim>& val_;
  };

  template<>
  template<typename F, int dim>
  void OutputTecplot<Point<F,dim> >::header
  (std::ostream& os, std::string name, uint cnt) {
    if (dim == 1)
      os << name << cnt;
    else {
      OutputTecplot<F> o(val_[0]);
      for(uint i = 0; i < dim; ) {
  o.header(os, name, cnt+i);
  if (++i < dim) os << ", ";
      }
    }
  }

} // namespace concepts

#endif // OutputTecplot_hh