LinsysFunctions.hpp

/*------------------------------------------------------------------------*/
/*                 Copyright 2010 Sandia Corporation.                     */
/*  Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive   */
/*  license for use of this work by or on behalf of the U.S. Government.  */
/*  Export of this program may require a license from the                 */
/*  United States Government.                                             */
/*------------------------------------------------------------------------*/

#ifndef stk_linsys_LinsysFunctions_hpp
#define stk_linsys_LinsysFunctions_hpp

#include <stk_linsys/LinearSystemInterface.hpp>

#include <stk_mesh/base/Field.hpp>
#include <stk_mesh/base/Selector.hpp>
#include <stk_mesh/base/BulkData.hpp>

#include <Teuchos_ParameterList.hpp>


namespace stk_classic {
namespace linsys {

void add_connectivities(stk_classic::linsys::LinearSystemInterface& ls,
                        stk_classic::mesh::EntityRank from_type,
                        stk_classic::mesh::EntityRank to_connected_type,
                        const stk_classic::mesh::FieldBase& field,
                        const stk_classic::mesh::Selector& selector,
                        const stk_classic::mesh::BulkData& mesh_bulk);

void dirichlet_bc(stk_classic::linsys::LinearSystemInterface& ls,
                  const stk_classic::mesh::BulkData& mesh,
                  const stk_classic::mesh::Part& bcpart,
                  stk_classic::mesh::EntityRank entity_rank,
                  const stk_classic::mesh::FieldBase& field,
                  unsigned field_component,
                  double prescribed_value);


int fei_solve(int & status, fei::LinearSystem &fei_ls, const Teuchos::ParameterList & params);

double compute_residual_norm2(fei::LinearSystem& fei_ls, fei::Vector& r);

void copy_vector_to_mesh( fei::Vector & vec,
                          const DofMapper & dof,
                          stk_classic::mesh::BulkData & mesh_bulk_data
                        );

void scale_matrix(double scalar, fei::Matrix& matrix);

void add_matrix_to_matrix(double scalar,
                          const fei::Matrix& src_matrix,
                          fei::Matrix& dest_matrix);

void scale_vector(double scalar,
                  fei::Vector& vec);

void add_vector_to_vector(double scalar,
                          const fei::Vector& src_vector,
                          fei::Vector& dest_vector);

}//namespace linsys
}//namespace stk_classic

#endif