ParseToolkitsOptionsx.h File Reference

header file for parsing petsc options file More...

#include "../../classes/classes.h"

Go to the source code of this file.

Functions
void	ParseToolkitsOptionsx (Parameters *parameters, FILE *petscoptionsfid)

Detailed Description

header file for parsing petsc options file

Definition in file ParseToolkitsOptionsx.h.

Function Documentation

ParseToolkitsOptionsx()

void ParseToolkitsOptionsx	(	Parameters *	parameters,
		FILE *	petscoptionsfid
	)

Definition at line 19 of file ParseToolkitsOptionsx.cpp.

                                                            {
 
   char line[1000];
   char word1[1000];
   char word2[1000];
 
   /*intermediary: */
   int   *analyses     = NULL;
   char **strings      = NULL;
   char  *string       = NULL;
   char **toolkits     = NULL;
   char  *toolkit      = NULL;
   char  *newstring    = NULL;
   char  *catstring    = NULL;
   int    numanalyses;
   int    stringlength,toolkitlength;
 
   /*Get my_rank:*/
   int my_rank=IssmComm::GetRank();
 
   if(my_rank==0){
 
      /*Now, go through lines and figure out how many analyses we have: */
      numanalyses=0;
      while(fgets(line, sizeof line, fid)){
         /*skip comments and empty lines: */
         if ((line[0]=='%') || (line[0]=='\n') || (line[0]==' ') || (line[0]=='\t') || (line[0]=='\r')) continue;
         /*ok, first time, we should get an analysis enum, starting with a +: */
         if (line[0]=='+'){
            numanalyses++;
         }
         else{
            continue;
         }
      }
 
      /*Now, allocate analyses and strings: */
      analyses = xNew<int>(numanalyses);
      strings  = xNew<char*>(numanalyses);
      toolkits = xNew<char*>(numanalyses);
      for(int i=0;i<numanalyses;i++) strings[i]  = NULL;
      for(int i=0;i<numanalyses;i++) toolkits[i] = NULL;
 
      /*Go back to beginning of file:*/
      fseek(fid,0,SEEK_SET);
      numanalyses=0;
      while( fgets(line, sizeof line, fid)){
         /*skip comments and empty lines: */
         if ((line[0]=='%') || (line[0]=='\n') || (line[0]==' ') || (line[0]=='\t') || (line[0]=='\r'))continue;
 
         /*Get rid of end of line: */
         line[strlen(line)-1]='\0';
 
         if(line[0]=='+'){ /*this is the analysis line: */
            analyses[numanalyses]=StringToEnumx(&line[1]);  //skip the '+'
            numanalyses++;
            /*Now get toolkit (should be following line)*/
            if(!fgets(line, sizeof line, fid)) _error_("could not read toolkit for analysis " << &line[1]);
            sscanf(line,"-%s %s",&word1[0],&word2[0]);
            if(strcmp(word1,"toolkit")!=0) _error_("toolkit option does not follow analysis");
            toolkit=xNew<char>((strlen(word2)+1));
            xMemCpy<char>(toolkit,word2,(strlen(word2)+1));
            toolkits[numanalyses-1]=toolkit;
            continue;
         }
         else{
            /*this is an option corresponding to analysis numanalyses-1. Add it to the already existing options*/
            if(strings[numanalyses-1]==NULL){
               string=xNew<char>((strlen(line)+1));
               xMemCpy<char>(string,line,(strlen(line)+1));
 
               strings[numanalyses-1]=string;
            }
            else{
               string=strings[numanalyses-1];
               newstring=xNew<char>((strlen(line)+1));
               xMemCpy<char>(newstring,line,(strlen(line)+1));
 
               /*concatenate:*/
               catstring=xNew<char>(strlen(string)+1+strlen(newstring)+1+1); //fit in a space " "
               xMemCpy<char>(catstring,string,(strlen(string)+1));
 
               strcat(catstring," ");
               strcat(catstring,newstring);
               strings[numanalyses-1]=catstring;
               xDelete<char>(newstring);
               xDelete<char>(string);
            }
         }
      }
   }
 
   /*Ok, broadcast to other cpus: */
   ISSM_MPI_Bcast(&numanalyses,1,ISSM_MPI_INT,0,IssmComm::GetComm());
   if(my_rank!=0){
      analyses  = xNew<int>(numanalyses);
      toolkits  = xNew<char*>(numanalyses);
      strings   = xNew<char*>(numanalyses);
   }
   ISSM_MPI_Bcast(analyses,numanalyses,ISSM_MPI_INT,0,IssmComm::GetComm());
   for(int i=0;i<numanalyses;i++){
      char* toolkit=toolkits[i];
      if(my_rank==0){
         if(toolkit==NULL) _error_("PETSc options for analysis " << EnumToStringx(reCast<int>(analyses[i])) << " have been declared but were not found");
      }
      if(my_rank==0)toolkitlength=(strlen(toolkit)+1)*sizeof(char);
      ISSM_MPI_Bcast(&toolkitlength,1,ISSM_MPI_INT,0,IssmComm::GetComm());
      if(my_rank!=0)toolkit=xNew<char>(toolkitlength);
      ISSM_MPI_Bcast(toolkit,toolkitlength,ISSM_MPI_CHAR,0,IssmComm::GetComm());
      if(my_rank!=0)toolkits[i]=toolkit;
 
      char* string=strings[i];
      if(my_rank==0){
         if(string==NULL) _error_("PETSc options for analysis " << EnumToStringx(reCast<int>(analyses[i])) << " have been declared but were not found");
      }
      if(my_rank==0)stringlength=(strlen(string)+1)*sizeof(char);
      ISSM_MPI_Bcast(&stringlength,1,ISSM_MPI_INT,0,IssmComm::GetComm());
      if(my_rank!=0)string=xNew<char>(stringlength);
      ISSM_MPI_Bcast(string,stringlength,ISSM_MPI_CHAR,0,IssmComm::GetComm());
      if(my_rank!=0)strings[i]=string;
   }
 
   /*Ok, out of strings and analyses and numanalyses, create parameters, and plug them into parameters container: */
   parameters->AddObject(new StringArrayParam(ToolkitsOptionsStringsEnum,strings,numanalyses));
   parameters->AddObject(new StringArrayParam(ToolkitsTypesEnum,toolkits,numanalyses));
   parameters->AddObject(new IntVecParam(ToolkitsOptionsAnalysesEnum,analyses,numanalyses));
 
   /*Clean up and return*/
   for(int i=0;i<numanalyses;i++) xDelete<char>(strings[i]);
   for(int i=0;i<numanalyses;i++) xDelete<char>(toolkits[i]);
   xDelete<char*>(strings);
   xDelete<char*>(toolkits);
   xDelete<int>(analyses);
   return;
}