Dear all, I need to set a pressure reference point in my parallel FSI code by modifying the matrix, and I realise the following PETSc function may allow me to do this: MatDiagonalSet (…) Has this been implemented in FreeFEM (seems I cannot directly call it in my FreeFEM code)? More generally, wh…

There is no real need for the PETSc interface. You can set the diagonal of your local matrix, and then put that into your parallel Mat. matrix Af = ...; real[int] diagofA(Af.n); diagofA = Af.diag; // get the diagonal of the matrix ...; // some operations on diagofA Af.diag = diagofA ; // set the di…

Thank you very much Pierre! I watched your online talk about parallel interpolation between two meshes, is there a apiece of template code of doing this? Best, Yongxing.

[image] Copy variables between different meshes with PETSc General Discussion Dear Sir or Madam, I would like to use different 2 meshes with PETSc. Additionally, variables want to be copied from each other. However, as shown following codes, the direct copy, such as “v=…

Thank pjr, I can see the interpolation of values from the post, but can we have the interpolation matrix like a serial code?

Search for the keyword transferMat.

Thank you very much! I have found an example of PtAP-2d-PETSc.edp which uses transferMat(). It seems this is only implemented in FreeFEM v4.7, does this version have a docker image? I can only use docker image on the university’s ARC system. Best, Yongxing.

https://hub.docker.com/r/freefem/freefem

It seems be working if the following is not an error: Mat Object: 2 MPI processes type: mpiaij rows=1681, cols=1681 total: nonzeros=18721, allocated nonzeros=18721 total number of mallocs used during MatSetValues calls=0 not using I-node (on process 0) routines Mat Object: 2 MPI processes ty…

Good morning. Do we have v4.7 or above for Windows? when I install FreeFem on Ubuntu (make -j), I came across the following errors: \n\n ****** mshmet ****** \n\n grep: WHERE-LD: No such file or directory make[5]: Entering directory ‘/home/yongxing/Downloads/FreeFem-sources-master/3rdparty/mshmet…

Call PETSc matrix operator MatDiagonalSet (…)

General Discussion

prj June 22, 2021, 3:36pm 4

Topic		Replies	Views
Lagrange multiplier for Dirichlet boundary condition General Discussion	30	1118	June 21, 2023
Saving and loading global/distibuted arrays General Discussion	39	2674	August 5, 2020
Changing Matrix components PETSc Parallel General Discussion	10	561	April 2, 2023
Include two different PDE problems in PETSc General Discussion	17	791	March 9, 2023
Copy variables between different meshes with PETSc General Discussion	14	1819	January 11, 2022

Call PETSc matrix operator MatDiagonalSet (…)

Related topics