Algebraic_Metrics_2.h

Go to the documentation of this file.

/* *******************************************************************
 * Rocstar Simulation Suite                                          *
 * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
 *                                                                   *
 * Illinois Rocstar LLC                                              *
 * Champaign, IL                                                     *
 * www.illinoisrocstar.com                                           *
 * sales@illinoisrocstar.com                                         *
 *                                                                   *
 * License: See LICENSE file in top level of distribution package or *
 * http://opensource.org/licenses/NCSA                               *
 *********************************************************************/
/* *******************************************************************
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
 * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
 * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
 * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
 * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
 *********************************************************************/
// $Id: Algebraic_Metrics_2.h,v 1.5 2008/12/06 08:45:23 mtcampbe Exp $

#ifndef __ALGEBRAIC_METRICS_2_H__
#define __ALGEBRAIC_METRICS_2_H__

#include "Metric.h"

MOP_BEGIN_NAMESPACE


class Alg_Metric_Base_2 : public Metric {
  friend std::ostream & operator<<(std::ostream &, 
                                   const Alg_Metric_Base_2 &);
public:
  
  Alg_Metric_Base_2() {}
  

  virtual ~Alg_Metric_Base_2() {}  


  virtual void initialize(Vector_3<double> n[], int type);

  virtual void initialize(Element_node_enumerator &ene);

  virtual double maxValue() const { return 1.0; }

  virtual double minValue() const { return 0.0; }

protected:


  double compute_size(double ref_area=1.) const;

  double compute_shape() const;

  double compute_skew() const;

protected:
  double   alpha[4];
  J_Matrix A[4];
  Matrix   L[4];
  int type_;
};

class Shape_Metric_2 : public Alg_Metric_Base_2 {
public:

  Shape_Metric_2() {}


  virtual void compute(double atts[]) const;
};

class Size_Metric_2 : public Alg_Metric_Base_2 {
public:


  explicit Size_Metric_2(double r_area) { ref_area=r_area; }


  virtual void compute(double atts[]) const;

private:
  double ref_area;
};

class Size_Shape_Metric_2 : public Alg_Metric_Base_2 {
public:


  explicit Size_Shape_Metric_2(double r_area) { ref_area=r_area; }


  virtual void compute(double atts[]) const;

  void set_ref_area(double r_area) { ref_area = r_area; }

private:
  double ref_area;
};

class Skew_Metric_2 : public Alg_Metric_Base_2 {
public:

  // Construct a Skew_Metric_2
  Skew_Metric_2() {}


  virtual void compute(double atts[]) const;
};

class Size_Skew_Metric_2 : public Alg_Metric_Base_2 {
public:


  explicit Size_Skew_Metric_2(double r_area) { ref_area=r_area; }
 

  virtual void compute(double atts[]) const;

private:
  double ref_area;
};

MOP_END_NAMESPACE

#endif