Hi,
> i want read an array of 10000 scalars from
/home/Desktop/project/psi.mat to petsc. i have called the following
routines:
ierr =
PetscViewerBinaryOpen(PETSC_COMM_WORLD,"psi.mat",FILE_MODE_READ,&fd);CHKERRQ(ierr);
PetscViewerBinaryGetDescriptor(fd,&REV);
PetscBinaryRead(REV,trial,10000,PETSC_SCALAR);
I recommend
a) replacing "psi.mat" with "/home/Desktop/project/psi.mat"
b) checking the error codes of PetscViewerBinaryGetDescriptor() and
PetscBinaryRead() by using ierr and CHKERRQ(ierr);
To provide the filename at runtime, simply pass the appropriate string
instead of a hard-coded string to PetscViewerBinaryOpen(). This is a
basic C and you can find enough tutorials on the internet on how to use
e.g. command line parameters.
Best regards,
Karli
