solver_utils.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: solver_utils.h 1354 2013-08-17 00:01:22Z 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

// Include deal.II classes
#include <lac/block_sparse_matrix.h>
#include <lac/full_matrix.h>
#include "lac/block_vector.h"
#include "appframe/base.h"

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

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

class SolverUtils
{
 public:
  SolverUtils();
  
  void check_diagonal(const BlockSparseMatrix<double>& A) const;

  void output_diagonal(const BlockSparseMatrix<double>& A) const;
  
  void print_diagonal(const BlockSparseMatrix<double>& A, 
                      const std::string& file = std::string("diag_matrix.dat")) const;
  void repair_diagonal (BlockSparseMatrix<double>& A);

  void repair_diagonal (BlockSparseMatrix<double>& A, AppFrame::FEVector& ,  const AppFrame::FEVector& );
};


#endif //_SOLVER_UTILS__H