JMatrix1S.hh

Go to the documentation of this file.

#ifndef __JMATRIX1S__
#define __JMATRIX1S__
 
#include <cmath>
#include <algorithm>
 
#include "JMath/JMatrix1D.hh"
#include "JLang/JException.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JMATH {}
namespace JPP { using namespace JMATH; }
 
namespace JMATH {
 
  using JLANG::JDivisionByZero;
 
 
  /**
   * 1 x 1 symmetric matrix
   */
  class JMatrix1S :
    public JMatrix1D
  {
  public:
    /**
     * Default constructor.
     */
    JMatrix1S() :
      JMatrix1D()
    {}
 
 
    /**
     * Contructor.
     *
     * \param   A      matrix
     */
    JMatrix1S(const JMatrix1D& A) :
      JMatrix1D(A)
    {}
 
 
    /**
     * Contructor.
     *
     * \param   __a00   (0,0)
     */
    JMatrix1S(const double __a00) :
      JMatrix1D(__a00)
    {}
 
 
    /**
     * Invert matrix
     */
    void invert()
    {
      if (a00 == 0) {
        throw JDivisionByZero("LDU decomposition zero pivot");
      }
 
      a00 = 1.0 / a00;
    }
  };
}
 
#endif