,MAT_IGNORE_ZERO_ENTRIES,PETSC_TRUE)
For each column call VecGetArray(), zero the "small entries", then call
MatSetValues() for that single column.
Barry
On Feb 2, 2024, at 12:28 PM, TARDIEU Nicolas via petsc-users
wrote:
Dear PETSc users,
I consider a sparse rectangular matrix B of size n x m, n
Dear PETSc users,
I consider a sparse rectangular matrix B of size n x m, n
verlap
On Thu, Jan 18, 2024 at 9:24 AM TARDIEU Nicolas via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Dear PETSc Team,
I am trying to understand how to increase the overlap of a matrix.
I wrote the attached petsc4py script where I build a simple matrix and play
with the inc
Envoyé : jeudi 18 janvier 2024 15:29
À : TARDIEU Nicolas
Cc : petsc-users@mcs.anl.gov
Objet : Re: [petsc-users] Undestanding how to increase the overlap
On Thu, Jan 18, 2024 at 9:24 AM TARDIEU Nicolas via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Dear PETSc Team,
I am
Dear PETSc Team,
I am trying to understand how to increase the overlap of a matrix.
I wrote the attached petsc4py script where I build a simple matrix and play
with the increaseOverlap method. Unfortunately, before and after the call,
nothing changes in the index set. FYI, I have tried to mimic
for a fieldsplit in petsc4py
> On 12 Jul 2023, at 6:04 PM, TARDIEU Nicolas via petsc-users
> wrote:
>
> Dear PETSc team,
>
> In the attached example, I set up a block pc for a saddle-point problem in
> petsc4py. The IS define the unknowns, namely some physical quantity (ph
Dear PETSc team,
In the attached example, I set up a block pc for a saddle-point problem in
petsc4py. The IS define the unknowns, namely some physical quantity (phys) and
a Lagrange multiplier (lags).
I would like to attach a near null space to the physical block, in order to get
the best
(only
adding the mass term in DMTSSetI2Function()).
Matthew Knepley writes:
> On Mon, Dec 19, 2022 at 6:12 AM TARDIEU Nicolas via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
>> Dear PETSc users,
>>
>> I plan to solve nonlinear quasi-static probl
Dear PETSc users,
I plan to solve nonlinear quasi-static problems with PETSc. More precisely,
these are solid mechanics problems with elasto-plasticity.
So they do not involve "physical time", rather "pseudo time", which is
mandatory to describe the stepping of the loading application.
In
row/column indices 8 and 9 being selected by the second process.
What is it that you want to do exactly? Play with ISEmbed(), or get A(u, u)
without using A but B instead?
Thanks,
Pierre
> On 19 Oct 2022, at 12:00 PM, TARDIEU Nicolas via petsc-users
> mailto:petsc-users@mcs.anl.gov>&
by the second process.
What is it that you want to do exactly? Play with ISEmbed(), or get A(u, u)
without using A but B instead?
Thanks,
Pierre
> On 19 Oct 2022, at 12:00 PM, TARDIEU Nicolas via petsc-users
> mailto:petsc-users@mcs.anl.gov>> wrote:
>
> Dear PETSc Team,
>
&
022, at 12:00 PM, TARDIEU Nicolas via petsc-users
> wrote:
>
> Dear PETSc Team,
>
> I am trying to use IS embeding in parallel.
> In order to (try to) understand how it works, I have built a simple example,
> attached to this email.
>
> I consider a 20X20 matrix. The dof
Dear PETSc Team,
I am trying to use IS embeding in parallel.
In order to (try to) understand how it works, I have built a simple example,
attached to this email.
I consider a 20X20 matrix. The dof (u, p, t) in global numbering are the
following :
u: 0..9 p: 10..14 t: 15..19
I have
racked when doing
DMPlexInterpolate. It would be great that you could solve this problem.
Yuan
2021年12月13日(月) 6:36 TARDIEU Nicolas via petsc-users
mailto:petsc-users@mcs.anl.gov>>:
Dear Patrick and Matthew,
Thank you very much for your answers. I am gonna try to set up such a test by
as
ill break, but I am guessing you would not
be interested in those anyway.
Thanks,
Matt
Am Fr., 10. Dez. 2021 um 17:04 Uhr schrieb TARDIEU Nicolas via petsc-users
mailto:petsc-users@mcs.anl.gov>>:
Dear PETSc Team,
Following a previous discussion on the mailing list, I'd like to experim
Dear PETSc Team,
Following a previous discussion on the mailing list, I'd like to experiment
with DMPLEX with a very simple non-manifold mesh as shown in the attached
picture : a cube connected to a square by an edge and to an edge by a point.
I have read some of the papers that Matthew et al.
-users@mcs.anl.gov
Objet : Re: [petsc-users] How to construct DMPlex of cells with different
topological dimension?
On Wed, Nov 3, 2021 at 4:20 AM TARDIEU Nicolas via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote:
Dear all,
As a practitioner of structural mechanics, I modestly wish to s
Dear all,
As a practitioner of structural mechanics, I modestly wish to support Yuan's
argument. Model mixing, and consequently finite element dimension mixing, is
(perhaps unfortunately) an inescapable practice of this engineering field. This
has consequences not only on the nature of the
De : knep...@gmail.com
Envoyé : mardi 28 septembre 2021 14:54
À : TARDIEU Nicolas
Cc : petsc-users@mcs.anl.gov
Objet : Re: [petsc-users] GMRES Breakdown
On Tue, Sep 28, 2021 at 8:39 AM TARDIEU Nicolas via petsc-users
mailto:petsc-users@mcs.anl.gov>> wrote
Dear PETSc Team,
We, code_aster's development team, are using PETSc in our application for years
now, mainly for the KSP et PC features. We run tests cases every week without
problems.
After upgrading from PETSc 3.12.3 to 3.15.3, a test case using GMRES failed
with KSP_DIVERGED_BREAKDOWN. By
le
-mat_block_size 3
-matptap_via scalable
-run_type 1
-pc_gamg_repartition false
-pc_gamg_threshold 0.0
-pc_gamg_threshold_scale .25
-pc_gamg_square_graph 1
-check_pointer_intensity 0
-snes_type ksponly
-ex56_dm_view
-options_left
On Wed, Jun 26, 2019 at 8:21 AM TARDIEU Nicolas via petsc-users
Dear PETSc team,
I have run a simple weak scalability test based on canonical 3D elasticity
problem : a cube, meshed with 8 nodes hexaedra, clamped on one of its face and
submited to a pressure load on the opposite face.
I am using the FGMRES ksp with GAMG as preconditioner. I have set the
22 matches
Mail list logo