source: issm/trunk-jpl/src/c/objects/ExternalResults/PetscVecExternalResult.cpp@ 12562

/*!\file PetscVecExternalResult.c
 * \brief: implementation of the PetscVecExternalResult object
 */

/*header files: */
/*{{{*/
#ifdef HAVE_CONFIG_H
        #include <config.h>
#else
#error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
#endif

#include <stdio.h>
#include <string.h>
#include "../objects.h"
#include "../../EnumDefinitions/EnumDefinitions.h"
#include "../../shared/shared.h"
#include "../../Container/Container.h"
#include "../../include/include.h"
/*}}}*/

/*PetscVecExternalResult constructors and destructor*/
/*FUNCTION PetscVecExternalResult::PetscVecExternalResult(){{{*/
PetscVecExternalResult::PetscVecExternalResult(){
        return;
}
/*}}}*/
/*FUNCTION PetscVecExternalResult::PetscVecExternalResult(int enum_type,IssmPetscVec value){{{*/
PetscVecExternalResult::PetscVecExternalResult(int in_id, int in_enum_type,Vector* in_value,int in_step, IssmDouble in_time){

        id=in_id;
        enum_type=in_enum_type;

        value=NULL;

        if(in_value){
                value=in_value->Duplicate();
                in_value->Copy(value);
        }
        else value=NULL;

        step=in_step;
        time=in_time;
}
/*}}}*/
/*FUNCTION PetscVecExternalResult::~PetscVecExternalResult(){{{*/
PetscVecExternalResult::~PetscVecExternalResult(){
        VecFree(&value);
}
/*}}}*/

/*Object virtual functions definitions:*/
/*FUNCTION PetscVecExternalResult::Echo {{{*/
void PetscVecExternalResult::Echo(void){

        _printLine_("PetscVecExternalResult:");
        _printLine_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")");

}
/*}}}*/
/*FUNCTION PetscVecExternalResult::DeepEcho{{{*/
void PetscVecExternalResult::DeepEcho(void){

        int i;
        _printLine_("PetscVecExternalResult:");
        _printLine_("   id: " << this->id);
        _printLine_("   enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")");
        _printLine_("   step: " << this->step);
        _printLine_("   time: " << this->time);
        VecView(value,PETSC_VIEWER_STDOUT_WORLD);
}
/*}}}*/
/*FUNCTION PetscVecExternalResult::Id{{{*/
int    PetscVecExternalResult::Id(void){ return -1; }
/*}}}*/
/*FUNCTION PetscVecExternalResult::MyRank{{{*/
int    PetscVecExternalResult::MyRank(void){ 
        extern int my_rank;
        return my_rank; 
}
/*}}}*/
/*FUNCTION PetscVecExternalResult::ObjectEnum{{{*/
int PetscVecExternalResult::ObjectEnum(void){

        return PetscVecExternalResultEnum;

}
/*}}}*/
/*FUNCTION PetscVecExternalResult::copy{{{*/
Object* PetscVecExternalResult::copy() {
        
        return new PetscVecExternalResult(this->id,this->enum_type,this->value,this->step,this->time);

}
/*}}}*/

/*PetscVecExternalResult management: */
/*FUNCTION PetscVecExternalResult::WriteData{{{*/
void   PetscVecExternalResult::WriteData(FILE* fid,bool io_gather){

        int     length;
        int     type;
        int     size;
        char   *name      = NULL;
        IssmPDouble *serialvec = NULL;
        extern int my_rank;
        IssmPDouble passiveDouble;

        /*serialize: */
        VecGetSize(this->value,&size);
        VecToMPISerial(&serialvec,this->value);

        /*now, exit if we are not on cpu 0: */
        if(my_rank)return;

        /*First write enum: */
        EnumToStringx(&name,this->enum_type);
        length=(strlen(name)+1)*sizeof(char);
        fwrite(&length,sizeof(int),1,fid);
        fwrite(name,length,1,fid);
        xDelete<char>(name);

        /*Now write time and step: */
        passiveDouble=reCast<IssmPDouble>(time);
        fwrite(&passiveDouble,sizeof(IssmPDouble),1,fid);
        fwrite(&step,sizeof(int),1,fid);

        /*writing a IssmDouble, type is 1, size is 1: */
        type=1;
        
        fwrite(&type,sizeof(int),1,fid);
        fwrite(&size,sizeof(int),1,fid);
        fwrite(serialvec,size*sizeof(IssmPDouble),1,fid);

        /*Free ressources:*/
        xDelete<IssmPDouble>(serialvec);
}
/*}}}*/
/*FUNCTION PetscVecExternalResult::GetResultName{{{*/
void PetscVecExternalResult::GetResultName(char**pname){
        EnumToStringx(pname,this->enum_type);
}
/*}}}*/
/*FUNCTION PetscVecExternalResult::GetStep{{{*/
int PetscVecExternalResult::GetStep(void){

        return this->step;
}
/*}}}*/