PICH_NUMVERSION: 3042
so you want to use that mpiexec. This is what will be used if you try
make check
in $PETSC_DIR.
Thanks,
Matt
> Thanks,
> Matt
>
>
>
>
> On 25 Jul 2024, at 11:26 AM, Matthew Knepley wrote:
>
> On Wed, Jul 24, 2024 at 8:37 PM Matthew T
petsc expects which could explain the memory usage.
>
> Do you have any ideas how to fix this?
>
Yes. First we determine what MPI you configured with. Send configure.log,
which has this information.
Thanks,
Matt
> Thanks,
> Matt
>
> On 24 Jul 2024, at 8:41 P
t sure what
is taking up the rest of the memory, but I do not
think it is PETSc from the log you included.
Thanks,
Matt
> I have attached the log view as well if that helps.
>
> Thanks,
> Matt
>
>
>
>
> On 23 Jul 2024, at 9:24 PM, Matthew Knepley wrote:
>
>
Also, you could run with
-mat_view ::ascii_info_detail
and send the output for both cases. The storage of matrix values is not
redundant, so something else is
going on. First, what communicator do you use for the matrix, and what
partitioning?
Thanks,
Matt
On Mon, Jul 22, 2024 at
On Thu, Jul 18, 2024 at 3:07 AM Ivan Luthfi wrote:
> Hi friend, I get many warning (but its ok for the warning). However I
> didnt get the result of my code when i compile it, is there any possible
> mistake in my makefile? can you please help me? The attached file is my
> error and warning
We cannot see your code. The warning says that you give a value to a
variable, but then never use it. We cannot tell if that is true without
looking at the code.
Thanks,
Matt
On Tue, Jul 16, 2024 at 9:54 AM Ivan Luthfi wrote:
> Hello guys, I am still trying to compile my old multigrid
On Sun, Jul 14, 2024 at 4:55 AM Ivan Luthfi wrote:
> Hi there, I have an issue in compiling my code. Here is the warning:
> MsFEM_poisson2D_DMDA. c: 159: 65: error: cannot convert ‘MatStructure’ to
> ‘Mat’ {aka ‘_p_Mat*’} 159 | ierr = KSPSetOperators(ksp_direct,up.
>
++){
> if(IdxPtr[ii] >= 0) pOwnedPtr[jj]++;
>}
>
> PetscCallMPI(MPI_Allreduce([jj],[jj],1,MPIU_INT,MPI_MAX,mpiComm));
>PetscCall(PetscMalloc1(pOwnedPtr[jj],[jj]));
>kk = 0;
>for(PetscInt ii = pStart;ii < pEnd; ii++){
> if(IdxPtr[ii] >= 0){
>
arator. Stop. " in line 24, which is in that first "$(LINK.C) -o
> ". What is wrong from my makefile?
>
You used spaces instead of a tab at the beginning of that line.
Thanks,
Matt
>
> Pada Jum, 12 Jul 2024 pukul 18.57 Matthew Knepley
> menulis:
>
On Fri, Jul 12, 2024 at 7:25 AM Blauth, Sebastian <
sebastian.bla...@itwm.fraunhofer.de> wrote:
> Dear Barry,
>
>
>
> thanks for the clarification. Oh, that’s unfortunate, that not all TAO
> algorithms use the supplied matrix for the norm (and then probably also not
> for computing inner products
On Fri, Jul 12, 2024 at 4:41 AM Ivan Luthfi wrote:
> I get a warning about an incompatible pointer type when compile a code,
> anyone know how to fix this? make bin_MsFEM_poisson2D_DMDA mpicc -Wall -c
> PCMsFEM. c -isystem/home/ivan/petsc/opt-3. 21. 2/include PCMsFEM. c: In
> function
On Fri, Jul 12, 2024 at 3:16 AM Ivan Luthfi wrote:
> I try to compile my code, but i get this error. Anyone can help me? Here
> is my terminal: $make bin_MsFEM_poisson2D_DMDA mpicc -o
> bin_MsFEM_poisson2D_DMDA MsFEM_poisson2D_DMDA. o UserParameter. o
> FormFunction. o MsFEM. o PCMsFEM. o
On Fri, Jul 12, 2024 at 2:26 AM Ivan Luthfi wrote:
> I have warning: FormatFunction. c: In function 'ComputeStiffnessMatrix':
> FormatFunction. c: 128: 46: warning: '%d' directive writing between 1 and
> 11 bytes into a region of size between 0 and 99 [-Wformat-overflow]
> (line128)
On Mon, Jul 8, 2024 at 10:28 PM Ferrand, Jesus A.
wrote:
> Dear PETSc team: Greetings. I keep working on mesh I/O utilities using
> DMPlex. Specifically for the output stage, I need a solid grasp on the
> global numbers and ideally how to set them into the DMPlex during an input
> operation and
> Best,
> Miguel
>
> On Sat, Jul 6, 2024 at 3:33 PM Matthew Knepley wrote:
>
>> On Fri, Jul 5, 2024 at 3:29 AM Miguel Angel Salazar de Troya <
>> miguel.sala...@corintis.com> wrote:
>>
>>> Hello, I have the Navier-Stokes equation coupled wi
On Fri, Jul 5, 2024 at 3:29 AM Miguel Angel Salazar de Troya <
miguel.sala...@corintis.com> wrote:
> Hello, I have the Navier-Stokes equation coupled with a
> convection-diffusion equation for the temperature. It is a two-way coupling
> because the viscosity depends on the temperature. One way to
gree that this is annoying. Currently, we are
> moving toward PyVista, which can read our HDF5 files directly (and also
> work directly with running PETSc), although this is not done yet.
>
> Thanks,
>
> Matt
>
>
>> Best,
>> Miguel
>>
>>
7eWS!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-Pas_5OTyRd$
>
>
>
>
> *From:* Matthew Knepley
> *Sent:* Monday, July 1, 2024 2:27 PM
> *To:* Blauth, Sebastian
> *Cc:* petsc-users@mcs.anl.gov
> *Subject:* Re: [petsc-users] Question re
On Fri, Jun 28, 2024 at 4:05 AM Blauth, Sebastian <
sebastian.bla...@itwm.fraunhofer.de> wrote:
> Hello everyone,
>
>
>
> I have a question regarding the naming convention using PETSc’s
> PCFieldsplit. I have been following
>
directly with running PETSc), although this is not done yet.
Thanks,
Matt
> Best,
> Miguel
>
> On Jun 27, 2024, at 4:59 AM, Matthew Knepley wrote:
>
> Do you get output when you run an example with that option? Is it possible
> that your current working directory is
Do you get output when you run an example with that option? Is it possible
that your current working directory is not what you expect? Maybe try
putting in an absolute path.
Thanks,
Matt
On Wed, Jun 26, 2024 at 5:30 PM MIGUEL MOLINOS PEREZ wrote:
> Sorry, I did not put in cc
On Tue, Jun 25, 2024 at 12:05 PM Julien BRUCHON via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hi, Based on 'cg. c', I'm trying to develop my own Krylov solver (a
> projected conjugate gradient). I want to integrate this into my C++ code,
> where I already have an interface for PETSC which
On Tue, Jun 25, 2024 at 4:50 AM onur.notonur via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hi, I'm trying to implement a Tetgen mesh importer for my
> Petsc/DMPlex-based solver. I am encountering some issues and suspect they
> might be due to my import process. The Tetgen mesh definitions
On Mon, Jun 24, 2024 at 2:50 PM MIGUEL MOLINOS PEREZ wrote:
> Dear all, I am trying to assemble a matrix A with coefficients which I
> need to assemble the RHS (F) and its Jacobian (J) in a TS type of problem.
> Determining each coefficient of A involves the resolution of a small
> non-linear
hat there were still no MKL outputs even I set
> MKL_VERBOSE to be 1. Although, I think it should be many spmv operations
> when doing KSPSolve(). Do you see the possible reasons?
>
>
>
> SPMV are not reported with MKL_VERBOSE (last I checked), only dense BLAS
> is.
>
>
>
>
with that solution as the initial guess.
Thanks,
Matt
> Thanks,
>
> Samar
>
> On Jun 24, 2024, at 12:15 PM, Matthew Knepley wrote:
>
> On Mon, Jun 24, 2024 at 4:24 AM Samar Khatiwala <
> samar.khatiw...@earth.ox.ac.uk> wrote:
>
>> This Message Is From an
ion since I am quite familiar with
> PETSc. Also, the combination seems to be a powerful pair.
>
>
>
> -Ling
>
>
>
> *From: *Matthew Knepley
> *Date: *Monday, June 24, 2024 at 6:12 AM
> *To: *Zou, Ling
> *Cc: *petsc-users@mcs.anl.gov
> *Subject: *Re: [petsc-
On Mon, Jun 24, 2024 at 4:24 AM Samar Khatiwala <
samar.khatiw...@earth.ox.ac.uk> wrote:
> Hi, Sorry to hijack this thread but I just want to add that this is a more
> general problem that I constantly face with PETSc. Not being able to
> checkpoint the complete state of a solver instance and
On Sun, Jun 23, 2024 at 5:04 PM Zou, Ling via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hi all, I am just curious … any effort trying to include PETSc as
> Modelica’s solution option?
>
> (Modelica forum or email list seem to be quite dead so asking here.)
>
I had not heard of it before. I
e close to the one with KSPGuess,
> specifically
>
> KSPGuess Object: 1 MPI process
>
> type: fischer
>
> Model 1, size 200
>
> However, I found at higher frequency, the # of iteration steps are
> significant higher than the one with KSPGuess, I have attahced bot
>
>
>
> On Jun 13, 2024, at 1:27 PM, Yongzhong Li
> wrote:
>
>
>
> This Message Is From an External Sender
>
> This message came from outside your organization.
>
> Hi Matt,
>
> I have rerun the program with the keys you provided. The system output
On Wed, Jun 12, 2024 at 6:36 PM Yongzhong Li
wrote:
> Dear PETSc’s developers, I hope this email finds you well. I am currently
> working on a project using PETSc and have encountered a performance issue
> with the KSPSolve function. Specifically, I have noticed that the time
> taken by KSPSolve
On Sun, Jun 2, 2024 at 10:27 AM Matthew Knepley wrote:
> On Sat, Jun 1, 2024 at 11:39 PM Carpenter, Mark H. (LARC-D302) via
> petsc-users wrote:
>
>> Mark Carpenter, NASA Langley. I am a novice PETSC user of about 10 years.
>> I’ve build a DG-FEM code with petsc as one o
On Sat, Jun 1, 2024 at 11:39 PM Carpenter, Mark H. (LARC-D302) via
petsc-users wrote:
> Mark Carpenter, NASA Langley. I am a novice PETSC user of about 10 years.
> I’ve build a DG-FEM code with petsc as one of the solver paths (I have my
> own as well). Furthermore, I use petsc for MPI
On Tue, May 28, 2024 at 11:04 AM Jose E. Roman via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> It should be: smoother=pc. getMGSmoother(0) The general rule is to drop
> the class name and move "Get" or "Set" to the front (in small letters). But
> sometimes this does not hold. It is better
I put in stuff to propagate the nullspace if you use DM.
Matt
On Thu, May 23, 2024 at 11:04 PM Barry Smith wrote:
> On May 23, 2024, at 3: 48 PM, Stefano Zampini com> wrote: the null space of the Schur complement is the restriction of
> the original null space. I guess if fieldsplit is
em (the poor
> nonlinear solver convergence being caused by missing Jacobian elements
> representing interaction between the sources)?
>
> - Adrian
> On 20/05/24 12:41 pm, Matthew Knepley wrote:
>
> On Sun, May 19, 2024 at 8:25 PM Barry Smith wrote:
>
>> You can call MatS
On Sun, May 19, 2024 at 8:25 PM Barry Smith wrote:
> You can call MatSetOption(mat,MAT_NEW_NONZERO_LOCATION_ERR) then insert
> the new values. If it is just a handful of new insertions the extra time
> should be small. Making a copy of the matrix won't give you a new matrix
> that is any faster
code in the
> GitHub, I would appreciate it if you could have a look at it and guide me
> how to fix it.
> I believe the problem is in the memory of LHS and res_vector because they
> are new vectors I created.
>
> Thank you in advance for your help, I really appreciate it
>
>
On Tue, May 14, 2024 at 9:07 AM Matthew Knepley wrote:
> On Mon, May 13, 2024 at 10:04 PM Adrian Croucher <
> a.crouc...@auckland.ac.nz> wrote:
>
>> On 14/05/24 1:44 pm, Matthew Knepley wrote:
>>
>> I wish GMsh was clearer about what is optional:
>> https
On Mon, May 13, 2024 at 10:04 PM Adrian Croucher
wrote:
> On 14/05/24 1:44 pm, Matthew Knepley wrote:
>
> I wish GMsh was clearer about what is optional:
> https://urldefense.us/v3/__https://gmsh.info/doc/texinfo/gmsh.html*MSH-file-format__;Iw!!
On Tue, May 14, 2024 at 12:14 AM Matteo Semplice <
matteo.sempl...@uninsubria.it> wrote:
> Dear petsc-users, I am playing with DMPlexGetCellCoordinates and observing
> that it returns correct periodic coordinates for cells, but not for faces.
> More precisely, adding
On Mon, May 13, 2024 at 9:33 PM Adrian Croucher
wrote:
> hi, We often create meshes in GMSH format using the meshio library. This
> works OK if we stick to GMSH file format 2. 2. If we use GMSH file format
> 4. 1, DMPlex can't read them because it expects the "Entities" section to
> be present:
On Mon, May 13, 2024 at 1:40 PM Sreeram R Venkat
wrote:
> I have a MatShell object that computes matrix-vector products of a dense
> symmetric matrix of size NxN. The MatShell does not actually form the dense
> matrix, so it is never in memory/storage. For my application, N ranges from
> 1e4 to
On Sun, May 12, 2024 at 10:42 PM Adrian Croucher
wrote:
> hi Matt,
> On 11/05/24 4:12 am, Matthew Knepley wrote:
>
> Thanks. I tried it out and the error below was raised, looks like it's
>> when it tries to destroy the viewer. It still runs ok when
>> DMGetOutputSeque
ide the nonzero pattern of the sparse Jacobian; that is
> what residual components f_i are coupled to what input variables in the
> array x_i. This information comes from your PDE and discretization and
> appears implicitly in your residual function.
>
> Barry
>
>
> I hope
On Fri, May 10, 2024 at 1:01 AM Adrian Croucher
wrote:
> hi Matt,
> On 10/05/24 12:15 am, Matthew Knepley wrote:
>
> I just tried to test it, but there doesn't seem to be a Fortran interface
>> for DMGetOutputSequenceLength().
>>
> Pushed.
>
> Thanks. I
On Wed, May 8, 2024 at 10:46 PM Adrian Croucher
wrote:
> hi Matt,
> On 9/05/24 3:00 am, Matthew Knepley wrote:
>
> Sorry about the delay. I had lost track of this. Can you look at branch
>
> knepley/feature-hdf5-seq-len
>
> I have not made a test yet, but if this wo
eature-hdf5-seq-len
I have not made a test yet, but if this works for you, I will make a test
and merge it in.
Thanks!
Matt
> Regards, Adrian
> On 11/10/21 2:08 pm, Adrian Croucher wrote:
>
> On 10/11/21 11:59 AM, Matthew Knepley wrote:
>
> On Sun, Oct 10, 2021 at
ow why it
would not work for your problem if you don't use it.
Thanks,
Matt
>
> Please find the attached file contains a draft of my code
>
> Thank you in advance for your time and help.
>
>
> Best regards,
>
> Sawsan
>
> ----------
hy not first permute back to the Poisson matrix? Then it
would be trivial.
Thanks,
Matt
> I'm not sure what you mean by "permuted randomly." A random permutation
> of the matrix can kill performance but not math.
>
> Mark
>
>
> On Tue, May 7, 2024 at 8:03 AM Matt
On Tue, May 7, 2024 at 5:12 AM Pierre Jolivet wrote:
> On 7 May 2024, at 9: 10 AM, Marco Seiz wrote: Thanks
> for the quick response! On 07. 05. 24 14: 24, Pierre Jolivet wrote: On 7
> May 2024, at 7: 04 AM, Marco Seiz wrote: This
> Message Is From an External
> ZjQcmQRYFpfptBannerStart
> This
On Mon, May 6, 2024 at 1:14 PM Mark Adams wrote:
> But that will hardwire disabling -options_left. right?
>
True. Your suggestion would only filter out those two in particular.
Thanks,
Matt
> On Mon, May 6, 2024 at 11:30 AM Matthew Knepley wrote:
>
>> On Mon, May
On Mon, May 6, 2024 at 11:15 AM Pierre Jolivet wrote:
> On 6 May 2024, at 3:14 PM, Matthew Knepley wrote:
>
> This Message Is From an External Sender
> This message came from outside your organization.
> On Mon, May 6, 2024 at 1:04 AM Adrian Croucher
> wrote:
>
&g
On Mon, May 6, 2024 at 1:04 AM Adrian Croucher
wrote:
> hi, My code has some optional command line arguments -v and -h for output
> of version number and usage help. These are processed using Fortran's
> get_command_argument(). Since updating PETSc to version 3. 21, I get some
> extra warnings
connectivity list.
>
Right, but that means that the connectivity you read in is not in local
numbering. You will have to renumber it, so making a copy is probably
necessary anyway. This is the same sort of processing I do for the parallel
read.
Thanks,
Matt
> ----
On Thu, May 2, 2024 at 12:53 PM Samar Khatiwala <
samar.khatiw...@earth.ox.ac.uk> wrote:
> Hello, I have a couple of questions about createWithArray in petsc4py: 1)
> What is the correct usage for creating a standard MPI Vec with it?
> Something like this seems to work but is it right?: On each
tivity lists.
>
You know the number of cells if you can read connectivity. Do you mean that
the format does not tell you
the number of vertices?
Thanks,
Matt
> --
> *From:* Ferrand, Jesus A.
> *Sent:* Wednesday, May 1, 2024 8:17 PM
> *To:* Matth
On Wed, May 1, 2024 at 3:34 PM Ferrand, Jesus A.
wrote:
> Dear PETSc team: For a project that I'm working on, I need to manually
> build a DMPlex. From studying the source code of the various APIs in which
> the plex is built from some supported file format, I get that the workflow
> is this:
On Fri, Apr 26, 2024 at 7:23 AM Karthikeyan Chockalingam - STFC UKRI via
petsc-users wrote:
> Hello, When PETSc is installed with GPU support, will it run only on GPUs
> or can it be run on CPUs (without GPUs)? Currently, PETSc crashes when run
> on CPUs. Thank you. Best regards, Karthik. --
On Tue, Apr 23, 2024 at 4:00 PM Yongzhong Li
wrote:
> Thanks Barry! Does this mean that the sparse matrix-vector products, which
> actually constitute the majority of the computations in my GMRES routine in
> PETSc, don’t utilize multithreading? Only basic vector operations such as
> VecAXPY and
it. It would make a fantastic
contribution to PETSc.
Thanks,
Matt
> Thanks,
> Miguel
>
> On 21 Apr 2024, at 14:53, Matthew Knepley wrote:
>
> On Thu, Apr 18, 2024 at 8:23 AM MIGUEL MOLINOS PEREZ
> wrote:
>
>> Dear all, I am working on the implementation of period
On Thu, Apr 18, 2024 at 8:23 AM MIGUEL MOLINOS PEREZ wrote:
> Dear all, I am working on the implementation of periodic bcc using a
> discretisation (PIC-style). I am working with a test case which consists on
> solving the advection of a set of particles inside of a box (DMDA mesh)
> with
On Fri, Apr 19, 2024 at 3:52 PM Ashish Patel wrote:
> Hi Jed, VmRss is on a higher side and seems to match what
> PetscMallocGetMaximumUsage is reporting. HugetlbPages was 0 for me. Mark,
> running without the near nullspace also gives similar results. I have
> attached the malloc_view and gamg
Smith
> >>> Sent: Saturday, July 15, 2023 8:36 AM
> >>> To: Ng, Cho-Kuen
> >>> Cc: petsc-users@mcs.anl.gov
> >>> Subject: Re: [petsc-users] Using PETSc GPU backend
> >>>
> >>> Cho,
> >>>
> >>>
On Wed, Apr 10, 2024 at 7:49 AM Khurana, Parv
wrote:
> Hello PETSc users community, Thank you in advance for your help as always.
> I am using BoomerAMG from Hypre via PETSc as a part of preconditioner in my
> software (Nektar++). I am trying to understand the time profiling
> information that
On Wed, Apr 3, 2024 at 2:28 PM Ferrand, Jesus A.
wrote:
> Dear PETSc team: (Hoping to get a hold of Jed and/or Matt for this one)
> (Also, sorry for the mouthful below) I'm developing routines that will
> read/write CGNS files to DMPlex and vice versa. One of the recurring
> challenges is the
re is a different idea
> behind it.
>
Yes, that is right.
Thanks,
Matt
> Thanks.
> Noam
> On Tuesday, March 26th, 2024 at 11:17 PM, Matthew Knepley <
> knep...@gmail.com> wrote:
>
> On Tue, Mar 26, 2024 at 2:23 PM Noam T. via petsc-users <
>
On Mon, Apr 1, 2024 at 1:57 PM Shatanawi, Sawsan Muhammad via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hello everyone, I hope this email finds you well. Is there a way we can
> check how the matrix looks like after setting it. I have tried debugging it
> with gdb- break points- and print
fT7DAXF6-0fOpeWMEOJcrVL5grSHp-QCF8PdvRl51$
>
> Este correo electrónico y, en su caso, cualquier fichero anexo al
> mismo, contiene información de carácter confidencial exclusivamente
> dirigida a su destinatario o destinatarios. Si no es UD. el
> destinatario del mensaje, le rue
On Sun, Mar 31, 2024 at 4:08 PM MIGUEL MOLINOS PEREZ wrote:
> Dear all, I am writing a function which store datasets (Vectors) coming
> from a DMSwarm structure into a hdf5 file. This step is done nicely
> write_function(){ PetscViewerHDF5Open(…) PetscViewerHDF5PushTimestepping(…)
>
On Tue, Mar 26, 2024 at 2:23 PM Noam T. via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hello, I am trying to understand the FE Tabulation data obtained from e. g
> . PetscFEComputeTabulation. Using a 2D mesh with a single triangle, first
> order, with vertices (0,0), (0,1), (1,0) (see msh
On Sat, Mar 23, 2024 at 8:03 AM Gong Yujie
wrote:
> Dear PETSc group, I'm reading the DS part for the discretization start
> from SNES ex17. c which is a demo for solving linear elasticity problem. I
> have two questions for the details. The first question is for the residual
> function. Is the
On Fri, Mar 15, 2024 at 9:55 AM Frank Bramkamp wrote:
> Dear PETSc Team, I am using the latest petsc version 3. 20. 5. I would
> like to create a matrix using MatCreateSeqAIJ To insert values, I use
> MatSetValues. It seems that the Fortran interface/stubs are missing for
> MatsetValues, as the
On Thu, Feb 29, 2024 at 4:22 PM Blondel, Sophie via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> I still get the same error when deactivating GPU-aware MPI. I also tried
> unloading spectrum MPI and using openMPI instead (recompiling everything)
> and I get a segfault in PETSc in that case
On Thu, Feb 29, 2024 at 11:03 AM Blondel, Sophie via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Thank you Junchao, Yes, I am using gpu-aware MPI. Is "-use_gpu_aware_mpi
> 0" a runtime option or a compile option? Best, Sophie From: Junchao Zhang
> Sent: Thursday, February 29, 2024 10: 50 To:
Right now, there is no particle identity.
Thanks,
Matt
> Best,
> Miguel
>
> On Feb 25, 2024, at 6:53 PM, Matthew Knepley wrote:
>
> On Sun, Feb 25, 2024 at 7:21 AM MIGUEL MOLINOS PEREZ
> wrote:
>
>> Dear Matthew,
>>
>> Thank you for answer
the residual equations of my problem (please,
> don’t judge me, publish or perish...)
>
> Many thanks,
> Miguel
>
>
> On Feb 24, 2024, at 4:24 PM, Matthew Knepley wrote:
>
> On Fri, Feb 23, 2024 at 5:14 PM MIGUEL MOLINOS PEREZ
> wrote:
>
>> Dear all, I am s
On Fri, Feb 23, 2024 at 5:14 PM MIGUEL MOLINOS PEREZ wrote:
> Dear all, I am struggling on how to include ghost points in DMSwarm local
> vectors. According to PETSc documentation it seems “straightforward” for a
> DMDA mesh. However, I am not so sure on how to do it for a DMSwarm. In
> fact, if
On Fri, Feb 23, 2024 at 9:54 AM 袁煕 wrote:
> Hello, I cannot fetch right values of a local vector after a second time
> calling DMSetLocalSection. My program runs like call
> PetscSectionSetup(section1,ierr) call
> DMSetLocalSection(dm_mesh,section1,ierr)call
>
On Thu, Dec 8, 2022 at 2:56 PM Barry Smith wrote:
>
>You would use *PETSC_NULL_DM*LABEL but Matt needs to customize the
> PETSc Fortran stub for DMAddField() for you to handle accepting the NULL
> from PETSc.
>
Sorry, I missed this message. The customization is done.
Thanks,
Matt
You are right. I was not clearing the caches when the section was reset.
Stefano ahs fixed this now:
https://urldefense.us/v3/__https://gitlab.com/petsc/petsc/-/merge_requests/7314__;!!G_uCfscf7eWS!ZPkq_k9vdzTZxeW53aWVTIcDgeGnywSm0tGKDJkJqdwgTomfYjw1lsTMyRtAjtruSslcgLIoKGg7IILlBxbJ$
On Mon, Feb 12, 2024 at 7:10 AM Jeremy Theler (External) <
jeremy.theler-...@ansys.com> wrote:
> Hi Barry
>
> > The bug fix for 2 is availabel in
> https://gitlab.com/petsc/petsc/-/merge_requests/7279
>
> Note that our code goes through MatZeroRowsColumns_MPIAIJ() not through
>
On Fri, Feb 9, 2024 at 7:04 AM Jeremy Theler (External) <
jeremy.theler-...@ansys.com> wrote:
> Hi all
>
> Because of a combination of settings, our code passes through this line:
>
>
> https://gitlab.com/petsc/petsc/-/blob/main/src/ksp/pc/impls/fieldsplit/fieldsplit.c?ref_type=heads#L692
>
>
t other bug with coordinates (finger
crossed). We are running the tests now.
Thanks,
Matt
> Thanks, best regards,
> Berend.
>
> p.s. If someone is interested in the code to transform a DMForest to a
> DMPlex of type DM_POLY, I am more than happy to share.
>
> On
On Wed, Feb 7, 2024 at 8:30 AM Mark Adams wrote:
> Wait for Matt to chime in, but I think you will need to start from scratch
> and use
> https://petsc.org/release/manualpages/DMPlex/DMPlexBuildFromCellList
>
Let me restate things to make sure I get it. You want to select a submesh
that retains
On Mon, Feb 5, 2024 at 10:10 AM GUBERNATIS Pierre
wrote:
> Hello,
>
> I would appreciate an advice on this task, if possible.
>
>
>
> So far, we use PETSc to ‘couple’ two computations of the same simulation
> code, in a single MPI run.
>
> I simplify but roughly, we have 2 Cartesian domains on
On Mon, Feb 5, 2024 at 8:48 AM TAY Wee Beng wrote:
> Hi all,
>
> Not in my code, nor in PETSc.
>
I would try this. Of course, it would be best to see where exactly the
overflow is happening.
THanks,
Matt
> On 5/2/2024 8:21 pm, Matthew Knepley wrote:
>
> On Mon, F
On Mon, Feb 5, 2024 at 4:29 AM TAY Wee Beng wrote:
> Hi,
>
> I have a code which uses PETSc and it was working fine. However, I
> haven't used the code for a while and now we've migrated to a new cray
> cluster so the compilers have been updated. I'm using gnu fortran.
>
> When running the code
Can you give us a simple code that is not working for you? My test should
work
src/dm/impls/plex/tests/ex26f90
Thanks,
Matt
On Sat, Feb 3, 2024 at 11:36 AM Barry Smith wrote:
>
> The Fortran "stubs" (subroutines) should be in
>
On Wed, Jan 31, 2024 at 10:10 AM Alain O' Miniussi
wrote:
> Hi,
>
> It is indicated in:
> https://petsc.org/release/manualpages/Sys/PetscInitialize/
> that the init function will call MPI_Init.
>
> What if MPI_Init was already called (as it is the case in my application)
>From the page: "
On Wed, Jan 31, 2024 at 8:21 AM Mark Adams wrote:
> Iterative solvers have to be designed for your particular operator.
> You want to look in your field to see how people solve these problems.
> (eg, zeros on the diagonal will need something like a block solver or maybe
> ILU with a particular
left preconditioning
> using NONE norm type for convergence test
> PC Object: (redistribute_sub_) 1 MPI process
> type: lu
> out-of-place factorization
> tolerance for zero pivot 2.22045e-14
>
> On Sat, Jan 27, 2024 at 10:24 AM Matthew Knepley
View the solver.
Matt
On Sat, Jan 27, 2024 at 9:43 AM Mark Adams wrote:
> I am not getting ksp_rtol 1e-12 into pc_redistribute correctly?
>
>
>
> * Linear redistribute_ solve converged due to CONVERGED_RTOL iterations 1
> 0 KSP Residual norm 2.182384017537e+02 1 KSP Residual norm
>
ith v.
Thanks,
Matt
> Thanks,
> Pantelis
> ------
> *Από:* Matthew Knepley
> *Στάλθηκε:* Τρίτη, 23 Ιανουαρίου 2024 5:21 μμ
> *Προς:* Pantelis Moschopoulos
> *Κοιν.:* Barry Smith ; petsc-users@mcs.anl.gov <
> petsc-users@mcs.anl.gov>
getting
MatLRC to do it.
Thanks,
Matt
> Thanks,
> Pantelis
> --
> *Από:* Matthew Knepley
> *Στάλθηκε:* Τρίτη, 23 Ιανουαρίου 2024 3:20 μμ
> *Προς:* Pantelis Moschopoulos
> *Κοιν.:* Barry Smith ; petsc-users@mcs.anl.gov <
> petsc-users@mcs.anl.gov>
> *
upply the IS using the command PCFieldSplitSetIS, I do not specify
> anything in the matrix set up right?
>
You should just need to specify the rows for each field as an IS.
Thanks,
Matt
> Thanks,
> Pantelis
> ------
> *Από:* Matthew Knepley
> *Σ
te the Field split:
> ufields(3) = [0, 1, 2]
> pfields(1) = [3]
>
> call PCSetType(pc, PCFIELDSPLIT, ierr)
> call PCFieldSplitSetBlockSize(pc, 4,ierr)
> call PCFieldSplitSetFields(pc, "0", 3, ufields, ufields,ierr)
> call PCFieldSplitSetFields(pc, "1"
makes sense. You could call DMLocatePoints(), but maybe
you are doing that.
Thanks,
Matt
> Best, Berend.
>
> On 1/22/24 20:03, Matthew Knepley wrote:
> > On Mon, Jan 22, 2024 at 1:57 PM Berend van Wachem <
> berend.vanwac...@ovgu.de <mailto:berend.vanwac...@ovgu.de
me above, you would have to write both DMs. Are you trying to
avoid this?
Thanks,
Matt
> Thanks,
>
> Berend.
>
> On 1/22/24 18:58, Matthew Knepley wrote:
> > On Mon, Jan 22, 2024 at 10:49 AM Berend van Wachem <
> berend.vanwac...@ovgu.de <mailto:berend.vanwac...
1 - 100 of 6176 matches
Mail list logo