d_uniform_point_locator.h

Go to the documentation of this file.

//
// Copyright (c) 2014 Limit Point Systems, Inc.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//   http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
//


#ifndef D_UNIFORM_POINT_LOCATOR_H
#define D_UNIFORM_POINT_LOCATOR_H

#ifndef SHEAF_DLL_SPEC_H
#include "SheafSystem/sheaf_dll_spec.h"
#endif

#ifndef POINT_LOCATOR_H
#include "SheafSystem/point_locator.h"
#endif

//#undef DIAGNOSTIC_OUTPUT
//#define DIAGNOSTIC_OUTPUT

namespace fiber_bundle
{
class section_evaluator;
class structured_block;
}

namespace geometry
{

template <int DC, int DB>
class SHEAF_DLL_SPEC d_uniform_point_locator : public point_locator
{

  // ===========================================================
  // ===========================================================

public:

  d_uniform_point_locator(sec_ed& xcoords);

  virtual ~d_uniform_point_locator();

  virtual bool invariant() const;


protected:

  d_uniform_point_locator() {};

  d_uniform_point_locator(const d_uniform_point_locator<DC, DB>& xother) {};

    virtual d_uniform_point_locator& operator=(const point_locator& xother) {return *this;};

      d_uniform_point_locator& operator=(const d_uniform_point_locator<DC, DB>& xother) {return *this;};

  sec_vd_value_type _clipped_value[DC];

  structured_block* _block;

  chart_point* _block_pt;

  sec_vd_dof_type* _dofs;

  size_type _dof_ct;

  section_evaluator* _evaluator;

  sec_vd_value_type clip_component(sec_vd_value_type xcomp,
                                   sec_vd_value_type xupper_bound);



  // ===========================================================
  // ===========================================================

public:

  virtual void update();

  virtual void point_at_value(const sec_vd_value_type* xvalue,
                              size_type xvalue_ub,
                              chart_point& xresult);

  virtual void all_points_at_value(const sec_vd_value_type* xvalue,
                                   size_type xvalue_ub,
                                   block<chart_point_3d>& xresult);

  virtual void branch_points_at_value(const sec_vd_value_type* xvalue,
                                      size_type xvalue_ub,
                                      block<branch_point_pair>& xresult);


};


// ============================================================================
// SPECIALIZATIONS
// ============================================================================

template <>
SHEAF_DLL_SPEC
d_uniform_point_locator<1, 1>::
d_uniform_point_locator(sec_ed& xcoords);

template <>
SHEAF_DLL_SPEC
d_uniform_point_locator<2, 2>::
d_uniform_point_locator(sec_ed& xcoords);

template <>
SHEAF_DLL_SPEC
d_uniform_point_locator<3, 3>::
d_uniform_point_locator(sec_ed& xcoords);

} // namespace geometry

#endif // D_UNIFORM_POINT_LOCATOR_H