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 POINT_LOCATOR_H
#define POINT_LOCATOR_H

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

#ifndef BLOCK_H
#include "SheafSystem/block.h"
#endif

#ifndef GEOMETRY_H
#include "SheafSystem/geometry.h"
#endif

namespace fiber_bundle
{
class chart_point;
class sec_ed;
}

namespace geometry
{

using namespace sheaf;
using namespace fiber_bundle;  

class SHEAF_DLL_SPEC point_locator
{

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

public:

  virtual bool invariant() const;

  virtual ~point_locator();



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

  sec_ed& coordinates() const;

  const block<sec_vd_value_type>& lb() const;

  const block<sec_vd_value_type>& ub() const;

  int dc() const;

  int db() const;

  bool domain_contains(sec_vd_value_type* xpt, size_type xpt_ub) const;

  virtual void update() = 0;

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

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

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

protected:

  point_locator();

  point_locator(sec_ed& xcoords);

  static sec_vd_value_type slightly_greater_than(sec_vd_value_type x);

  static sec_vd_value_type slightly_less_than(sec_vd_value_type x);


  block<sec_vd_value_type> _lb;

  block<sec_vd_value_type> _ub;

  int _dc;

  int _db;

  void update_domain();


private:

  point_locator(const point_locator& xother);

  virtual point_locator& operator=(const point_locator& xother);

  sec_ed& _coordinates;

};


SHEAF_DLL_SPEC 
std::ostream&
operator<<(std::ostream& xos, const point_locator& xpl);
 
} // namespace geometry

#endif // POINT_LOCATOR_H