IssmVec.h File Reference

Main Vector class for the Issm toolkit. More...

#include "../../shared/Enum/Enum.h"
#include "../../shared/Exceptions/exceptions.h"
#include "../../shared/MemOps/MemOps.h"
#include "./IssmToolkitUtils.h"
#include <math.h>

Go to the source code of this file.

Data Structures
class	IssmVec< doubletype >

Functions
int	IssmVecTypeFromToolkitOptions (void)

Detailed Description

Main Vector class for the Issm toolkit.

Definition in file IssmVec.h.

Function Documentation

IssmVecTypeFromToolkitOptions()

int IssmVecTypeFromToolkitOptions

(

void

)

Definition at line 56 of file IssmToolkitUtils.cpp.

                                       { /*{{{*/
 
   char* vec_type=NULL;
   int   vec_type_enum;
   int   num_procs=0;
   bool  isparallel=false;
 
   /*first, figure out if we are running in parallel: */
   num_procs=IssmComm::GetSize();
   if(num_procs>1)isparallel=true;
 
   /*retrieve vector type as a string, from the Toolkits Options database, similar to what Petsc does. Actually, 
    *we try and stick with the Petsc vector types: */
   vec_type=ToolkitOptions::GetToolkitOptionValue("vec_type");
 
   if (strcmp(vec_type,"mpi")==0){
      vec_type_enum=MpiEnum;
   }
   else if (strcmp(vec_type,"seq")==0){
      if (isparallel) _error_("Dense vector type not supported for parallel runs with num_procs>1");
      else vec_type_enum=SeqEnum;
   }
   else _error_("vector type not supported yet!");
 
   /*free ressources: */
   xDelete<char>(vec_type);
 
   /*return: */
   return vec_type_enum;
} /*}}}*/