jcb_e33.cc

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


#include "SheafSystem/jcb_e33.impl.h"

#include "SheafSystem/abstract_poset_member.impl.h"
#include "SheafSystem/assert_contract.h"
#include "SheafSystem/at1_space.h"
#include "SheafSystem/e3.h"
#include "SheafSystem/jcb_space.h"
#include "SheafSystem/fiber_bundles_namespace.h"
#include "SheafSystem/schema_poset_member.h"
#include "SheafSystem/wsv_block.h"

using namespace std;
using namespace fiber_bundle; // Workaround for MS C++ bug.


//==============================================================================
// CLASS JCB_E33_LITE
//==============================================================================


//============================================================================
// JCB_E33 FACET OF CLASS JCB_E33_LITE
//============================================================================

fiber_bundle::jcb_e33_lite::
jcb_e33_lite()
{

  // Preconditions:

  // Body:

  // Postconditions:

  ensure(invariant());

  // Exit:
}

fiber_bundle::jcb_e33_lite::
jcb_e33_lite(const jcb_e33_lite& xother)
{
  // Preconditions:

  // Body:

  *this = xother;

  // Postconditions:

  ensure(invariant());

  // Exit:
}

fiber_bundle::jcb_e33_lite&
fiber_bundle::jcb_e33_lite::
operator=(const jcb_e33_lite& xother)
{

  // Preconditions:


  // Body:

  if(this == &xother)
    return *this;

  _row_dofs = xother._row_dofs;
  

  // Postconditions:

  ensure(invariant());
  ensure_for_all(i, 0, d(), component(i) == xother[i]); 

  // Exit:

  return *this;
}

fiber_bundle::jcb_e33_lite::
~jcb_e33_lite()
{
  // Preconditions:

  // Body:

  // Postconditions:

  // Exit:

}

fiber_bundle::jcb_e33_lite::
jcb_e33_lite(const row_dofs_type& xrow_dofs)
{
  // Preconditions:

  // Body:

  *this = xrow_dofs;

  // Postconditions:

  ensure(invariant());

  // Exit:
}

fiber_bundle::jcb_e33_lite::
jcb_e33_lite(const matrix_type& xmatrix)
{
  // Preconditions:

  // Body:

  *this = xmatrix;

  // Postconditions:

  ensure(invariant());

  // Exit:
}

fiber_bundle::jcb_e33_lite&
fiber_bundle::jcb_e33_lite::
operator=(const row_dofs_type& xrow_dofs)
{
  // Preconditions:

  // Body:

  _row_dofs = xrow_dofs;
  
  // Postconditions:

  ensure(invariant());

  // Exit:

  return *this;

}

fiber_bundle::jcb_e33_lite& 
fiber_bundle::jcb_e33_lite::
operator=(const matrix_type& xmatrix)
{
  // Preconditions:

  // Body:

  //_row_dofs = xmatrix;
  _row_dofs = reinterpret_cast<const row_dofs_type&>(xmatrix);

  // Postconditions:

  ensure(invariant());

  // Exit:

  return *this;

}

fiber_bundle::jcb_e33_lite::
operator jcb_e33_lite::row_dofs_type& ()
{
  // Preconditions:

  // Body:

  row_dofs_type& result = _row_dofs;

  // Postconditions:

  // Exit:

  return result;
}

fiber_bundle::jcb_e33_lite::
operator const jcb_e33_lite::row_dofs_type& () const
{
  // Preconditions:

  // Body:

  const row_dofs_type& result = _row_dofs;

  // Postconditions:

  // Exit:

  return result;
}

fiber_bundle::jcb_e33_lite::
operator jcb_e33_lite::matrix_type& ()
{
  // Preconditions:

  // Body:

  matrix_type& result = _row_dofs;

  // Postconditions:

  // Exit:

  return result;
}

fiber_bundle::jcb_e33_lite::
operator const jcb_e33_lite::matrix_type& () const
{
  // Preconditions:

  // Body:

  const matrix_type& result = _row_dofs;

  // Postconditions:

  // Exit:

  return result;
}


