PN_patch.h File Reference

#include "mopbasic.h"

Include dependency graph for PN_patch.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
MOP_BEGIN_NAMESPACE Vector_3 < double >	PN_project (std::vector< int * > v_ids, Vector_3< double > bcoords, const Vector_3< double > *pnts, const Vector_3< double > *vnrms, const Vector_3< double > *evects, const Vector_3< double > *evals, const Vector_3< double > *bs, const int *tranks, std::vector< bool > is_ridge[4])
void	project_edge (int v_id, int t_id, int id3, Vector_3< double > v_crd, Vector_3< double > t_crd, Vector_3< double > &p_crd, Vector_3< double > &p2_crd, const Vector_3< double > *pnts, const Vector_3< double > *vnrms, const Vector_3< double > *evects, const Vector_3< double > *evals, const Vector_3< double > *bs, const int *tranks, bool is_ridge)
Vector_3< double >	one_sided_normal (int id1, int id2, int id3, const Vector_3< double > *evects, const Vector_3< double > *pnts, const Vector_3< double > *evals, const Vector_3< double > *vnorms, const Vector_3< double > *bs)

Function Documentation

Vector_3<double> one_sided_normal	(	int	id1,
		int	id2,
		int	id3,
		const Vector_3< double > *	evects,
		const Vector_3< double > *	pnts,
		const Vector_3< double > *	evals,
		const Vector_3< double > *	vnorms,
		const Vector_3< double > *	bs
	)

Definition at line 156 of file PN_patch.C.

References Vector_3< Type >::cross_product(), i, j, Vector_3< Type >::normalize(), sqrt(), and y.

Referenced by project_edge().

                                                             {
  int i = ((evects[3*id1-3]*bs[id1-1]/-evals[id1-1][0]) >
    (evects[3*id1-2]*bs[id1-1]/-evals[id1-1][1])) ? 0 : 1;
  int j = 1-i;
  Vector_3<double> os = std::sqrt(evals[id1-1][i]) * vnorms[id1-1];
  Vector_3<double> fnormal = Vector_3<double>::cross_product(pnts[id2-1]-pnts[id1-1],
                                                             pnts[id3-1]-pnts[id1-1]);
  Vector_3<double> y = Vector_3<double>::cross_product(vnorms[id1-1],evects[3*id1-1]);
  if(fnormal*y >= 0)
    os += std::sqrt(evals[id1-1][j])*y;
  else
    os -= std::sqrt(evals[id1-1][j])*y;
  os.normalize();
  return os;
}

Here is the call graph for this function:

Here is the caller graph for this function:

MOP_BEGIN_NAMESPACE Vector_3<double> PN_project	(	std::vector< int * >	v_ids,
		Vector_3< double >	bcoords,
		const Vector_3< double > *	pnts,
		const Vector_3< double > *	vnrms,
		const Vector_3< double > *	evects,
		const Vector_3< double > *	evals,
		const Vector_3< double > *	bs,
		const int *	tranks,
		std::vector< bool >	is_ridge[4]
	)

void project_edge	(	int	v_id,
		int	t_id,
		int	id3,
		Vector_3< double >	v_crd,
		Vector_3< double >	t_crd,
		Vector_3< double > &	p_crd,
		Vector_3< double > &	p2_crd,
		const Vector_3< double > *	pnts,
		const Vector_3< double > *	vnrms,
		const Vector_3< double > *	evects,
		const Vector_3< double > *	evals,
		const Vector_3< double > *	bs,
		const int *	tranks,
		bool	is_ridge
	)

Definition at line 103 of file PN_patch.C.

References COM_assertion_msg, Vector_3< Type >::normalize(), and one_sided_normal().

Referenced by PN_project().

{
  // convert from ids to indices
  v_id -=1;
  t_id -=1;
  Vector_3<double> vp_vect = (v_crd-p_crd);
  Vector_3<double> f(0.0,0.0,0.0);
  switch(tranks[v_id]) {
  case 2:{ // V is smooth    
    f = (vp_vect)*evects[3*v_id]*evects[3*v_id];
    break;
  }
  case 1:{ // V is a ridge
    if(is_ridge) // edge VT is on the ridge
      f = vp_vect - vp_vect*evects[3*v_id+2]*evects[3*v_id+2];
    else{ // edge VT is not part of the ridge
      Vector_3<double> n_os = 
        one_sided_normal(v_id+1,t_id+1,id3,
                         evects,
                         pnts,
                         evals,
                         vnrms,
                         bs);
      // one_sided_normal(v_id,n_os);
      f = vp_vect*n_os*n_os;
    }
    break;
  }
  case 0:{ // V is a corner, do nothing.
    return;
  }
  default :
    COM_assertion_msg(0, "Invalid tangent space size");
  }
  if(tranks[t_id]==0){// Other point is a corner, handle it here.
    Vector_3<double> tangent = p_crd - v_crd;
    tangent.normalize();
    p2_crd += f - 2.0*f*tangent*tangent;
  }
  p_crd += f;
}

Here is the call graph for this function:

Here is the caller graph for this function: