Re: [petsc-users] PETSc initialization error

2020-06-25 Thread Junchao Zhang
Sam, The MPI_Comm_create_keyval() error was fixed in maint/master. From the error message, it seems you need to configure --with-log=1 Otherwise, please send your full error stack trace and configure.log. Thanks. --Junchao Zhang On Thu, Jun 25, 2020 at 2:18 PM Sam Guo wrote: > Hi Juncha

Re: [petsc-users] Problem coupling Petsc into OpenFOAM

2020-06-25 Thread Barry Smith
Mark, Thanks for letting us know. We have added it to the PETSc package list. Should any issues come up we'll be glad to help support it. Barry > On Jun 25, 2020, at 4:53 AM, Mark Olesen wrote: > > Late followup, but people may want to know that a PETSc/OpenFOAM interface is > in

Re: [petsc-users] Inaccurate Matrix Inversion using MatMatSolve

2020-06-25 Thread Barry Smith
Please send the code. This should work fine. Barry > On Jun 25, 2020, at 2:43 PM, Rakesh Halder wrote: > > Hi all, > > I'm using PETSc for matrix calculations as part of a model order reduction > code. An algorithm I'm using requires that I compute the explicit inverse of > a matrix, a

Re: [petsc-users] Inaccurate Matrix Inversion using MatMatSolve

2020-06-25 Thread Zhang, Hong via petsc-users
Rakesh, Can you find the condition number of your matrix A? Hong From: petsc-users on behalf of Rakesh Halder Sent: Thursday, June 25, 2020 2:43 PM To: petsc-users@mcs.anl.gov Subject: [petsc-users] Inaccurate Matrix Inversion using MatMatSolve Hi all, I'm usi

[petsc-users] Inaccurate Matrix Inversion using MatMatSolve

2020-06-25 Thread Rakesh Halder
Hi all, I'm using PETSc for matrix calculations as part of a model order reduction code. An algorithm I'm using requires that I compute the explicit inverse of a matrix, as it needs to be used in matrix-matrix products. The matrix is small (I'll show a 5x5 example), and I first find the LU factor

Re: [petsc-users] PETSc initialization error

2020-06-25 Thread Sam Guo
Hi Junchao, I now encountered the same error with parallel. I am wondering if there is a need for parallel fix as well. [1]PETSC ERROR: #1 PetscInitialize() line 969 in ../../../petsc/src/sys/objects/pinit.c PETSC ERROR: Logging has not been enabled. You might have forgotten to call PetscInitial

Re: [petsc-users] Gather and Broadcast Parallel Vectors in k-means algorithm

2020-06-25 Thread Matthew Knepley
On Thu, Jun 25, 2020 at 11:10 AM Eda Oktay wrote: > Dear Richard, > > From now on, I am actually trying to write a parallel k-means > algorithm by using petsc routines (I don't have to use petsc but I > believe it will be easier) and I used the algorithm you mentioned > before about finding clust

Re: [petsc-users] Gather and Broadcast Parallel Vectors in k-means algorithm

2020-06-25 Thread Eda Oktay
Dear Richard, >From now on, I am actually trying to write a parallel k-means algorithm by using petsc routines (I don't have to use petsc but I believe it will be easier) and I used the algorithm you mentioned before about finding cluster centroids. However, something is bothering me: You said th

Re: [petsc-users] [Ext] Re: Question on SLEPc + computing SVD with a "matrix free" matrix

2020-06-25 Thread Ernesto Prudencio via petsc-users
Thanks, Dave. Ernesto. From: Dave May [mailto:dave.mayhe...@gmail.com] Sent: Thursday, June 25, 2020 2:12 AM To: Ernesto Prudencio Cc: Jose E. Roman ; petsc-users@mcs.anl.gov Subject: Re: [petsc-users] [Ext] Re: Question on SLEPc + computing SVD with a "matrix free" matrix On Thu 25. Jun 202

Re: [petsc-users] Problem coupling Petsc into OpenFOAM

2020-06-25 Thread Mark Olesen
Late followup, but people may want to know that a PETSc/OpenFOAM interface is included in the pending OpenFOAM-v2006 release as a git submodule: https://develop.openfoam.com/modules/external-solver Cheers, /mark On 2019-05-24 12:06, Matthew Knepley via petsc-users wrote: On Thu, May 23, 2

Re: [petsc-users] [Ext] Re: Question on SLEPc + computing SVD with a "matrix free" matrix

2020-06-25 Thread Dave May
On Thu 25. Jun 2020 at 08:23, Ernesto Prudencio via petsc-users < petsc-users@mcs.anl.gov> wrote: > Thank you, Jose. > > However, in the case of a "matrix free" matrix, the APIs on PETSc seem to > allow just the implementation of A.v, not of A' . w > > One could create another "matrix free" matrix