section_pusher.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 SECTION_PUSHER_H
#define SECTION_PUSHER_H

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

#ifndef ANY_H
#include "SheafSystem/any.h"
#endif

#ifndef FIELDS_H
#include "SheafSystem/fields.h"
#endif

#ifndef SECTION_PUSHER_PULLBACK_ACTION_H
#include "SheafSystem/section_pusher_pullback_action.h"
#endif

#ifndef SEC_ED_INVERTIBLE_H
#include "SheafSystem/sec_ed_invertible.h"
#endif

namespace fiber_bundle
{
class eval_iterator;
class sec_rep_space;
}

namespace fields
{

using namespace sheaf;
using namespace fiber_bundle;
 
class field_vd;
class field_vd;
class section_pusher_push_action;

class SHEAF_DLL_SPEC section_pusher : public any
{

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

public:

  //@todo: Delete the field_vd versions of the next 2 constructors when
  //       field_vd is removed from the build (remove the implementations also).

  section_pusher(const field_vd& xsrc, const field_vd& xdst, bool xauto_access);


//   ///
//   /// Creates the push-forward from the property space of
//   /// xsrc to a property space with the same fiber on the
//   /// base space of xdst_coords. Creates the range space
//   /// if it doesn't already exist.
//   ///
//   section_pusher(const field_vd& xsrc, const sec_ed& xdst_coords,
//                  bool xauto_access);

  section_pusher(const sec_rep_space& xdomain,
                 const sec_rep_space& xrange,
                 const sec_ed& xdomain_coord,
                 const sec_ed& xrange_coord,
                 bool xauto_access);

  const sec_rep_space& domain() const;

  const geometry::sec_ed_invertible& domain_coords() const;

  const sec_rep_space& range() const;

  const sec_ed& range_coords() const;

  sec_vd* push(const sec_vd& xinput, bool xauto_access);

  void push_pa(const sec_vd& xinput, sec_vd& result, bool xauto_access);

  bool same_base(const sec_ed& xcoord, const sec_rep_space& xspace) const;

  bool same_discretization(const sec_ed& xcoord, const sec_rep_space& xspace) const;

  bool same_fiber(const sec_rep_space& xspace1, const sec_rep_space& xspace2) const;

  section_pusher_push_action& overlap_action() const;

  section_pusher_push_action& default_overlap_action() const;

  section_pusher_push_action& extension_action() const;

  section_pusher_push_action& default_extension_action() const;

  void put_overlap_action(section_pusher_push_action& xoverlap_action);

  void put_extension_action(section_pusher_push_action& xextension_cction);

protected:

  section_pusher();

private:

  sec_rep_space* _domain;

  sec_rep_space* _range;

  geometry::sec_ed_invertible _domain_coords;

  sec_ed _range_coords;

  typedef section_pusher_pullback_action::pb_type pb_type;

  typedef section_pusher_pullback_action::pb_map_type pb_map_type;

  typedef section_pusher_pullback_action::pb_map_itr_type pb_map_itr_type;

  typedef section_pusher_pullback_action::pb_map_range_type pb_map_range_type;

  block<pb_map_type*> _pb_maps;

  virtual void make_range_monovalent_subposet(bool xis_ext);

  void pull_back_range_disc_mono_mono();

  void pull_back_range_disc_mono_multi();

  void pull_back_range_disc_multi_mono();

  typedef unordered::unordered_multimap<pod_index_type, scoped_index> chart_to_branch_map_type;

  chart_to_branch_map_type _chart_to_branch_map;

  scoped_index _first_branch;

  void make_maps(const sec_vd& xinput);

  void clear_maps();

  int _prop_df;

  int _src_db;

  block<sec_vd_dof_type> _src_dofs;

  block<sec_vd_dof_type> _src_value;

  block<sec_vd_dof_type> _dst_dofs;

  void define_overlap(const sec_vd& xinput, sec_vd& xresult);

  void push_eval_mbr(const eval_iterator& xeval_itr, const sec_vd& xinput,
                     sec_vd& xresult);

  void define_extension(sec_vd& xresult);

  section_pusher_push_action* _overlap_action;

  section_pusher_push_action* _default_overlap_action;

  section_pusher_push_action* _extension_action;

  section_pusher_push_action* _default_extension_action;

  void initialize();



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

public:

  section_pusher(const section_pusher& xother);

  virtual ~section_pusher();

  virtual bool invariant() const;


};


// ===========================================================
// NONMEMBER FUNCTIONS
// ===========================================================
 
} // namespace fields

#endif // ifndef SECTION_PUSHER_H