packages/Alea/mc_solvers/GmresPolynomial.cc
//---------------------------------------------------------------------------//
// Build FOM coefficients from Arnoldi decomposition
//---------------------------------------------------------------------------//
Build FOM coefficients from Arnoldi decomposition
Build FOM coefficients from Arnoldi decomposition
void GmresPolynomial::buildFomPolyFromArnoldi()
{
    if( b_verbosity >= LOW )
        std::cout << "Creating FOM polynomial coefficients"
            << " of order " << b_m << " from Arnoldi process" << std::endl;
    Teuchos::RCP<Teuchos::SerialDenseMatrix<LO,SCALAR> > arnoldi_mat =
        computeArnoldiMatrix(b_m+1);
    CHECK( arnoldi_mat != Teuchos::null );
    Teuchos::ArrayRCP<const Anasazi::Value<SCALAR> > fom_roots =
        computeFomRoots(arnoldi_mat);
    CHECK( fom_roots.size() == b_m+1 );
    b_coeffs = computeCoeffsFromRoots(fom_roots);
}
void GmresPolynomial :: buildFomPolyFromArnoldi ( ) { if ( b_verbosity >= LOW ) std :: cout << " " << " " << b_m << " " << std :: endl ; Teuchos :: RCP < Teuchos :: SerialDenseMatrix < LO , SCALAR > > arnoldi_mat = computeArnoldiMatrix ( b_m + 1 ) ; CHECK ( arnoldi_mat != Teuchos :: null ) ; Teuchos :: ArrayRCP < const Anasazi :: Value < SCALAR > > fom_roots = computeFomRoots ( arnoldi_mat ) ; CHECK ( fom_roots . size ( ) == b_m + 1 ) ; b_coeffs = computeCoeffsFromRoots ( fom_roots ) ; }
Build FOM coefficients from Arnoldi decomposition
Build FOM coefficients from Arnoldi decomposition