fiber_bundle::jcb_e33_lite::
jcb_e33_lite(
  const value_type& dxdu, const value_type& dxdv, const value_type& dxdw,
  const value_type& dydu, const value_type& dydv, const value_type& dydw,
  const value_type& dzdu, const value_type& dzdv, const value_type& dzdw)

{
  // Preconditions:

  // Body:

  put_components(dxdu, dxdv, dxdw, dydu, dydv, dydw,  dzdu, dzdv, dzdw);

  // Postconditions:

  ensure(invariant());
  ensure(isunordered_or_equals(component(0), dxdu));
  ensure(isunordered_or_equals(component(1), dxdv));
  ensure(isunordered_or_equals(component(2), dxdw));
  ensure(isunordered_or_equals(component(3), dydu));
  ensure(isunordered_or_equals(component(4), dydv));
  ensure(isunordered_or_equals(component(5), dydw));
  ensure(isunordered_or_equals(component(6), dzdu));
  ensure(isunordered_or_equals(component(7), dzdv));
  ensure(isunordered_or_equals(component(8), dzdw));

  // Exit:
}

void
fiber_bundle::jcb_e33_lite::
put_components(const value_type& dxdu,
               const value_type& dxdv,
               const value_type& dxdw,
               const value_type& dydu,
               const value_type& dydv,
               const value_type& dydw,
               const value_type& dzdu,
               const value_type& dzdv,
               const value_type& dzdw)
{
  // Preconditions:

  // Body:

  put_component(0, dxdu);
  put_component(1, dxdv);
  put_component(2, dxdw);
  put_component(3, dydu);
  put_component(4, dydv);
  put_component(5, dydw);
  put_component(6, dzdu);
  put_component(7, dzdv);
  put_component(8, dzdw);

  // Postconditions:

  ensure(invariant());

  ensure(isunordered_or_equals(component(0), dxdu));
  ensure(isunordered_or_equals(component(1), dxdv));
  ensure(isunordered_or_equals(component(2), dxdw));
  ensure(isunordered_or_equals(component(3), dydu));
  ensure(isunordered_or_equals(component(4), dydv));
  ensure(isunordered_or_equals(component(5), dydw));
  ensure(isunordered_or_equals(component(6), dzdu));
  ensure(isunordered_or_equals(component(7), dzdv));
  ensure(isunordered_or_equals(component(8), dzdw));

  // Exit:
}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// JCB_ED FACET OF CLASS JCB_E33_LITE
//============================================================================

// PUBLIC MEMBER FUNCTIONS

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// JACOBIAN (JCB) FACET OF CLASS JCB_E33_LITE
//============================================================================

// PUBLIC MEMBER FUNCTIONS

int
fiber_bundle::jcb_e33_lite::
dd() const
{
  // Preconditions:

  // Body:

  int result = 3;

  // Postconditions:

  ensure(invariant());
  ensure(result == 3);

  // Exit:

  return result;
}

int
fiber_bundle::jcb_e33_lite::
dr() const
{
  // Preconditions:

  // Body:

  int result = 3;

  // Postconditions:

  ensure(invariant());
  ensure(result == 3);

  // Exit:

  return result;
}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// VECTOR ALGEBRA (VD) FACET OF CLASS JCB_E33_LITE
//============================================================================

// PUBLIC MEMBER FUNCTIONS


// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// CARTESIAN ALGEBRA (TUPLE) FACET OF CLASS JCB_E33_LITE
//============================================================================

// PUBLIC MEMBER FUNCTIONS

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// ABSTRACT POSET MEMBER FACET OF CLASS JCB_E33_LITE
//============================================================================

// PUBLIC MEMBER FUNCTIONS

const std::string&
fiber_bundle::jcb_e33_lite::
class_name() const
{
  // Preconditions:

  // Body:

  const string& result = static_class_name();

  // Postconditions:

  ensure(!result.empty());

  // Exit:

  return result;
}

