Common.cxx

// Copyright (C) 2001-2012 Vivien Mallet
//
// This file is part of the linear-algebra library Seldon,
// http://seldon.sourceforge.net/.
//
// Seldon is free software; you can redistribute it and/or modify it under the
// terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 2.1 of the License, or (at your option)
// any later version.
//
// Seldon is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
// FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
// more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with Seldon. If not, see http://www.gnu.org/licenses/.
 
 
#ifndef SELDON_FILE_COMMON_CXX
 
#include "Common.hxx"
 
#ifdef SELDON_WITH_SLEPC
#include <slepceps.h>
#endif
 
namespace Seldon
{
  string GetExtension(const string& nom)
  {
    size_t index = nom.find_last_of('.');
    if (index == string::npos)
      return string("");
    
    string extension = nom.substr(index+1, nom.size()-index);  
    return extension;
  }
 
  
  string GetBaseString(const string& nom)
  {
    size_t index = nom.find_last_of('.');
    if (index == string::npos)
      return nom;
    
    string base = nom.substr(0, index);
    return base;
  }
  
    
#ifdef SELDON_WITH_HDF5
 
  template <class T>
  hid_t GetH5Type(T& input)
  {
    double d;
    float f;
    int i;
    long l;
    char c;
    unsigned char uc;
    long long ll;
    unsigned int ui;
    unsigned short us;
    unsigned long ul;
    unsigned long long ull;
 
    if (typeid(input) == typeid(d))
      return H5T_NATIVE_DOUBLE;
    if (typeid(input) == typeid(f))
      return H5T_NATIVE_FLOAT;
    if (typeid(input) == typeid(i))
      return H5T_NATIVE_INT;
    if (typeid(input) == typeid(l))
      return H5T_NATIVE_LONG;
    if (typeid(input) == typeid(c))
      return H5T_NATIVE_CHAR;
    if (typeid(input) == typeid(uc))
      return H5T_NATIVE_UCHAR;
    if (typeid(input) == typeid(ll))
      return H5T_NATIVE_LLONG;
    if (typeid(input) == typeid(ui))
      return H5T_NATIVE_UINT;
    if (typeid(input) == typeid(us))
      return H5T_NATIVE_USHORT;
    if (typeid(input) == typeid(ul))
      return H5T_NATIVE_ULONG;
    if (typeid(input) == typeid(ull))
      return H5T_NATIVE_ULLONG;
    else
      throw Error("hid_t GetH5Type(T& input)",
                  "Type has no corresponding native HDF5 datatype.");
  }
#endif
 
 
  void InitSeldon(int argc, char** argv)
  {
    // double precision for writing
    cout.precision(16);
 
    // initialization of MPI
#ifdef SELDON_WITH_MPI
 
#ifdef MONTJOIE_WITHOUT_THREAD
    MPI_Init(&argc, &argv);
#else
    int required = MPI_THREAD_MULTIPLE;
    int provided = -1;
    int rank, print_level = 1;
    MPI_Init_thread(&argc, &argv, required, &provided);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    if ((rank == 0) && (print_level > 0))
      {
        switch (provided)
          {
          case MPI_THREAD_SINGLE:
            cout << "MPI_Init_thread level = MPI_THREAD_SINGLE" << endl;
            break;
          case MPI_THREAD_FUNNELED:
            cout << "MPI_Init_thread level = MPI_THREAD_FUNNELED" << endl;
            break;
          case MPI_THREAD_SERIALIZED:
            cout << "MPI_Init_thread level = MPI_THREAD_SERIALIZED" << endl;
            break;
          case MPI_THREAD_MULTIPLE:
            cout << "MPI_Init_thread level = MPI_THREAD_MULTIPLE" << endl;
            break;
          default:
            cout << "MPI_Init_thread level = ???" << endl;
          }
      }
#endif
 
#ifdef SELDON_WITH_SLEPC
    string help("Slepc interface");  
    SlepcInitialize(&argc, &argv, (char*)0, help.data());
#endif
    
#endif
 
    // initialization of random generator
    // (in order to have the same init for all processors)
    srand(0);
  }
 
  
  int FinalizeSeldon()
  {
#ifdef SELDON_WITH_SLEPC
    SlepcFinalize();
#endif
    
#ifdef SELDON_WITH_MPI
    MPI_Finalize();
#endif
 
    return 0;
  }
 
}  // namespace Seldon.
 
#define SELDON_FILE_COMMON_CXX
#endif