I've seen that !... But can't get BVOrthogonalize to work.
I tried:
Mat Z; MatCreateSeqDense(PETSC_COMM_SELF, n, m, NULL, &Z);
...; // MatSetValues(Z, ...)
BVCreate(PETSC_COMM_SELF, &bv);
BVCreateFromMat(Z, &bv); // Z is tall-skinny
Mat R; MatCreateSeqDense(PETSC_COMM_SELF, n, m, NULL, &R); // Sam
Matrix R must be mxm.
BVOrthogonalize computes Z=Q*R, where Q overwrites Z.
Jose
> El 28 oct 2017, a las 13:11, Franck Houssen
> escribió:
>
> I've seen that !... But can't get BVOrthogonalize to work.
>
> I tried:
> Mat Z; MatCreateSeqDense(PETSC_COMM_SELF, n, m, NULL, &Z);
> ...; // MatSetVa
The two matrices are the same ...
> El 28 oct 2017, a las 13:11, Franck Houssen
> escribió:
>
> I just added that before EPSSetOperators:
> PetscViewer viewerA;
> PetscViewerBinaryOpen(PETSC_COMM_WORLD,"Atau.out",FILE_MODE_WRITE,&viewerA);
> MatView(A,viewerA);
> PetscViewer viewerB;
> Petsc
On Sat, Oct 28, 2017 at 7:11 AM, Franck Houssen
wrote:
> I just added that before EPSSetOperators:
> PetscViewer viewerA; PetscViewerBinaryOpen(PETSC_
> COMM_WORLD,"Atau.out",FILE_MODE_WRITE,&viewerA); MatView(A,viewerA);
> PetscViewer viewerB; PetscViewerBinaryOpen(PETSC_
> COMM_WORLD,"Btau.out"