---

CellJacobian.h

Go to the documentation of this file.

#ifndef CELLJACOBIAN_H
#define CELLJACOBIAN_H

#include "SundanceDefs.h"
#include "TSFError.h"
#include <string>


namespace Sundance
{

  using namespace TSF;
  using std::string;

  using std::ostream;

  /** \ingroup LowLevelGeometry
   * Jacobian of physical <---> reference coordinate transformations for
   * a cell.
   */

  class CellJacobian
    {
    public:
      /** empty ctor */
      inline CellJacobian();
      /** construct a jacobian for a 1D cell  */
      inline CellJacobian(const double& Jxu);
      /** construct a jacobian for a 2D cell  */
      inline CellJacobian(const double& Jxu, const double& Jxv,
                          const double& Jyu, const double& Jyv);
      /** construct a jacobian for a 3D cell  */
      inline CellJacobian(const double& Jxu, const double& Jxv, const double& Jxw,
                          const double& Jyu, const double& Jyv, const double& Jyw,
                          const double& Jzu, const double& Jzv, const double& Jzw);
      /** copy ctor */
      inline CellJacobian(const CellJacobian& other);
      /** assignment op */
      inline CellJacobian& operator=(const CellJacobian& other);

      /** get an element */
      inline const double& operator()(int phys, int ref) const ;
      /** return dimension of cell */
      inline int dim() const {return dim_;}

      /** run a pass-fail test */
      static bool unitTest() ;
    private:
      inline void boundsCheck(int i, int j) const ;
      double J_[9];
      int dim_;
    };

  ostream& operator<<(ostream& os, CellJacobian& J);

  inline CellJacobian::CellJacobian()
    : dim_(0)
    {
      ;
    }

  inline CellJacobian::CellJacobian(const double& Jxu)
    : dim_(1)
    {
      J_[0] = Jxu;
    }


  inline CellJacobian::CellJacobian(const double& Jxu, const double& Jxv,
                                    const double& Jyu, const double& Jyv)
    : dim_(2)
    {
      J_[0] = Jxu;
      J_[1] = Jxv;
      J_[2] = Jyu;
      J_[3] = Jyv;
    }



  inline CellJacobian::CellJacobian(const double& Jxu, const double& Jxv,
                                    const double& Jxw,
                                    const double& Jyu, const double& Jyv,
                                    const double& Jyw,
                                    const double& Jzu, const double& Jzv,
                                    const double& Jzw)
    : dim_(3)
    {
      J_[0] = Jxu;
      J_[1] = Jxv;
      J_[2] = Jxw;
      J_[3] = Jyu;
      J_[4] = Jyv;
      J_[5] = Jyw;
      J_[6] = Jzu;
      J_[7] = Jzv;
      J_[8] = Jzw;
    }

  inline CellJacobian::CellJacobian(const CellJacobian& other)
    : dim_(other.dim_)
    {
      int n = dim_*dim_;
      for (int i=0; i<n; i++) J_[i] = other.J_[i];
    }

  inline CellJacobian& CellJacobian::operator=(const CellJacobian& other)
    {
      if (this == &other) return *this;

      dim_ = other.dim_;

      int n = dim_*dim_;
      for (int i=0; i<n; i++) J_[i] = other.J_[i];

      return *this;
    }

  inline void CellJacobian::boundsCheck(int phys, int ref) const
    {
      if (phys >= dim_ || ref >= dim_) TSFError::raise("CellJacobian bounds error");
    }

  inline const double& CellJacobian::operator()(int phys, int ref) const
    {
#ifndef NOBOUNDSCHECK
      boundsCheck(phys, ref);
#endif
      return J_[ref + dim_*phys];
    }







}
#endif

---

Contact:

Kevin Long (krlong@ca.sandia.gov)


Documentation generated by