eval_iterator.h

//
// 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.
//

// Interface for class eval_iterator


#ifndef EVAL_ITERATOR_H
#define EVAL_ITERATOR_H

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

#ifndef DEPTH_FIRST_ITERATOR_H
#include "SheafSystem/depth_first_iterator.h"
#endif

#ifndef BASE_SPACE_POSET_H
#include "SheafSystem/base_space_poset.h"
#endif

#ifndef STD_LIST_H
#include "SheafSystem/std_list.h"
#endif

#ifndef SEC_VD_H
#include "SheafSystem/sec_vd.h" // just to get typedefs
#endif

namespace sheaf
{
class index_space_handle;
}

namespace fiber_bundle
{

using namespace sheaf;


class eval_family;
class section_space_schema_member;
class section_evaluator;
class sec_rep_descriptor;

class SHEAF_DLL_SPEC eval_iterator : public depth_first_iterator
{
public:

  // CANONICAL MEMBERS

  eval_iterator();

  eval_iterator(const eval_iterator& xother);

  ~eval_iterator();

  virtual bool is_ancestor_of(const any* xother) const;

  virtual eval_iterator* clone() const;

  bool invariant() const;

  // OTHER CONSTRUCTORS

  eval_iterator(const section_space_schema_member& xschema_anchor, bool xvisit_once = true);

  // ITERATOR FACET

  virtual bool is_initialized() const;

  virtual void force_is_done();

  inline void next()
  {
    next(false);
  }

  void next(bool xtruncate);

  // EVALUATION ITERATOR FACET

  section_space_schema_member& schema_anchor();

  const section_space_schema_member& schema_anchor() const;

  void put_schema_anchor(const section_space_schema_member& xschema_anchor);

  void put_schema_anchor(pod_index_type xschema_anchor_index);

  void put_schema_anchor(const scoped_index& xschema_anchor_index);

  base_space_poset& base_space();

  const base_space_poset& base_space() const;

  int db() const;

  pod_index_type type_id() const;

  scoped_index client_id() const;

  section_evaluator& evaluator();

  const section_evaluator& evaluator() const;

  block<scoped_index>& discretization_members();

  const block<scoped_index>& discretization_members() const;

  scoped_index discretization_client_id(size_type xi);

  void gather_dofs(const sec_vd& xsec, block<sec_vd::dof_type>& xdofs);

  void scatter_dofs(sec_vd& xsec, const block<sec_vd::dof_type>& xdofs, size_type xindex);

protected:

  void initialize_schema_anchor(const section_space_schema_member& xanchor);


private:

  section_space_schema_member* _schema_anchor;

  base_space_poset* _base_space;

  bool _down_set_collection_enabled;

  eval_family* _evaluators;

  block<scoped_index> _discretization_members;

  block<scoped_index> _non_discretization_members;

  scoped_index _evaluator_id;

  scoped_index _maximal_eval_id;

  bool _vertex_discretization;

  bool _disc_reversal_enabled;

  index_space_handle* _eval_id_space;

  index_space_handle* _disc_id_space;

  subposet* _eval_sp;

  subposet* _disc_sp;

  subposet* _mono_sp;

  void put_schema_variables();

};

} // namespace fiber_bundle

#endif // ifndef EVAL_ITERATOR_H