Ejemplos de RCP::GetProblem en C++ (Cpp)

Ejemplo n.º 1

Archivo: BlockKrylovSchurEpetraExBuckling.cpp Proyecto: Tech-XCorp/Trilinos

AmesosBucklingOp::AmesosBucklingOp( Epetra_LinearProblem& problem,
				    const Teuchos::RCP<Amesos_BaseSolver>& solver,
				    const Teuchos::RCP<Epetra_Operator>& stiffMtx
				    )
  : solver_(solver),
    stiffMtx_(stiffMtx)
{
  problem_ = const_cast<Epetra_LinearProblem*>( solver->GetProblem() );
}

Ejemplo n.º 2

Archivo: N_LAS_MOROperators.C Proyecto: Xycedev/Xyce

//-----------------------------------------------------------------------------
// Function      : N_LAS_AmesosGenOp::N_LAS_AmesosGenOp
// Purpose       : Constructor 
// Special Notes :
// Scope         : Public
// Creator       : Heidi Thornquist, SNL 
// Creation Date : 06/04/12
//-----------------------------------------------------------------------------
N_LAS_AmesosGenOp::N_LAS_AmesosGenOp( const Teuchos::RCP<Amesos_BaseSolver>& solver,
                                      const Teuchos::RCP<Epetra_Operator>& B,
                                      bool useTranspose )
  : useTranspose_(useTranspose),
    solver_(solver),
    B_(B)
{
  problem_ = Teuchos::rcp( const_cast<Epetra_LinearProblem*>( solver->GetProblem() ), false);

  if ( solver_->UseTranspose() )
    solver_->SetUseTranspose(!useTranspose);
  else
    solver_->SetUseTranspose(useTranspose);

  if ( B_->UseTranspose() )
    B_->SetUseTranspose(!useTranspose);
  else
    B_->SetUseTranspose(useTranspose);
}

Ejemplo n.º 3

Archivo: BlockKrylovSchurEpetraExGenAmesos.cpp Proyecto: cakeisalie/oomphlib_003

AmesosGenOp::AmesosGenOp( Epetra_LinearProblem& problem,
                          const Teuchos::RCP<Amesos_BaseSolver>& solver,
                          const Teuchos::RCP<Epetra_Operator>& massMtx,
                          bool useTranspose )
  : useTranspose_(useTranspose),
    solver_(solver),
    massMtx_(massMtx)
{
  problem_ = const_cast<Epetra_LinearProblem*>( solver->GetProblem() );
  
  if ( solver_->UseTranspose() )
    solver_->SetUseTranspose(!useTranspose); 
  else
    solver_->SetUseTranspose(useTranspose);
  
  if ( massMtx_->UseTranspose() )
    massMtx_->SetUseTranspose(!useTranspose);
  else
    massMtx_->SetUseTranspose(useTranspose);    
}

Ejemplo n.º 4

Archivo: Anasazi_Block_KrylovSchur_Amesos.cpp Proyecto: 00liujj/trilinos

AmesosGenOp::
AmesosGenOp (const Teuchos::RCP<Amesos_BaseSolver>& solver,
             const Teuchos::RCP<Epetra_Operator>& massMtx,
             const bool useTranspose)
  : solver_ (solver),
    massMtx_ (massMtx),
    problem_ (NULL),
    useTranspose_ (useTranspose)
{
  if (solver.is_null ()) {
    throw std::invalid_argument ("AmesosGenOp constructor: The 'solver' "
                                 "input argument is null.");
  }
  if (massMtx.is_null ()) {
    throw std::invalid_argument ("AmesosGenOp constructor: The 'massMtx' "
                                 "input argument is null.");
  }

  Epetra_LinearProblem* problem = const_cast<Epetra_LinearProblem*> (solver->GetProblem ());
  if (problem == NULL) {
    throw std::invalid_argument ("The solver's GetProblem() method returned "
                                 "NULL.  This probably means that its "
                                 "LinearProblem has not yet been set.");
  }
  problem_ = problem;

  if (solver_->UseTranspose ()) {
    solver_->SetUseTranspose (! useTranspose);
  } else {
    solver_->SetUseTranspose (useTranspose);
  }

  if (massMtx_->UseTranspose ()) {
    massMtx_->SetUseTranspose (! useTranspose);
  } else {
    massMtx_->SetUseTranspose (useTranspose);
  }
}