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

#ifndef LINEAR_3D_H
#define LINEAR_3D_H

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

#ifndef LINEAR_FCN_SPACE_H
#include "SheafSystem/linear_fcn_space.h"
#endif

namespace fiber_bundle
{

class SHEAF_DLL_SPEC linear_3d : public linear_fcn_space
{

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

public:

  linear_3d();

  linear_3d(const linear_3d& xother);

  virtual ~linear_3d();


  //   ///
  //   ///
  //   double inverse_jacobian(const double xjacobian[3][3],
  //                           double xinverse_jacobian[3][3]);

  //   ///
  //   ///
  //   double determinant(const dof_type xcoord_dofs[],
  //                      size_type xcoord_dofs_ub,
  //                      const coord_type xlocal_coords[],
  //                      size_type xlocal_coords_ub);

protected:

  double determinant_3x3(double a00, double a01, double a02,
                         double a10, double a11, double a12,
                         double a20, double a21, double a22);

private:

  enum static_const_int {DL = 4};  

  value_type _basis_value_buffer[DL];

  value_type _basis_deriv_value_buffer[3*DL];

  //value_type _jacobian_value_buffer[4];
  value_type _jacobian_value_buffer[128]; 

  //block<int[2]> _jacobian_value_block;
  //vector<int[2]> _jacobian_value_vector;


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

public:

protected:

private:


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

public:

  virtual int dl() const;

  virtual void basis_at_coord(const dof_type xlocal_coord[], size_type xlocal_coord_ub);

  void basis_derivs_at_coord(const dof_type xlocal_coords[],
                             size_type xlocal_coords_ub);

protected:

private:


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

public:

  value_type jacobian_determinant(const dof_type xcoord_dofs[],
                                  size_type xcoord_dofs_ub,
                                  size_type xdf,
                                  const coord_type xlocal_coords[],
                                  size_type xlocal_coords_ub);

  value_type volume(const dof_type xcoord_dofs[],
                    size_type xcoord_dofs_ub, size_type xdf);

  virtual void integrate(const dof_type xcoord_dofs[],
                         size_type xcoord_dofs_ub,
                         size_type xdf,
                         const dof_type xintegrands[],
                         size_type xintegrands_ub,
                         value_type xresult_integrals[],
                         size_type xresult_integrals_ub);

  virtual void integrate(const dof_type xcoord_dofs[],
                         size_type xcoord_dofs_ub,
                         size_type xdf,
                         const dof_type xintegrand,
                         value_type xresult_integrals[],
                         size_type xresult_integrals_ub);

protected:

  void gauss_point(pod_index_type xindex,
                   coord_type xresult[],
                   size_type xresult_ub);

private:


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

public:

  virtual void dxi_local(size_type xcrd,
                         const dof_type xsource_dofs[],
                         size_type xsource_dofs_ub,
                         dof_type xresult_dofs[],
                         size_type xresult_dofs_ub) const;

  void jacobian(const dof_type xcoord_dofs[],
                size_type xcoord_dofs_ub,
                size_type xdf,
                const dof_type xlocal_coords[],
                size_type xlocal_coords_ub);

  //   ///
  //   ///
  //   double jacobian(const dof_type xcoord_dofs[],
  //                   size_type xcoord_dofs_ub,
  //                   const coord_type xlocal_coords[],
  //                   size_type xlocal_coords_ub,
  //                   double xjacobian[3][3] = NULL);

protected:

private:


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

public:

  virtual int db() const;

  virtual void local_coordinates(pod_index_type xindex,
                                 coord_type xresult[],
                                 size_type xresult_ub) const;
  virtual void center(coord_type xresult[], size_type xresult_ub) const;

  virtual bool in_standard_domain(const dof_type xlocal_coords[],
                                  size_type xlocal_coords_ub) const;

protected:

private:


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

public:

  virtual void coord_at_value(const dof_type xdofs[],
                              size_type xdofs_ub,
                              const dof_type xglobal_coords[],
                              size_type xglobal_coord_ub,
                              dof_type xlocal_coords[],
                              size_type xlocal_coords_ub) const;

protected:

private:


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

public:

  virtual linear_3d* clone() const;

  virtual linear_3d& operator=(const section_evaluator& xother);

  linear_3d& operator=(const linear_3d& xother);

  virtual bool invariant() const;

  virtual bool is_ancestor_of(const any* xother) const;

protected:

private:


};

} // namespace fiber_bundle

#endif // ifndef LINEAR_3D_H