const std::string&
fiber_bundle::jcb_e33_lite::
static_class_name()
{
  // Preconditions:

  // Body:

  static const string result("jcb_e33_lite");

  // Postconditions:

  ensure(!result.empty());

  // Exit:

  return result;
}

fiber_bundle::jcb_e33_lite*
fiber_bundle::jcb_e33_lite::
clone() const
{
  jcb_e33_lite* result = 0;

  // Preconditions:

  // Body:

  result = new jcb_e33_lite();

  // Postconditions:

  ensure(result != 0);
  ensure(is_same_type(*result));

  // Exit:

  return result;
}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// ANY FACET OF CLASS JCB_E33_LITE
//============================================================================

// PUBLIC MEMBER FUNCTIONS

bool
fiber_bundle::jcb_e33_lite::
is_ancestor_of(const any_lite& xother) const
{
  // Preconditions:

  require(&xother != 0);

  // Body:

  // True if other conforms to this.

  bool result = dynamic_cast<const jcb_e33_lite*>(&xother) != 0;

  // Postconditions:

  return result;
}

bool
fiber_bundle::jcb_e33_lite::
invariant() const
{
  bool result = true;

  if(invariant_check())
  {
    // Prevent recursive calls to invariant.

    disable_invariant_check();

    // Must satisfy base class invariant.

    invariance(jcb_ed_lite::invariant());

    // Invariances for this class:

    // Finished, turn invariant checking back on.

    enable_invariant_check();
  }

  // Exit

  return result;
}

void*
fiber_bundle::jcb_e33_lite::
row_dofs()
{
  return &_row_dofs;
}

const void*
fiber_bundle::jcb_e33_lite::
row_dofs() const
{
  return &_row_dofs;
}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//==============================================================================
// CLASS JCB_E33
//==============================================================================

// ===========================================================
// HOST FACTORY FACET
// ===========================================================

// PUBLIC MEMBER FUNCTIONS

const sheaf::poset_path&
fiber_bundle::jcb_e33::
standard_schema_path()
{
  // Preconditions:


  // Body:

  static const poset_path result(standard_schema_poset_name(), "jcb_e33_schema");

  // Postconditions:

  // Exit:

  return result;
}

void
fiber_bundle::jcb_e33::
make_standard_schema(namespace_poset& xns)
{
  // Preconditions:

  require(xns.state_is_read_write_accessible());
  require(xns.contains_poset(standard_schema_poset_name()));
  require(!xns.contains_poset_member(standard_schema_path()));

  // Body:

  string lmember_names = "dxdu DOUBLE false ";
  lmember_names       += "dxdv DOUBLE false ";
  lmember_names       += "dxdw DOUBLE false ";
  lmember_names       += "dydu DOUBLE false ";
  lmember_names       += "dydv DOUBLE false ";
  lmember_names       += "dydw DOUBLE false ";
  lmember_names       += "dzdu DOUBLE false ";
  lmember_names       += "dzdv DOUBLE false ";
  lmember_names       += "dzdw DOUBLE false";

  schema_poset_member lschema(xns,
                              standard_schema_path().member_name(),
                              jcb_ed::standard_schema_path(),
                              lmember_names,
                              false);

  lschema.detach_from_state();

  // Postconditions:

  ensure(xns.contains_poset_member(standard_schema_path()));

  // Exit:

  return;
}

