Propagation_3.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: Propagation_3.h,v 1.17 2008/12/06 08:45:27 mtcampbe Exp $

#ifndef _PROPAGATION_3_H_
#define _PROPAGATION_3_H_

#include <iostream>
#include <vector>
#include <map>
#include "propbasic.h"

PROP_BEGIN_NAMESPACE

const int MAX_LINES = 6;
const int BOUND_LEN = 15;

class Propagation_3 {
public: 
  Propagation_3()
    : _surf(NULL), _mode( SURF::ACROSS_PANE), _buf(NULL), _rank(0), 
      _verb(0), _cnstr_set(0), _cnstr_nodes(0), _cnstr_faces(0), 
      _cnstr_bndry_edges(0), _cnstr_bndry_nodes(), _cnstr_bound(0) {}

  virtual ~Propagation_3() {}

  Propagation_3( Manifold *wm, COM::Window *buf);

  virtual double time_stepping( const COM::Attribute *spd, double dt,
                                COM::Attribute *disps, int *smoothed=NULL)=0;

  virtual void set_constraints( const COM::Attribute *cnstr_types);

  void set_bounds( const COM::Attribute *bnd);

  // Bound a face based on nodal constraints. If all the nodes are beyond 
  // the bounded, then set the value to zero.
  void bound_facial_speed( COM::Attribute *fa);

  virtual void enforce_nodal_constraints( COM::Attribute *du);

  virtual void bound_nodal_motion( COM::Attribute *disps);

  void set_verbose( bool b) { _verb = b; }

protected:
  void convert_constraints( const COM::Attribute *ctypes_faces,
                            COM::Attribute *ctypes_nodes);
  void determine_constraint_boundary( const COM::Attribute *ctypes_faces,
                                      COM::Attribute *ctypes_bndry_edges,
                                      COM::Attribute *ctypes_bndry_nodes);

  static void get_constraint_directions( int type, int &ndirs, 
                                         Vector_3 dirs[2]);

  static void enforce_nodal_constraint( int type, Vector_3 &du);

  static void bound_nodal_motion( const Point_3 &pnt, const double *bnd, 
                                  Vector_3 &du, double eps=0);

  static bool check_spherical_bound( const Point_3 &pnt, const Point_3 &org, 
                                     const double rad_min, const double rad_max,
                                     double eps=0);
  static void bound_spherical_disp( const Point_3 &pnt, const Point_3 &org, 
                                    const double rad_min, const double rad_max, 
                                    Vector_3 &disp, double eps=0);

  static bool check_radial_bound( const double x, const double y,
                                  const double bnd_min, const double bnd_max,
                                  double eps=0);
  static void bound_radial_disp( const double x, const double y, 
                                 const double bnd_min, const double bnd_max,
                                 double &dx, double &dy, double eps=0);

  static bool check_axial_bound( const double x, const double bnd_min, 
                                 const double bnd_max, double eps=0);
  static void bound_axial_disp( const double x, const double bnd_min, 
                                const double bnd_max, double &dx, double eps=0);

  static bool reached_nodal_bound( const Point_3 &pnt, const double *bnd, 
                                   double eps=0);

  static bool in_bounding_box( const Point_3 &pnt, 
                               const Point_3 &lb, const Point_3 &ub) {
    return pnt[0] >= lb[0] && pnt[0] <= ub[0] &&
      pnt[1] >= lb[1] && pnt[1] <= ub[1] &&
      pnt[2] >= lb[2] && pnt[2] <= ub[2];
  }

protected:
  static double    square(double x) { return x*x; }
  Manifold         *_surf;
  SURF::Access_Mode _mode;      // Mode for accessing manifold
  
  COM::Window      *_buf;       // The buffer window created by Rocprop.

  int             _rank;        // Process rank.
  bool            _verb;        // Whether to constraints were set
  bool            _cnstr_set;   // Whether to constraints were set
  bool            _bnd_set;     // Whether to bound was set
  COM::Attribute *_cnstr_nodes; // Stores types of nodal constaints.
  COM::Attribute *_cnstr_faces; // Stores facial constaints
  COM::Attribute *_cnstr_bndry_edges;   // Stores edges along constraint boundaries
  COM::Attribute *_cnstr_bndry_nodes;   // Stores nodes along constraint boundaries
  COM::Attribute *_cnstr_bound; // Stores cylindrical constraint

  std::vector< COM::Pane*>      _panes; // panes of buffer window
};

PROP_END_NAMESPACE

#endif