app-radu/tree.hh

Go to the documentation of this file.

#ifndef appRaduTree_hh
#define appRaduTree_hh

#include "basics/typedefs.hh"
#include "basics/vectorsMatricesForward.hh"
#include "basics/exceptions.hh"
#include "space/space.hh"
#include "app-radu/element.hh"
#include "cluster/tree.hh"

// ***************************************************************** cebysev **

void cebysev(const concepts::Real3d** vtx, uint nvtx, concepts::Real3d& c,
       concepts::Real& r, const concepts::Real3d* (&e)[4], uint& ne);

// **************************************************************** BBall001 **

template<class F>
class BBall001 : public cluster::BBall<F> {
public:
  void operator()(const concepts::Element<F>& elm, concepts::Real3d& c,
      concepts::Real& r) const;
};

template<class F>
void BBall001<F>::operator()(const concepts::Element<F>& elm,
           concepts::Real3d& c, concepts::Real& r) const {

  const Constant3d000<F>* elmA =
    dynamic_cast<const Constant3d000<F>*>(&elm);
  if (elmA != 0) {
    const concepts::Real3d *v[4], *e[4];
    v[0] = &elmA->vertex(0);  v[1] = &elmA->vertex(1);
    v[2] = &elmA->vertex(2);  v[3] = &elmA->vertex(3);
    uint ne;
    cebysev(v, 4, c, r, e, ne);
    return;
  }

  const Constant3d001<F>* elmB =
    dynamic_cast<const Constant3d001<F>*>(&elm);
  if (elmB != 0) {
    const concepts::Real3d *v[8], *e[4];
    v[0] = &elmB->vertex(0);  v[1] = &elmB->vertex(1);
    v[2] = &elmB->vertex(2);  v[3] = &elmB->vertex(3);
    v[4] = &elmB->vertex(4);  v[5] = &elmB->vertex(5);
    v[6] = &elmB->vertex(6);  v[7] = &elmB->vertex(7);
    uint ne;
    cebysev(v, 8, c, r, e, ne);
    return;
  }

  const Constant3d002<F>* elmC =
    dynamic_cast<const Constant3d002<F>*>(&elm);
  if (elmC != 0) {
    const concepts::Real3d *v[8], *e[4];
    v[0] = &elmC->vertex(0);  v[1] = &elmC->vertex(1);
    v[2] = &elmC->vertex(2);  v[3] = &elmC->vertex(3);
    v[4] = &elmC->vertex(4);  v[5] = &elmC->vertex(5);
    v[6] = &elmC->vertex(6);  v[7] = &elmC->vertex(7);
    uint ne;
    cebysev(v, 8, c, r, e, ne);
    return;
  }

  throw conceptsException(concepts::MissingFeature("cell not supported"));
}


#endif // appRaduTree_hh