Re: [petsc-users] Forming a matrix from vectors

[Jose E. Roman]

I would use MatDenseGetColumn() and VecGetArrayRead() to get the two pointers 
and then copy the values with a loop.

Jose

> El 17 feb 2020, a las 9:35, Eda Oktay <eda.ok...@metu.edu.tr> escribió:
> 
> Hello all,
> 
> I am trying to form a matrix whose columns are eigenvectors I have calculated 
> before U = [v1,v2,v3]. Is there any easy way of forming this matrix? My 
> matrix should be parallel and I have created vectors as below, where nev i s 
> the number of requested eigenvalues. So each V[i] represents an eigenvector 
> and I should form a matrix by using V.
> 
> Vec *V;
>   VecDuplicateVecs(vr,nev,&V); 
>   for (i=0; i<nev;i++){
>   ierr = EPSGetEigenpair(eps,i,&kr,NULL,V[i],NULL);  
>   }
> 
> Thanks!
> 
> Eda