Dear Barry and Matt,
Thanks very much for the reply.
I am a little bit confused about binary write and read as an input and
output method.
For example, I have code1 and code2, where code1 will take code2's
output as an input.
In code2, the following routine is used to save the solution as a binary file:
#undef __FUNCT__
#define __FUNCT__ "DumpSolutionToMatlab"
PetscErrorCode DumpSolutionToMatlab (DMMG*dmmg, char * fn)
{
DALocalInfo info;
PetscViewer viewer;
PetscFunctionBegin;
ierr = DAGetLocalInfo (DMMGGetDA(dmmg),&info);CHKERRQ(ierr);
sprintf(fileName,
"ff_atwqt2ff_tx%i_ty%i_x%i_y%i_nl%i_gt%i_ot%i_s%i.dat",info.mx,info.my,parameters->mxfield,parameters->myfield,parameters->numberOfLevels,(PetscInt)(parameters->currentTime*10000),parameters->timeAccuracyOrder,(PetscInt)(parameters->smoothFactor*100));
ierr =
PetscViewerBinaryOpen(PETSC_COMM_WORLD,fileName,FILE_MODE_WRITE,&viewer);CHKERRQ(ierr);
ierr =
PetscObjectSetName((PetscObject)dmmg[DMMGGetLevels(dmmg)-1]->x,"ff");
ierr = VecView(dmmg[DMMGGetLevels(dmmg)-1]->x,viewer);CHKERRQ(ierr);
ierr = PetscViewerDestroy (viewer); CHKERRQ (ierr);
PetscFunctionReturn(0);
}
In code1, the following routine is used to read the solution to a Vec
named FIELD:
#undef __FUNCT__
#define __FUNCT__ "FormInitialGuess_physical"
PetscErrorCode FormInitialGuess_physical(DMMG dmmg, Vec X)
{
Vec FIELD;
Field **field;
DA da,da2_4;
PetscViewer viewer;
DALocalInfo info,info2;
PetscFunctionBegin;
ierr =
DMCompositeGetEntries((DMComposite)(dmmg->dm),&da,PETSC_IGNORE);CHKERRQ(ierr);
ierr = DAGetLocalInfo(da,&info);CHKERRQ(ierr);
ierr = DACreate2d(PETSC_COMM_WORLD,DA_NONPERIODIC,DA_STENCIL_BOX,
(info.mx-1-parameters->abandonNumber),
(info.my-1-parameters->abandonNumber), PETSC_DECIDE, PETSC_DECIDE, 4,
2, PETSC_NULL, PETSC_NULL, &da2_4);CHKERRQ(ierr);
ierr = DAGetLocalInfo(da2_4,&info2);CHKERRQ(ierr);
sprintf(fileName,
"ff_twqt2ff_tx%i_ty%i_x%i_y%i_nl%i_gt%i_ot%i_s%i.dat",info2.mx,info2.my,parameters->mxgrid-1,parameters->mygrid-1,parameters->numberOfLevels,(PetscInt)(parameters->timeToGenerateGrid*10000),parameters->timeAccuracyOrder,(PetscInt)(parameters->smoothFactor*100));
PetscViewerBinaryOpen(PETSC_COMM_WORLD,fileName,FILE_MODE_READ,&viewer);
VecLoad(viewer,PETSC_NULL,&FIELD);
VecLoadIntoVector(viewer,FIELD);
ierr = DAVecGetArray(da2_4,FIELD,&field);CHKERRQ(ierr);
ierr = DAVecRestoreArray(da2_4,FIELD,&field);CHKERRQ(ierr);
ierr = VecDestroy(FIELD);CHKERRQ(ierr);
ierr = DADestroy(da2_4);CHKERRQ(ierr);
ierr = PetscViewerDestroy(viewer);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
When I call
VecLoad(viewer,PETSC_NULL,&FIELD);
VecLoadIntoVector(viewer,FIELD);
there is an error in proc0:
[0]PETSC ERROR: PetscBinaryRead() line 251 in src/sys/fileio/sysio.c
Read past end of file
In gdb, it shows that
Program received signal SIGABRT, Aborted.
[Switching to Thread 0xb7c7c6b0 (LWP 598)]
0xb7f4c410 in __kernel_vsyscall ()
(gdb) where
#0 0xb7f4c410 in __kernel_vsyscall ()
#1 0xb7ccc085 in raise () from /lib/tls/i686/cmov/libc.so.6
#2 0xb7ccda01 in abort () from /lib/tls/i686/cmov/libc.so.6
#3 0x08733fe9 in PetscAbortErrorHandler (line=251,
fun=0x886a9c9 "PetscBinaryRead", file=0x886a942 "sysio.c",
dir=0x886a94a "src/sys/fileio/", n=66, p=1,
mess=0xbffdbe64 "Read past end of file", ctx=0x0) at errabort.c:62
#4 0x086a8f5a in PetscError (line=251, func=0x886a9c9 "PetscBinaryRead",
file=0x886a942 "sysio.c", dir=0x886a94a "src/sys/fileio/", n=66, p=1,
mess=0x886a9e6 "Read past end of file") at err.c:482
#5 0x086b73a2 in PetscBinaryRead (fd=9, p=0xbffdc764, n=1, type=PETSC_INT)
at sysio.c:251
#6 0x085d8dee in VecLoadIntoVector_Binary (viewer=0x89f2870, vec=0x89f4a00)
at vecio.c:445
#7 0x085d9e76 in VecLoadIntoVector_Default (viewer=0x89f2870, vec=0x89f4a00)
at vecio.c:514
#8 0x085e9d2c in VecLoadIntoVector (viewer=0x89f2870, vec=0x89f4a00)
at vector.c:1031
#9 0x0804f158 in FormInitialGuess_physical (dmmg=0x8999b30, X=0x898b3c0)
at vecviewload_out.c:391
#10 0x08052d05 in DMMGSolve (dmmg=0x89999e0) at damg.c:307
#11 0x0804d479 in main (argc=Cannot access memory at address 0x256
) at vecviewload_out.c:186
When I call
VecLoad(viewer,PETSC_NULL,&FIELD);
// VecLoadIntoVector(viewer,FIELD);
[0]PETSC ERROR: [1]PETSC ERROR: DAVecGetArray() line 53 in
src/dm/da/src/dagetarray.c Vector local size 32 is not compatible with
DA local sizes 36 100
[2]PETSC ERROR: [3]PETSC ERROR: DAVecGetArray() line 53 in
src/dm/da/src/dagetarray.c Vector local size 28 is not compatible with
DA local sizes 24 80
DAVecGetArray() line 53 in src/dm/da/src/dagetarray.c Vector local
size 32 is not compatible with DA local sizes 36 100
DAVecGetArray() line 53 in src/dm/da/src/dagetarray.c Vector local
size 28 is not compatible with DA local sizes 24 80
in gdb:
Program received signal SIGABRT, Aborted.
[Switching to Thread 0xb7c4f6b0 (LWP 840)]
0xb7f1f410 in __kernel_vsyscall ()
(gdb) where
#0 0xb7f1f410 in __kernel_vsyscall ()
#1 0xb7c9f085 in raise () from /lib/tls/i686/cmov/libc.so.6
#2 0xb7ca0a01 in abort () from /lib/tls/i686/cmov/libc.so.6
#3 0x08733fd1 in PetscAbortErrorHandler (line=53,
fun=0x883be60 "DAVecGetArray", file=0x883be6e "dagetarray.c",
dir=0x883be7b "src/dm/da/src/", n=75, p=1,
mess=0xbffcbd54 "Vector local size 32 is not compatible with DA
local sizes 36 100\n", ctx=0x0) at errabort.c:62
#4 0x086a8f42 in PetscError (line=53, func=0x883be60 "DAVecGetArray",
file=0x883be6e "dagetarray.c", dir=0x883be7b "src/dm/da/src/", n=75, p=1,
mess=0x883be8c "Vector local size %D is not compatible with DA
local sizes %D %D\n") at err.c:482
#5 0x0820d20a in DAVecGetArray (da=0x89f1870, vec=0x89f4770,
array=0xbffcc770)
at dagetarray.c:53
#6 0x0804f162 in FormInitialGuess_physical (dmmg=0x898a560, X=0x899b070)
at vecviewload_out.c:392
#7 0x08052ced in DMMGSolve (dmmg=0x898a400) at damg.c:307
#8 0x0804d479 in main (argc=Cannot access memory at address 0x348
) at vecviewload_out.c:186
If I call
// VecLoad(viewer,PETSC_NULL,&FIELD);
VecLoadIntoVector(viewer,FIELD);
[1]PETSC ERROR: VecLoadIntoVector() line 1016 in
src/vec/vec/interface/vector.c Null Object: Parameter # 2
[0]PETSC ERROR: VecLoadIntoVector() line 1016 in
src/vec/vec/interface/vector.c Null Object: Parameter # 2
[2]PETSC ERROR: VecLoadIntoVector() line 1016 in
src/vec/vec/interface/vector.c Null Object: Parameter # 2
[3]PETSC ERROR: VecLoadIntoVector() line 1016 in
src/vec/vec/interface/vector.c Null Object: Parameter # 2
Program received signal SIGABRT, Aborted.
[Switching to Thread 0xb7c256b0 (LWP 1111)]
0xb7ef5410 in __kernel_vsyscall ()
(gdb) where
#0 0xb7ef5410 in __kernel_vsyscall ()
#1 0xb7c75085 in raise () from /lib/tls/i686/cmov/libc.so.6
#2 0xb7c76a01 in abort () from /lib/tls/i686/cmov/libc.so.6
#3 0x08733fc9 in PetscAbortErrorHandler (line=1016,
fun=0x885dc12 "VecLoadIntoVector", file=0x885d6b0 "vector.c",
dir=0x885d6b9 "src/vec/vec/interface/", n=85, p=1,
mess=0xbfc7f9d4 "Null Object: Parameter # 2", ctx=0x0) at errabort.c:62
#4 0x086a8f3a in PetscError (line=1016, func=0x885dc12 "VecLoadIntoVector",
file=0x885d6b0 "vector.c", dir=0x885d6b9 "src/vec/vec/interface/", n=85,
p=1, mess=0x885d6ed "Null Object: Parameter # %d") at err.c:482
#5 0x085e985a in VecLoadIntoVector (viewer=0x8a08d10, vec=0x0)
at vector.c:1016
#6 0x0804f138 in FormInitialGuess_physical (dmmg=0x89a2880, X=0x89b34f0)
at vecviewload_out.c:391
#7 0x08052ce5 in DMMGSolve (dmmg=0x89a2720) at damg.c:307
#8 0x0804d479 in main (argc=Cannot access memory at address 0x457
) at vecviewload_out.c:186
So I am not sure how this binary write and read working?
Did I miss sth?
Thanks a lot!
Rebecca