fiber_bundle::jcb_e33::host_type&
fiber_bundle::jcb_e33::
new_host(namespace_type& xns, 
         const poset_path& xhost_path, 
         const poset_path& xschema_path, 
     const poset_path& xdomain_path,
     const poset_path& xrange_path,
         bool xauto_access)
{
  // cout << endl << "Entering jcb_e33::new_host." << endl;

  // Preconditions:

  require(xns.state_is_auto_read_write_accessible(xauto_access));

  require(!xhost_path.empty());
  require(!xns.contains_path(xhost_path, xauto_access));

  require(xschema_path.full());
  require(xns.path_is_auto_read_accessible(xschema_path, xauto_access));
  require(schema_poset_member::conforms_to(xns, xschema_path, standard_schema_path()));
  require(schema_poset_member::row_dof_ct(xns, xschema_path, xauto_access) == 9);

  require(xns.path_is_auto_read_accessible(xdomain_path, xauto_access));
  require(xns.contains_poset<domain_type::host_type>(xdomain_path, xauto_access));
  require(xns.member_poset(xdomain_path, xauto_access).schema(xauto_access).conforms_to(domain_type::standard_schema_path()));
  require(xns.member_poset<domain_type::host_type>(xdomain_path, xauto_access).d(xauto_access) == 3);

  require(xns.path_is_auto_read_accessible(xrange_path, xauto_access));
  require(xns.contains_poset<range_type::host_type>(xrange_path, xauto_access));
  require(xns.member_poset(xrange_path, xauto_access).schema(xauto_access).conforms_to(range_type::standard_schema_path()));
  require(xns.member_poset<range_type::host_type>(xrange_path, xauto_access).d(xauto_access) == 3);

  require(host_type::d(xns, xschema_path, xauto_access) == host_type::d(xns, xdomain_path, xrange_path, xauto_access));

  require(xns.member_poset<domain_type::host_type>(xdomain_path, xauto_access).scalar_space_path(xauto_access) ==
      xns.member_poset<range_type::host_type>(xrange_path, xauto_access).scalar_space_path(xauto_access));

  // Body:

  host_type& result =
    host_type::new_table(xns, xhost_path, xschema_path, xdomain_path, xrange_path, xauto_access); 

  // Postconditions:

  ensure(xns.owns(result, xauto_access));
  ensure(result.path(true) == xhost_path);
  ensure(result.state_is_not_read_accessible());
  ensure(result.schema(true).path(xauto_access) == xschema_path);

  ensure(result.factor_ct(true) == 9);
  ensure(result.d(true) == 9);
  ensure(result.dd(true) == 3);
  ensure(result.dr(true) == 3);
  ensure(result.domain_path(true) == xdomain_path);
  ensure(result.range_path(true) == xrange_path);

  ensure(result.scalar_space_path(true) ==
     xns.member_poset<domain_type::host_type>(xdomain_path, xauto_access).scalar_space_path(xauto_access) );
  ensure(result.scalar_space_path(true) ==
     xns.member_poset<range_type::host_type>(xrange_path, xauto_access).scalar_space_path(xauto_access) );
  
  // Exit:

  // cout << "Leaving jcb_e33::new_host." << endl;
  return result;
}

fiber_bundle::jcb_e33::host_type&
fiber_bundle::jcb_e33::
standard_host(namespace_type& xns, const std::string& xsuffix, bool xauto_access)
{
  // cout << endl << "Entering jcb_e33::new_host." << endl;

  // Preconditions:

  require(xns.state_is_auto_read_write_accessible(xauto_access));

  require(xsuffix.empty() || poset_path::is_valid_name(xsuffix));
  require(standard_host_is_available<jcb_e33>(xns, xsuffix, xauto_access));

  require(xns.path_is_auto_read_accessible(standard_schema_path(), xauto_access));
  
  require(xns.path_is_auto_read_available(standard_host_path<domain_type>(xsuffix), xauto_access));
  require(xns.path_is_auto_read_available(standard_host_path<range_type>(xsuffix), xauto_access));

  // Body:

  // Create the domain space if necessary.

  poset_path ldomain_space_path = domain_type::standard_host(xns, xsuffix, xauto_access).path(true);

  // Create the range space if necessary.

  poset_path lrange_space_path = range_type::standard_host(xns, xsuffix, xauto_access).path(true);

  poset_path lpath(standard_host_path<jcb_e33>(xsuffix));

  host_type* result_ptr;
  if(xns.contains_path(lpath, xauto_access))
  {
    result_ptr = &xns.member_poset<host_type>(lpath, xauto_access);
  }
  else
  {
    result_ptr = &new_host(xns, lpath, standard_schema_path(), ldomain_space_path, lrange_space_path, xauto_access);
  }

  host_type& result = *result_ptr;

  // Postconditions:

  ensure(xns.owns(result, xauto_access));
  ensure(result.path(true) == standard_host_path<jcb_e33>(xsuffix));
  ensure(result.state_is_not_read_accessible());
  ensure(result.schema(true).path(xauto_access) == standard_schema_path());

  ensure(result.factor_ct(true) == 9);
  ensure(result.d(true) == 9);
  ensure(result.dd(true) == 3);
  ensure(result.dr(true) == 3);
  ensure(result.domain_path(true) == standard_host_path<domain_type>(xsuffix));
  ensure(result.range_path(true) == standard_host_path<range_type>(xsuffix));

  ensure(result.scalar_space_path(true) == standard_host_path<domain_type::scalar_type>(xsuffix));
  ensure(result.scalar_space_path(true) == standard_host_path<domain_type::scalar_type>(xsuffix));

  // Exit:

  // cout << "Leaving jcb_e33::new_host." << endl;
  return result;
}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS
 

