solver_utils.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: solver_utils.h 2605 2014-08-15 03:36:44Z secanell $
//
//    Copyright (C) 2008 by Marc Secanell
//
//    This file is subject to QPL and may not be  distributed
//    without copyright and license information. Please refer
//    to the file deal.II/doc/license.html for the  text  and
//    further information on this license.
//
//---------------------------------------------------------------------------


#ifndef _SOLVER_UTILS__H
#define _SOLVER_UTILS__H

//-- deal.II 
#include <deal.II/lac/block_sparse_matrix.h>
#include <deal.II/lac/full_matrix.h>
#include <deal.II/lac/block_vector.h>
#include <deal.II/lac/petsc_parallel_block_sparse_matrix.h>

//-- OpenFCST
#include <application_core/application_wrapper.h>

// Include STD classes
#include <vector>
#include <fstream>
#include <iostream>
#include <sstream>

// Use namespace of deal.II
using namespace dealii;

namespace SolverUtils
{
    void check_diagonal(const BlockSparseMatrix<double>& A);

    void output_diagonal(const BlockSparseMatrix<double>& A);

    void print_diagonal(const BlockSparseMatrix<double>& A,
            const std::string& file = std::string("diag_matrix.dat"));
    void repair_diagonal (BlockSparseMatrix<double>& A);
#ifdef OPENFCST_WITH_PETSC
    void repair_diagonal (PETScWrappers::MPI::SparseMatrix & A);
#endif

    void repair_diagonal (BlockSparseMatrix<double>& A, FuelCell::ApplicationCore::FEVector& ,  const FuelCell::ApplicationCore::FEVector& );


};


#endif //_SOLVER_UTILS__H