Re: [petsc-dev] QR factorization of dense matrix

2017-10-28 Thread Franck Houssen
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

Re: [petsc-dev] QR factorization of dense matrix

2017-10-28 Thread Jose E. Roman
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

Re: [petsc-dev] SLEPc failure

2017-10-28 Thread Jose E. Roman
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

Re: [petsc-dev] SLEPc failure

2017-10-28 Thread Matthew Knepley
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"