Rocsurf.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: Rocsurf.h,v 1.12 2008/12/06 08:43:23 mtcampbe Exp $

/*
 *  Las Modified:  Dec. 13, 2002
 */

#ifndef __ROCSURF_H_
#define __ROCSURF_H_

#include "surfbasic.h"
#include "roccom_devel.h"

SURF_BEGIN_NAMESPACE

class Window_manifold_2;

class Rocsurf : public COM_Object {
public:
  // protected:
  Rocsurf() : _wm(NULL), _cookie(SURF_COOKIE) {}
  virtual ~Rocsurf();

  static void load( const std::string &mname);
  static void unload( const std::string &mname);

  static void interpolate_to_centers( const COM::Attribute *x, 
                                      COM::Attribute *z);

  static void integrate( const COM::Attribute *x, double *z);

  static void compute_element_areas( COM::Attribute *areas, 
                                     const COM::Attribute *pnts=NULL);

  static void compute_element_normals( COM::Attribute *nrm,
                                       const int *to_normalize=NULL, 
                                       const COM::Attribute *pnts=NULL);

  static void compute_bounded_volumes( const COM::Attribute *old_location,
                                       const COM::Attribute *new_location,
                                       COM::Attribute *volumes,
                                       void *flag=NULL);

  static void compute_swept_volumes( const COM::Attribute *location,
                                     const COM::Attribute *disps,
                                     COM::Attribute *volumes,
                                     void *flag=NULL);

  static void compute_center( const COM::Attribute *mesh,
                              Vector_3<double> &cnt);

  static void compute_signed_volumes( const COM::Attribute *mesh,
                                      double *vol);

  virtual Window_manifold_2 *manifold() { return _wm; }

  void initialize( const COM::Attribute *pmesh);

  void compute_normals( const COM::Attribute *mesh,
                        COM::Attribute *nrm,
                        const int *scheme=NULL);

  void compute_mcn( COM::Attribute *mcn, COM::Attribute *lbmcn);

  void serialize_mesh( const COM::Attribute *inmesh, COM::Attribute *outmesh);

  void compute_edge_lengths( double *lave, double *lmin, double *lmax);
  
  void elements_to_nodes( const COM::Attribute *elem_vals,
                          COM::Attribute *nodal_vals,
                          const COM::Attribute *mesh=NULL,
                          const int *scheme=NULL,
                          const COM::Attribute *elem_weights=NULL,
                          COM::Attribute *nodal_weights=NULL);

protected:

  template <class T>
  static void normalize( T *a, int size) throw(int) {
    T tmp(0);
    for (int i=0; i<size; ++i) tmp += a[i]*a[i]; 
    
    if ( tmp == 0) return;
    tmp = std::sqrt(tmp);
    for (int i=0; i<size; ++i) a[i] /= tmp;
  }

  int validate_object() const {
    if ( _cookie != SURF_COOKIE) return -1;
    else return COM_Object::validate_object();
  }
  
protected:
  enum { SURF_COOKIE=7627873};
  Window_manifold_2   *_wm;
  static const int scheme_vals[];
  int _cookie;
};

SURF_END_NAMESPACE

#endif