FEI_tester.hpp

/*--------------------------------------------------------------------*/
/*    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.                    */
/*--------------------------------------------------------------------*/
 
#ifndef _FEI_tester_h_
#define _FEI_tester_h_
 
#include <fei_macros.hpp>
#include <fei_SharedPtr.hpp>
#include <fei_defs.h>
 
class LibraryWrapper;
class LinearSystemCore;
class FiniteElementData;
 
#include <test_utils/feitester.hpp>
#include <test_utils/DataReader.hpp>
 
#include <FEI_Implementation.hpp>
 
class FEI_tester : public feitester {
 public:
  FEI_tester(fei::SharedPtr<DataReader> data_reader,
       MPI_Comm comm, int localProc, int numProcs, bool useNewFEI=false);
  ~FEI_tester();
 
  const char* getName()
    {
      static const char name[] = "FEI_tester";
      return((const char*)name);
    }
 
  int testInitialization();
 
  int testLoading();
 
  int testSolve();
 
  int testCheckResult();
 
  void dumpMatrixFiles();
 
  void setParameter(const char* param);
 
 private:
  int createFEIinstance(const char* solverName);
  int setIDlists();
  int initializationPhase();
  int normalLoadPhase();
  int aggregateLoadPhase();
  int exerciseResidualNorm();
  int exercisePutFunctions();
 
  int save_block_node_soln(DataReader& data, FEI& fei,
         const char* solnFileName, int numProcs,
         int localProc, int solveCounter);
 
  int save_block_elem_soln(DataReader& data, FEI& fei,
         const char* solnFileName,
         int numProcs, int localProc, int solveCounter);
 
  int save_multiplier_soln(DataReader& data, FEI& fei,
         const char* solnFileName,
         int numProcs, int localProc, int solveCounter);
 
  int checkSolution(int localProc, int numProcs,
      const char* solnFileName, const char* checkFileName,
      const char* extension, int solveCounter);
 
  int lsc_matrix_check();
 
  MPI_Comm comm_;
 
  fei::SharedPtr<FEI> fei_;
 
  fei::SharedPtr<LibraryWrapper> wrapper_;
 
  fei::SharedPtr<DataReader> data_;
 
  int localProc_, numProcs_;
 
  int numMatrices;
  int* matrixIDs;
  int numRHSs;
  int* rhsIDs;
  bool useNewFEI_;
};
 
#endif // _FEI_tester_h_