//============================================================================
// JCB_E33 FACET OF CLASS JCB_E33
//============================================================================

// PUBLIC MEMBER FUNCTIONS

fiber_bundle::jcb_e33::
jcb_e33()
{

  // Preconditions:

  // Body:

  // Postconditions:

  ensure(invariant());
}

fiber_bundle::jcb_e33::
jcb_e33(const poset_state_handle* xhost, pod_index_type xhub_id)
{
  // Preconditions:

  require(xhost != 0);
  require(xhost->state_is_read_accessible());
  require(xhost->contains_member(xhub_id));

  // Body:

  attach_to_state(xhost, xhub_id);

  // Postconditions:

  ensure(invariant());
  // ensure(host() == xhost);
  ensure(index() == xhub_id);
  ensure(is_attached());
}

fiber_bundle::jcb_e33::
jcb_e33(const poset_state_handle* xhost, const scoped_index& xid)
{
  // Preconditions:

  require(xhost != 0);
  require(xhost->state_is_read_accessible());
  require(xhost->contains_member(xid));

  // Body:

  attach_to_state(xhost, xid.hub_pod());

  // Postconditions:

  ensure(invariant());
  // ensure(host() == xhost);
  ensure(index() ==~ xid);
  ensure(is_attached());
}

fiber_bundle::jcb_e33::
jcb_e33(const poset_state_handle* xhost, const std::string& xname)
{

  // Preconditions:

  require(xhost != 0);
  require(xhost->state_is_read_accessible());
  require(!xname.empty());
  require(xhost->contains_member(xname));

  // Body:

  attach_to_state(xhost, xname);

  // Postconditions:

  ensure(invariant());
  // ensure(host() == xhost);
  ensure(name() == xname);
  ensure(is_attached());

}

fiber_bundle::jcb_e33::
jcb_e33(abstract_poset_member* xother)
{

  // Preconditions:

  require(xother != 0);

  // Body:

  attach_to_state(xother);

  // Postconditions:

  ensure(invariant());
  ensure(is_attached());
  ensure(is_same_state(xother));

}

fiber_bundle::jcb_e33::
jcb_e33(poset_state_handle* xhost, bool xauto_access)
{

  // Preconditions:

  require(precondition_of(new_jim_state(xhost, 0, false, xauto_access)));

  // Body:

  new_jim_state(xhost, 0, false, xauto_access);

  // Postconditions:

  ensure(postcondition_of(new_jim_state(xhost, 0, false, xauto_access)));

  // Exit:

  return;
}

fiber_bundle::jcb_e33::
jcb_e33(poset_state_handle& xhost,
       const row_dofs_type& xrdt,
       bool xauto_access)
{

  // Preconditions:

  require(precondition_of(new_jim_state(&xhost, 0, false, xauto_access)));

  // Body:

  new_jim_state(&xhost, 0, false, xauto_access);

  if(xauto_access)
  {
    xhost.get_read_write_access();
  }

  *this = xrdt;

  if(xauto_access)
  {
    xhost.release_access();
  }

  // Postconditions:

  ensure(postcondition_of(new_jim_state(&xhost, 0, false, xauto_access)));

  // Exit:

  return;
}

