test_Factory.cpp

Go to the documentation of this file.

/*--------------------------------------------------------------------*/
/*    Copyright 2005 Sandia Corporation.                              */
/*    Under the 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.                    */
/*--------------------------------------------------------------------*/
 
#include <fei_macros.hpp>
#include <fei_mpi.h>
 
#include <test_utils/test_Factory.hpp>
 
#include <test_utils/LibraryFactory.hpp>
 
#include <test_utils/test_VectorSpace.hpp>
#include <test_utils/test_MatrixGraph.hpp>
 
#include <fei_Factory.hpp>
#include <snl_fei_Factory.hpp>
 
#include <test_utils/test_Factory_helper.hpp>
 
#ifdef HAVE_FEI_AZTECOO
#include <fei_Aztec_LinSysCore.hpp>
#endif
 
#include <fei_Factory_Trilinos.hpp>
 
#undef fei_file
#define fei_file "test_Factory.cpp"
#include <fei_ErrMacros.hpp>
 
test_Factory::test_Factory(MPI_Comm comm)
  : tester(comm)
{
}
 
test_Factory::~test_Factory()
{
}
 
int test_Factory::runtests()
{
  {
  if (localProc_==0) FEI_COUT << "constructing Factory_Trilinos...";
 
  fei::SharedPtr<fei::Factory> factory(new Factory_Trilinos(comm_)); 
 
  if (localProc_==0) FEI_COUT << "ok" << FEI_ENDL;
 
  factory_test1(factory);
 
  if (localProc_==0) FEI_COUT << "testing fei::Factory::clone..." << FEI_ENDL;
 
  fei::SharedPtr<fei::Factory> clone(factory->clone());
 
  factory_test1(clone);
 
  FEI_COUT << FEI_ENDL;
  }
 
#ifdef HAVE_FEI_AZTECOO
  {
  if (localProc_==0) FEI_COUT << "constructing snl_fei::Factory(Aztec)...";
 
  fei::SharedPtr<LinearSystemCore> az_lsc(new fei_trilinos::Aztec_LinSysCore(comm_));
 
  fei::SharedPtr<fei::Factory> factory(new snl_fei::Factory(comm_, az_lsc));
 
  if (localProc_==0) FEI_COUT << "ok" << FEI_ENDL;
 
  factory_test1(factory);
 
  if (localProc_==0) FEI_COUT << "testing fei::Factory::clone..." << FEI_ENDL;
 
  fei::SharedPtr<fei::Factory> clone(factory->clone());
 
  factory_test1(clone);
 
  FEI_COUT << FEI_ENDL;
  }
 
#endif
  return(0);
}
 
void test_Factory::factory_test1(fei::SharedPtr<fei::Factory> factory)
{
  if (localProc_==0) FEI_COUT << "  testing factory->createVectorSpace...";
 
  fei::SharedPtr<fei::VectorSpace> vecspace =
    factory->createVectorSpace(comm_, "dummy_Name");
 
  if (vecspace.get() == 0) {
    FEI_COUT << "no"<<FEI_ENDL;
    throw std::runtime_error("factory failed to create a fei::VectorSpace");
  }
 
  //do an extremely simple test to make sure the vector-space
  //is 'alive'.
  int fieldID = 0;
  int fieldSize = 3;
  vecspace->defineFields(1, &fieldID, &fieldSize);
 
  if (vecspace->getNumFields() != 1) {
    FEI_COUT << "no"<<FEI_ENDL;
    throw std::runtime_error("vecspace->defineFields/getNumFields failed.");
  }
 
  if (localProc_==0) FEI_COUT << "ok"<<FEI_ENDL;
 
  if (localProc_==0) FEI_COUT << "  testing factory->createFEI...";
 
  fei::SharedPtr<FEI> fei = factory->createFEI(comm_);
 
  //again, do a simple test to make sure the FEI instance is alive...
 
  int err = fei->initFields(1, &fieldSize, &fieldID);
  if (err != 0) {
    FEI_COUT << "failed"<<FEI_ENDL;
    throw std::runtime_error("fei->initFields() failed.");
  }
 
  int testFieldSize = -1;
  err = fei->getFieldSize(fieldID, testFieldSize);
  if (err != 0 || testFieldSize != fieldSize) {
    FEI_COUT << "failed"<<FEI_ENDL;
    throw std::runtime_error("fei->getFieldSize() failed.");
  }
 
  if (localProc_==0) FEI_COUT << "ok"<<FEI_ENDL;
}