JTable1D.hh

Go to the documentation of this file.

#ifndef __JTOOLS__JTABLE1D__
#define __JTOOLS__JTABLE1D__
 
#include "JMath/JMath.hh"
#include "JMath/JZero.hh"
#include "JIO/JSerialisable.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JTOOLS {}
namespace JPP { using namespace JTOOLS; }
 
namespace JTOOLS {
 
  using JMATH::JMath;
  using JIO::JReader;
  using JIO::JWriter;
 
 
  /**
   * 1D table with arithmetic capabilities.
   */
  template<unsigned int NX, class JData_t = double>
  struct JTable1D :
    public JMath< JTable1D<NX, JData_t> >
  {
    typedef JData_t      data_type;
 
    /**
     * Default constructor.
     */
    JTable1D()
    {
      for (int i = 0; i != NX; ++i) {
        data[i] = JMATH::zero;
      }
    }
 
 
    /**
     * Get number of entries.
     *
     * \return                 number of entries
     */
    static int getNX()
    { 
      return NX;
    }
 
 
    /**
     * Get entry data.
     *
     * \param  entry           entry number
     * \return                 entry data
     */
    const data_type& operator[](int entry) const
    { 
      return data[entry];
    };
 
 
    /**
     * Get entry data.
     *
     * \param  entry           entry number
     * \return                 entry data
     */
    data_type& operator[](int entry)
    {
      return data[entry];
    };
  
 
    /**
     * Negate table.
     *
     * \return                 this table
     */
    JTable1D& negate()
    {
      for (int i = 0; i != NX; ++i) {
        data[i] = -data[i];
      }
    
      return *this;
    }
  
  
    /**
     * Add table.
     *
     * \param  table           table
     * \return                 this table
     */
    JTable1D& add(const JTable1D& table)
    {
      for (int i = 0; i != NX; ++i) {
        data[i] += table.data[i];
      }
 
      return *this;
    }
  
  
    /**
     * Subtract table.
     *
     * \param  table           table
     * \return                 this table
     */
    JTable1D& sub(const JTable1D& table)
    {
      for (int i = 0; i != NX; ++i) {
        data[i] -= table.data[i];
      }
 
      return *this;
    } 
 
    /**
     * Scale table.
     *
     * \param  factor          multiplication factor
     * \return                 this table
     */
    JTable1D& mul(const double factor)
    {
      for (int i = 0; i != NX; ++i) {
        data[i] *= factor;
      }
    
      return *this;
    }
 
 
    /**
     * Scale table.
     *
     * \param  factor          division factor
     * \return                 this table
     */
    JTable1D& div(const double factor)
    {
      for (int i = 0; i != NX; ++i) {
        data[i] /= factor;
      }
    
      return *this;
    }
 
 
    /**
     * Read table from input.
     *
     * \param  in              reader
     * \param  table           table
     * \return                 reader
     */
    friend inline JReader& operator>>(JReader& in, JTable1D& table)
    {
      for (int i = 0; i != NX; ++i) {
        in >> table.data[i];
      }
 
      return in;
    }
 
 
    /**
     * Write table to output.
     *
     * \param  out             writer
     * \param  table           table
     * \return                 writer
     */
    friend inline JWriter& operator<<(JWriter& out, const JTable1D& table)
    {
      for (int i = 0; i != NX; ++i) {
        out << table.data[i];
      }
 
      return out;
    }
 
 
  protected:
    data_type data[NX];
  };
}
 
#endif