fiber_bundle::jcb_e33&
fiber_bundle::jcb_e33::
operator=(const row_dofs_type& xrdt)
{
  // Preconditions:

  require(state_is_read_write_accessible());

  // Body:

  int ld = d();
  for(int i=0; i<ld; ++i)
  {
    put_component(i, xrdt.components[i]);
  }

  // Postconditions:

  ensure_for_all(i, 0, d(), component(i) == xrdt.components[i]);

  // Exit:

  return *this;

}

fiber_bundle::jcb_e33&
fiber_bundle::jcb_e33::
operator=(const abstract_poset_member& xother)
{
  // Preconditions:

  require(is_ancestor_of(&xother));
  require(precondition_of(attach_to_state(&xother)));

  // Body:

  attach_to_state(&xother);

  // Postconditions:

  ensure(postcondition_of(attach_to_state(&xother)));

  // Exit:

  return *this;
}

fiber_bundle::jcb_e33&
fiber_bundle::jcb_e33::
operator=(const jcb_e33& xother)
{
  // Preconditions:

  require(precondition_of(attach_to_state(&xother)));

  // Body:

  attach_to_state(&xother);

  // Postconditions:

  ensure(postcondition_of(attach_to_state(&xother)));

  // Exit:

  return *this;
}

fiber_bundle::jcb_e33::
~jcb_e33()
{

  // Preconditions:

  // Body:

  // Postconditions:

}

const fiber_bundle::jcb_e33::volatile_type&
fiber_bundle::jcb_e33::
lite_prototype() const
{
  // Preconditions:

  // Body:

  static const volatile_type result;

  // Postconditions:

  // Exit:

  return result;
}

fiber_bundle::jcb_e33::volatile_type*
fiber_bundle::jcb_e33::
lite_type() const
{
  // Preconditions:

  // Body:

  volatile_type* result = new volatile_type(sheaf::row_dofs(*this));

  // Postconditions:

  // Exit:

  return result;
}


//============================================================================
// ROW DOF RELATED MEMBER FUNCTIONS
//============================================================================

fiber_bundle::jcb_e33::
operator jcb_e33::row_dofs_type& ()
{
  // Preconditions:

  // Body:

  row_dofs_type& result = sheaf::row_dofs(*this);

  // Postconditions:

  // Exit:

  return result;
}

fiber_bundle::jcb_e33::
operator const jcb_e33::row_dofs_type& () const
{
  // Preconditions:

  // Body:

  const row_dofs_type& result = sheaf::row_dofs(*this);

  // Postconditions:

  // Exit:

  return result;
}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// JCB_ED FACET OF CLASS JCB_E33
//============================================================================

// PUBLIC MEMBER FUNCTIONS

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// JCB FACET OF CLASS JCB_E33
//============================================================================

// PUBLIC MEMBER FUNCTIONS

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// VD FACET OF CLASS JCB_E33
//============================================================================

// PUBLIC MEMBER FUNCTIONS

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// TUPLE FACET OF CLASS JCB_E33
//============================================================================

// PUBLIC MEMBER FUNCTIONS

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// ABSTRACT POSET MEMBER FACET OF CLASS JCB_E33
//============================================================================

// PUBLIC MEMBER FUNCTIONS

const std::string&
fiber_bundle::jcb_e33::
class_name() const
{
  // Preconditions:

  // Body:

  const string& result = static_class_name();

  // Postconditions:

  ensure(!result.empty());

  // Exit:

  return result;
}

const std::string&
fiber_bundle::jcb_e33::
static_class_name()
{
  // Preconditions:

  // Body:

  static const string result("jcb_e33");

  // Postconditions:

  ensure(!result.empty());

  // Exit:

  return result;
}

fiber_bundle::jcb_e33*
fiber_bundle::jcb_e33::
clone() const
{

  // Preconditions:

  // Body:

  // Create new handle of the current class.

  jcb_e33* result = new jcb_e33();

  // Postconditions:

  ensure(result != 0);
  ensure(result->invariant());

  // Exit:

  return result;

}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//============================================================================
// ANY FACET OF CLASS JCB_E33
//============================================================================

