void
Write_Error( int* pid,
int* stepno,
int* nshg,
int* numVars,
double* array1 ) {
int isize, nitems;
int iarray[10];
int nfields;
int numparts;
int nprocs;
(void)*pid; /*silence compiler warning*/
nfields = outpar.nsynciofieldswriterestart;
numparts = workfc.numpe;
nprocs = workfc.numpe;
assert(numparts/nprocs == 1);/* Number of parts per proc ...*/
field_flag++;
if(*pid==0) {
printf("\n");
printf("The %d/%d th field to be written is 'errors'\n",field_flag,nfields);
}
isize = (*nshg)*(*numVars);
nitems = 3;
iarray[ 0 ] = (*nshg);
iarray[ 1 ] = (*numVars);
iarray[ 2 ] = (*stepno);
phio_writeheader(f_descriptor, "errors", (void*)iarray, &nitems,
&isize, "double", phasta_iotype);
phio_writedatablock(f_descriptor, "errors", (void*)array1, &isize,
"double", phasta_iotype );
if (field_flag==nfields){
phio_closefile(f_descriptor);
if (*pid==0) {
printf("Last field %d 'errors' finished! \n",nfields);
printf("\n");
}
}
}
char read_solution(double** solutiono, int* size, int* nshgo, int* ndofo,
int nump, int rank, int timestep, int nSyncFiles, char* casedir)
{
int iarray[10];
const char* iformat = "binary";
int ithree=3;
int igeombc;
char* fn;
int nshg;
int ndof;
double* solution;
phio_fp fp;
if( nSyncFiles == 0 )
posixio_setup(&fp, 'r');
else if( nSyncFiles > 0 )
syncio_setup_read(nSyncFiles, &fp);
char rname[1024];
phio_constructName(fp,"restart",rname);
asprintf(&fn,"%s/%s%d.",casedir,rname,timestep);
phio_openfile(fn, fp);
phio_readheader(fp, "solution", (void*) iarray, &ithree, "integer", iformat);
nshg = iarray[0];
ndof = iarray[1];
if(size != NULL)
*size = nshg*ndof;
solution = (double*) malloc(sizeof(double)*nshg*ndof);
phio_readdatablock(fp, "solution", solution, size, "double", iformat);
phio_closefile(fp);
if(solutiono != NULL)
*solutiono = solution;
if(nshgo != NULL)
*nshgo = nshg;
if(ndofo != NULL)
*ndofo = ndof;
free(fn);
return(0);
}
void
Write_PhAvg2( int* pid,
char* filemode,
char* fieldtag,
int* tagsize,
int* iphase,
int* nphasesincycle,
void* array,
char* arraytype,
int* nshg,
int* numvars,
int* stepno) {
int addtagsize=0; /* phase number is added to the name of the field*/
int tagsize2 = 0;
char *fieldlabel = NULL;
char straddtagsize[10] = "\0";
int isize, nitems;
int iarray[10];
char datatype[10];
int nfields;
int numparts;
int nprocs;
(void)*filemode; /*silence compiler warning*/
(void)*nphasesincycle; /*silence compiler warning*/
if(*iphase<10)
addtagsize=1;
else if(*iphase<100)
addtagsize=2;
else if(*iphase<1000)
addtagsize=3;
tagsize2=*tagsize+addtagsize;
fieldlabel = (char *)malloc((tagsize2+1)*sizeof(char));
strncpy(fieldlabel, fieldtag, *tagsize);
fieldlabel[tagsize2] = '\0';
sprintf(straddtagsize,"%d",*iphase);
if(*iphase<10) {
fieldlabel[tagsize2-1]=straddtagsize[0];
}
else if(*iphase<100) {
fieldlabel[tagsize2-2]=straddtagsize[0];
fieldlabel[tagsize2-1]=straddtagsize[1];
}
else if(*iphase<1000) {
fieldlabel[tagsize2-3]=straddtagsize[0];
fieldlabel[tagsize2-2]=straddtagsize[1];
fieldlabel[tagsize2-1]=straddtagsize[2];
}
if(!strncmp(arraytype,"i",1))
strcpy(datatype,"int");
else /* default is double*/
strcpy(datatype,"double");
nfields = outpar.nsynciofieldswriterestart;
numparts = workfc.numpe;
nprocs = workfc.numpe;
assert(numparts/nprocs == 1);/* Number of parts per proc ...*/
field_flag++;
if(*pid==0) {
printf("\n");
printf("The %d/%d th field to be written is '%s'\n",field_flag,nfields,fieldlabel);
}
/* Write solution field ...*/
isize = (*nshg)*(*numvars);
nitems = 3;
iarray[ 0 ] = (*nshg);
iarray[ 1 ] = (*numvars);
iarray[ 2 ] = (*stepno);
phio_writeheader(f_descriptor, fieldlabel, (void*)iarray, &nitems,
&isize, "double", phasta_iotype);
nitems = (*nshg)*(*numvars);
phio_writedatablock(f_descriptor, fieldlabel, array, &isize,
"double", phasta_iotype );
if (field_flag==nfields){
phio_closefile(f_descriptor);
if (*pid==0) {
printf("\n");
}
}
free(fieldlabel);
}
void
Write_Field( int *pid,
char* filemode,
char* fieldtag,
int* tagsize,
void* array,
char* arraytype,
int* nshg,
int* numvars,
int* stepno) {
char *fieldlabel = NULL;
int isize, nitems;
int iarray[10];
char datatype[10];
int nfields;
int numparts;
int nprocs;
(void)*filemode; /*silence compiler warning*/
if(!strncmp(arraytype,"i",1))
strcpy(datatype,"int");
else /* default is double*/
strcpy(datatype,"double");
nfields = outpar.nsynciofieldswriterestart;
numparts = workfc.numpe;
nprocs = workfc.numpe;
assert(numparts/nprocs == 1);/* Number of parts per proc ...*/
fieldlabel = (char *)malloc((*tagsize+1)*sizeof(char));
strncpy(fieldlabel, fieldtag, *tagsize);
fieldlabel[*tagsize] = '\0';
field_flag++;
if(*pid==0) {
printf("\n");
printf("The %d/%d th field to be written is '%s'\n",field_flag,nfields,fieldlabel);
}
/* Write solution field ...*/
isize = (*nshg)*(*numvars);
nitems = 3;
iarray[ 0 ] = (*nshg);
iarray[ 1 ] = (*numvars);
iarray[ 2 ] = (*stepno);
phio_writeheader(f_descriptor, fieldlabel, (void*)iarray, &nitems,
&isize, datatype, phasta_iotype);
nitems = (*nshg)*(*numvars);
phio_writedatablock(f_descriptor, fieldlabel, array, &isize,
datatype, phasta_iotype );
if (field_flag==nfields){
phio_closefile(f_descriptor);
if (*pid==0) {
printf("Last field %d '%s' finished! \n",nfields, fieldtag);
printf("\n");
}
}
free(fieldlabel);
}
void
Write_Restart( int* pid,
int* stepno,
int* nshg,
int* numVars,
double* array1,
double* array2 ) {
const char* magic_name = "byteorder magic number";
int magic_number = 362436;
int isize, nitems;
int iarray[10];
int nfiles;
int nfields;
int numparts;
int nprocs;
int ione = 1;
double iotime = 0;
char filename[255];
/* First, count the number of fields to write and store the result in*/
countfieldstowriterestart();
/* Retrieve and compute the parameters required for SyncIO*/
nfiles = outpar.nsynciofiles;
nfields = outpar.nsynciofieldswriterestart;
numparts = workfc.numpe;
nprocs = workfc.numpe;
assert(numparts/nprocs == 1);/* Number of parts per proc ...*/
bzero((void*)filename,255);
iotime = TMRC();
if(outpar.output_mode == -1 )
streamio_setup_write(&f_descriptor, streamio_get_r());
else if(outpar.output_mode == 0 )
posixio_setup(&f_descriptor, 'w');
else if(outpar.output_mode > 0 )
syncio_setup_write(nfiles, nfields, numparts/nfiles, &f_descriptor);
else
exit(EXIT_FAILURE);
phio_constructName(f_descriptor,"restart",filename);
phstr_appendInt(filename, *stepno);
phstr_appendStr(filename, ".");
phio_openfile(filename, f_descriptor);
field_flag=0;
/* write the magic number*/
phio_writeheader(f_descriptor, magic_name, (void*)&magic_number, &ione,
&ione, "integer", phasta_iotype);
phio_writedatablock(f_descriptor, magic_name, (void*)&magic_number,
&ione, "integer", phasta_iotype );
field_flag++;
/* Write solution field ...*/
isize = (*nshg)*(*numVars);
nitems = 3;
iarray[ 0 ] = (*nshg);
iarray[ 1 ] = (*numVars);
iarray[ 2 ] = (*stepno);
phio_writeheader(f_descriptor, "solution", (void*)iarray, &nitems,
&isize, "double", phasta_iotype);
nitems = (*nshg)*(*numVars);
phio_writedatablock(f_descriptor, "solution", (void*)(array1),
&isize, "double", phasta_iotype );
field_flag++;
/* Write solution field ...*/
isize = (*nshg)*(*numVars);
nitems = 3;
iarray[ 0 ] = (*nshg);
iarray[ 1 ] = (*numVars);
iarray[ 2 ] = (*stepno);
phio_writeheader(f_descriptor, "time derivative of solution",
(void*)iarray, &nitems, &isize, "double", phasta_iotype);
nitems = (*nshg)*(*numVars);
phio_writedatablock(f_descriptor, "time derivative of solution",
(void*)(array2), &isize, "double", phasta_iotype );
field_flag++;
if (field_flag==nfields){
phio_closefile(f_descriptor);
if (*pid==0) {
printf("\n");
}
}
iotime = TMRC() - iotime;
if (workfc.master == workfc.myrank)
printf("time to write restart (seconds) %f\n",iotime);
}