// PUBLIC MEMBER FUNCTIONS

bool
fiber_bundle::jcb_e33::
is_ancestor_of(const any* xother) const
{

  // Preconditions:

  require(xother != 0);

  // Body:

  // If other may be dynamically cast to the type of this then this is an
  // ancestor of other.

  bool result = dynamic_cast<const jcb_e33*>(xother) != 0;

  // Postconditions:

  //  ensure(invariant());

  // Exit:

  return result;

}

bool
fiber_bundle::jcb_e33::
invariant() const
{
  bool result = true;

  // Preconditions:

  // Body:

  // Must satisfy base class invariant.

  invariance(jcb_ed::invariant());

  if (invariant_check())
  {
    // Prevent recursive calls to invariant.

    disable_invariant_check();

    invariance(is_attached() ? (dr() == 3) : true);
    invariance(is_attached() ? (dd() == 3) : true);

    // Finished, turn invariant checking back on.

    enable_invariant_check();
  }

  // Postconditions:

  ensure(is_derived_query);

  // Exit:

  return result;
}

// PROTECTED MEMBER FUNCTIONS

// PRIVATE MEMBER FUNCTIONS


//==============================================================================
// NON-MEMBER FUNCTIONS
//==============================================================================

void
fiber_bundle::jcb_algebra::
push(const jcb_e33& xjcb, const e3& xvector, e3& xresult, bool xauto_access)
{
  // Preconditions:

  require(xjcb.state_is_auto_read_accessible(xauto_access));
  require(xvector.state_is_auto_read_accessible(xauto_access));
  require(xresult.state_is_auto_read_write_accessible(xauto_access));
  require(xvector.d(xauto_access) == xjcb.dd(xauto_access));
  require(xresult.d(xauto_access) == xjcb.dr(xauto_access));


  require(xvector.is_contravariant(0, xauto_access)); // !is_dual

  // Body:

  if(xauto_access)
  {
    xjcb.get_read_access();
    xvector.get_read_access();
    xresult.get_read_write_access(true);
  }

  // Left matrix multiply vector with jacobian.

  int num_rows = xjcb.dr();
  int num_cols = xjcb.dd();

  for(int i=0; i<num_rows; ++i)
  {
    vd_value_type sumi = 0.0;
    for(int j=0; j<num_cols; ++j)
    {
      vd_value_type vj = xvector.component(j);

      //@todo Implement xjcb.component(i, j, false) and/or
      //      index_for_row_column(i,j) as in jcb_lite.

      vd_value_type cij = xjcb.component(i, j);

      //int lindex = i*num_cols + j;
      //vd_value_type cij = xjcb.component(lindex);

      sumi += cij*vj;
    }
    xresult.put_component(i, sumi);
  }   

  // Make sure the result is contravariant/p_vector.

  xresult.put_is_p_vector(false);

  if(xauto_access)
  {
    xjcb.release_access();
    xvector.release_access();
    xresult.release_access();
  }

  // Postconditions:

  ensure(xresult.is_contravariant(0, xauto_access)); // !is_dual

  // Exit:
}

void
fiber_bundle::jcb_algebra::
push(const jcb_e33_lite& xjcb, const e3_lite& xvector, e3_lite& xresult)
{
  // Preconditions:

  require(xvector.d() == xjcb.dd());
  require(xresult.d() == xjcb.dr());


  require(unexecutable("xvector.is_contravariant())")); // !is_dual
  require(unexecutable("xresult.is_contravariant())")); // !is_dual

  // Body:

  // Left matrix multiply vector with jacobian.

  int num_rows = xjcb.dr();
  int num_cols = xjcb.dd();

  for(int i=0; i<num_rows; ++i)
  {
    vd_value_type sumi = 0.0;
    for(int j=0; j<num_cols; ++j)
    {
      vd_value_type vj = xvector.component(j);
      vd_value_type cij = xjcb.component(i, j);
      sumi += cij*vj;
    }
    xresult.put_component(i, sumi);
  }    

  // Postconditions:

  //ensure(unexecutable("???"));

}

fiber_bundle::e3_lite*
fiber_bundle::jcb_algebra::
push(const jcb_e33_lite& xjcb, const e3_lite& xvector)
{
  // Preconditions:

  require(precondition_of(push(xjcb, xvector, *result)));

  // Body:

  e3_lite* result = new e3_lite();

  push(xjcb, xvector, *result);

  // Postconditions:

  ensure(postcondition_of(push(xjcb, xvector, *result)));

  // Exit:

  return result;

}

void
fiber_bundle::jcb_algebra::
pull(const jcb_e33& xjcb, const e3& xcovector, e3& xresult, bool xauto_access)
{
  // Preconditions:

  require(xjcb.state_is_auto_read_accessible(xauto_access));
  require(xcovector.state_is_auto_read_accessible(xauto_access));
  require(xresult.state_is_auto_read_write_accessible(xauto_access));
  require(xcovector.d(xauto_access) == xjcb.dr(xauto_access));
  require(xresult.d(xauto_access) == xjcb.dd(xauto_access));


  require(xcovector.is_covariant(0, xauto_access)); // is_dual

  // Body:

  if(xauto_access)
  {
    xjcb.get_read_access();
    xcovector.get_read_access();
    xresult.get_read_write_access(true);
  }

  // Right matrix multiply covector with jacobian.

  int num_rows = xjcb.dd();
  int num_cols = xjcb.dr();

  for(int i=0; i<num_rows; ++i)
  {
    vd_value_type sumi = 0.0;
    for(int j=0; j<num_cols; ++j)
    {
      vd_value_type vj = xcovector.component(j);

      //@todo Implement xjcb.component(i, j, false) and/or
      //      index_for_row_column(i,j) as in jcb_lite.

      vd_value_type cji = xjcb.component(j, i);

      //int lindex = j*num_rows + i;
      //vd_value_type cji = xjcb.component(lindex);

      sumi += vj*cji;
    }
    xresult.put_component(i, sumi);
  }    

  // Make sure the result is covariant/p_form.

  xresult.put_is_p_form(false);

  if(xauto_access)
  {
    xjcb.release_access();
    xcovector.release_access();
    xresult.release_access();
  }

  // Postconditions:

  //ensure(unexecutable("???"));
  ensure(xresult.is_covariant(0, xauto_access));   // is_dual

  // Exit:
}

void
fiber_bundle::jcb_algebra::
pull(const jcb_e33_lite& xjcb, const e3_lite& xcovector, e3_lite& xresult)
{
  // Preconditions:

  require(xcovector.d() == xjcb.dr());
  require(xresult.d() == xjcb.dd());


  require(unexecutable("xcovector.is_covariant())")); // is_dual
  require(unexecutable("xresult.is_covariant())"));   // is_dual

  // Body:

  // Right matrix multiply covector with jacobian.

  int num_rows = xjcb.dd();
  int num_cols = xjcb.dr();

  for(int i=0; i<num_rows; ++i)
  {
    vd_value_type sumi = 0.0;
    for(int j=0; j<num_cols; ++j)
    {
      vd_value_type vj = xcovector.component(j);
      vd_value_type cji = xjcb.component(j, i);
      sumi += vj*cji;
    }
    xresult.put_component(i, sumi);
  }    

  // Postconditions:

  //ensure(unexecutable("???"));

  // Exit:
}

fiber_bundle::e3_lite*
fiber_bundle::jcb_algebra::
pull(const jcb_e33_lite& xjcb, const e3_lite& xcovector)
{
  // Preconditions:

  require(precondition_of(pull(xjcb, xcovector, *result)));

  // Body:

  e3_lite* result = new e3_lite();

  pull(xjcb, xcovector, *result);

  // Postconditions:

  ensure(postcondition_of(pull(xjcb, xcovector, *result)));

  // Exit:

  return result;

}


//==============================================================================
// TEMPLATE SPECIALIZATIONS
//==============================================================================