On Sun, Nov 6, 2022 at 5:31 PM Alexander Lindsay
wrote:
> We sometimes overallocate our sparsity pattern. Matrix assembly will
> squeeze out allocations that we never added into/set. Is there a convenient
> way to determine the size of the densest row post-assembly? I know that we
> could
We need to see configure.log to see what is going on. Can you send it?
Thanks,
Matt
On Sun, Nov 6, 2022 at 4:29 AM Mohammad Ali Yaqteen
wrote:
> Dear Sir/Madam,
>
>
>
> I am installing PETSc on windows but it keeps giving me unexpected errors.
> I want to use it on MS Visual Studio or
On Fri, Nov 4, 2022 at 7:46 AM Matteo Semplice <
matteo.sempl...@uninsubria.it> wrote:
> On 04/11/2022 02:43, Matthew Knepley wrote:
>
> On Thu, Nov 3, 2022 at 8:36 PM Matthew Knepley wrote:
>
>> On Thu, Oct 27, 2022 at 11:57 AM Semplice Matteo <
>> mat
On Fri, Nov 4, 2022 at 4:32 AM Edoardo alinovi
wrote:
> It is working like a charm now!
>
> Is it mandatory to use VecSetValuesBlocked to assemble the rhs? Does the
> Vec need to be of any other type than VECMPI?
>
SetValuesBlocked() is never required. You can always use the normal
versions,
On Thu, Nov 3, 2022 at 8:36 PM Matthew Knepley wrote:
> On Thu, Oct 27, 2022 at 11:57 AM Semplice Matteo <
> matteo.sempl...@uninsubria.it> wrote:
>
>> Dear Petsc developers,
>> I am trying to use a DMSwarm to locate a cloud of points with respect
>>
On Thu, Oct 27, 2022 at 11:57 AM Semplice Matteo <
matteo.sempl...@uninsubria.it> wrote:
> Dear Petsc developers,
> I am trying to use a DMSwarm to locate a cloud of points with respect
> to a background mesh. In the real application the points will be loaded
> from disk, but I have created a
On Thu, Nov 3, 2022 at 7:52 PM Alexander Lindsay
wrote:
> I have errors on quite a few (but not all) processes of the like
>
> [1]PETSC ERROR: - Error Message
> --
> [1]PETSC ERROR: Nonconforming object sizes
>
On Thu, Nov 3, 2022 at 12:16 PM Edoardo alinovi
wrote:
> Hello Jed/Barry/Petsc friends
>
> I am trying to assemble a block matrix with 3x3 in 2D and 4x4 blocks in 3D
> coming from the fully coupled NS equation.
>
> I am not sure I am understanding the example provided here:
>
On Wed, Nov 2, 2022 at 8:57 AM Gong Yujie wrote:
> Dear development team,
>
> Now I'm doing a project about visualization. In the process of
> visualization, the surface mesh is preferred. I have two questions about
> the DMPlex mesh.
>
>
>1. Can I output the 3D volume mesh in DMPlex as a
>
> > Thank you
> >
> > Ali
> >
> > -Original Message-
> > From: Satish Balay
> > Sent: Monday, October 31, 2022 10:56 PM
> > To: Satish Balay via petsc-users
> > Cc: Matthew Knepley ; Mohammad Ali Yaqteen
> >
> > Subject
Thanks,
Matt
> Your help is highly appreciated.
>
> Thank you
> Ali
>
> -Original Message-
> From: Satish Balay
> Sent: Saturday, October 29, 2022 2:11 PM
> To: Mohammad Ali Yaqteen
> Cc: Matthew Knepley ; petsc-users@mcs.anl.gov
> Subject: Re
On Sun, Oct 30, 2022 at 3:52 PM Barry Smith wrote:
>
>In general you should expect similar but not identical conference
> behavior.
>
> I suggest running with all the monitoring you can.
> -ksp_monitor_true_residual
> -fieldsplit_0_monitor_true_residual
On Fri, Oct 28, 2022 at 9:11 AM Mohammad Ali Yaqteen
wrote:
> Dear Sir,
>
>
>
> During the Installation of PETSc in windows, I installed Cygwin and the
> required libraries as mentioned on your website:
>
> However, when I install PETSc using the configure commands present on the
> petsc
>
> Gang
>
>
> -原始邮件-
> *发件人:*"Xiaoye S. Li"
> *发送时间:*2022-10-25 05:47:15 (星期二)
> *收件人:* "Matthew Knepley"
> *抄送:* "赵刚" , petsc-users@mcs.anl.gov
> *主题:* Re: [petsc-users] An issue of extraction of factorization matrices
> in sp
On Thu, Oct 27, 2022 at 11:57 AM Semplice Matteo <
matteo.sempl...@uninsubria.it> wrote:
> Dear Petsc developers,
> I am trying to use a DMSwarm to locate a cloud of points with respect
> to a background mesh. In the real application the points will be loaded
> from disk, but I have created a
On Thu, Oct 27, 2022 at 1:20 PM Yann Jobic wrote:
> There is this factor also in front of the rectangle, which should 1/4
> for a unit rectangle of dim 2, but it's also 0.5.
>
Hi Yann,
The reference element in Petsc is [-1, -1] to [1, 1]. This matches
FEniCS/Firedrake and also makes
it easier
On Sun, Oct 23, 2022 at 2:58 AM 赵刚 wrote:
> Dear developers,
>
> I have another question. How can I get the L and U matrices and store them
> in a file when I call SuperLU through PETSc? Thanks.
SuperLU stores these matrices in its own format. If you want to do I/O with
them, you would
On Sun, Oct 23, 2022 at 12:48 AM 赵刚 wrote:
> Dear PETSc team,
>
> I want to call the multithreading sparse direct solver SuperLU_MT in
> PETSc, Could I download by "--download-superlu-mt"? Or what is a good way
> to support calling SuperLU_MT interface in PETSc? Thank you.
Is it a separate
On Thu, Oct 20, 2022 at 10:27 PM Barry Smith wrote:
>
> Some of the valgrind information does not appear to make sense
>
> PetscMemcpy() is not calling SNESSolve() so I suspect there must be some
> serious corruption of something to this impossible stack trace
>
I ran Valgrind on it, and it
On Thu, Oct 20, 2022 at 6:48 PM Sepideh Kavousi wrote:
> Hello,
>
> I want to solve my 5 PDEs based on finite difference method using
> periodic BC in x-direction and non-periodic in y-direction but I run into
> error (Segmentation Violation, probably memory access out of range).
>
> For this,
On Thu, Oct 20, 2022 at 9:42 AM Gong Yujie
wrote:
> Dear development team,
>
> I'm trying to write a code to solve a partial differential equation. I
> didn't use PETSc's discretization. I have a question about the
> implementation of the boundary condition (Dirichlet boundary condition).
> When
On Wed, Oct 19, 2022 at 1:04 PM Jackie Chan
wrote:
> Dear All,
>
> I hope you're having a nice day.
> In finite element problems, the stiffness matrix and load vector are
> constructed to calculate the displacement vector using DMCreateMatrix and
> DMCreateGlobalVector, respectively.
On Wed, Oct 19, 2022 at 9:26 AM 김성익 wrote:
> Dear PETSc users,
>
>
> I have a question about structure of programming.
>
> My blueprint of Finite Element programming with PETSc solver is below.
> [image: image.png]
> The blue box is whole loop for my FE program.
> There is a loop B that performs
On Wed, Oct 19, 2022 at 9:13 AM 袁煕 wrote:
> Hello,
>
> I am using Mumps to solve a problem with multiple time steps. The matrix
> structure does not change but its value may or may not change during
> those steps. That means I should reuse the symbolic factorization but
> recall numeric
: ''}
>
Can you show the whole output?
Also, did you remember to 'export' it so that it goes to subshells?
Matt
> Best regards,
> Peng Sun
> --
> *From:* Stefano Zampini
> *Sent:* Thursday, October 13, 2022 1:57 PM
> *To:* Matthew Knepley
regards,
> Peng Sun
>
>
> ----------
> *From:* Matthew Knepley
> *Sent:* Thursday, October 13, 2022 6:34 AM
> *To:* Peng Sun
> *Cc:* petsc-users@mcs.anl.gov
> *Subject:* Re: [petsc-users] Issue with single precision complex numbers
> in petsc4
First send configure.log so we can see the setup.
Thanks,
Matt
On Thu, Oct 13, 2022 at 12:53 AM Peng Sun wrote:
> Dear PETSc community,
>
>
> I have a question regarding the single precision complex numbers of
> petsc4py. I configured PETSc with the “--with-scalar-type=complex
>
On Wed, Oct 12, 2022 at 1:38 PM Alfredo J Duarte Gomez
wrote:
> Good morning PETSC users,
>
> I have a current solver that requires the solution of a Laplace equation
> to be reused for all future time steps.
>
> The configuration is axisymmetric with Dirichlet BCs at the top and bottom
>
On Wed, Oct 12, 2022 at 11:04 AM Takahashi, Tadanaga wrote:
> Hi. I am using the snes nasm for the global solver and snes newtonls for
> the local subdomain solver. I am trying to get the total number of Newton
> iterations for just one subdomain. I've tried:
>
> SNESNASMGetSNES(snes,0,);
>
On Tue, Oct 11, 2022 at 6:04 PM Mike Michell wrote:
> Thank you for the reply and checking.
> Indeed, it seems that local-to-local halo is still not implemented for
> DMPlex. But I believe this is a very required feature for large 3D
> simulations with DMPlex.
>
It is possible that it will make
On Mon, Oct 10, 2022 at 10:41 PM Mike Michell wrote:
> Hi, I was wondering if there is any comment on the example file that I can
> refer to.
>
I see the problem. Local2Local is not implemented for Plex. I thought we
had this automated, but it was only
coded for DMDA. It is a fairly mechanical
On Tue, Oct 11, 2022 at 8:52 AM Barry Smith wrote:
>
>
> https://petsc.org/release/faq/#what-does-the-message-hwloc-linux-ignoring-pci-device-with-non-16bit-domain-mean
Ah. Obscure environment variables, not present in the error message, are
always the best way to deal with this. Stay classy,
hwloc is giving this warning on your machine:
hwloc/linux: Ignoring PCI device with non-16bit domain.
Pass --enable-32bits-pci-domain to configure to support such devices
(warning: it would break the library ABI, don't enable unless really
needed).
The PETSc results are fine, so you can use your
On Mon, Oct 10, 2022 at 5:37 PM Alexander Lindsay
wrote:
> I know that PETSc has native support for ASPIN. Has anyone tried MSPIN? I
> wouldn't be surprised if someone has implemented it in user code. Wondering
> what the barriers would be to creating an option like `-snes_type mspin` ?
>
David
umber of
> iterations are different to reach the specified "-eps_tol" and the relative
> residuals are also slightly different. Is this normal? I am using the
> default Krylov-Schur solver and double precision.
>
> Thanks,
> Feng
> --
> *Fr
will fill the values of the halo cells of *x_g *for
> each partition. Then *x_g* has local and halo cells assigned correctly
> and I can use *x_g* to do my computation. Is this what you mean?
>
>
> Yes
>
> Matt
>
>
> Thanks,
> Feng
>
> ---
On Thu, Oct 6, 2022 at 6:47 AM Gong Yujie wrote:
> Dear development team,
>
> I'm trying to write a code to deal with a multi-field problem. Currently I
> find that the vector ordering for the field is (a1,a2,a3,b1,b2,b3), here
> assume a1,a2,a3 belongs to one field and bs for another field.
On Thu, Oct 6, 2022 at 9:16 AM Zongze Yang wrote:
> Hi, everyone,
>
> I am trying to run some test cases with x window, but the x window never
> showed up with command `make -f ./gmakefile test ...`. It seems a default
> option `-nox` is set. How to disable this option for `make test`?
>
Yes,
On Wed, Oct 5, 2022 at 7:40 AM Zongze Yang wrote:
> Matthew Knepley 于2022年10月5日周三 00:33写道:
>
>> On Tue, Oct 4, 2022 at 3:19 PM Zongze Yang wrote:
>>
>>> Hi everyone,
>>>
>>> I am learning how to use the `DMAdaptLabel` for `DMPlex`, and found th
On Tue, Oct 4, 2022 at 3:19 PM Zongze Yang wrote:
> Hi everyone,
>
> I am learning how to use the `DMAdaptLabel` for `DMPlex`, and found the
> example `src/dm/impls/plex/tests/ex20.c` which label one cell to refine.
>
> 1. This example is just a uniform refinement when using the following
>
On Sat, Oct 1, 2022 at 8:51 PM Mike Michell wrote:
> Thank you for the reply. There is that file in src/dm/interface/ftn-auto/
> for me, instead of the path you mentioned.
>
> After "make allfortranstubs" was done and, PETSc reconfigured and
> reinstalled.
>
> However, I still have the same
On Fri, Sep 30, 2022 at 4:14 PM Mike Michell wrote:
> Hi,
>
> As a follow-up to this email thread,
> https://www.mail-archive.com/petsc-users@mcs.anl.gov/msg44070.html
>
> Are DMLocalToLocalBegin() and DMLocalToLocalEnd() really available for
> DMPlex with Fortran on the latest version of PETSc
On Wed, Sep 28, 2022 at 2:27 PM wangzj997 wrote:
> Dear PETSc development team:
>
> Currently, I am learning and trying to use PETSc's KSP to solve
> large-scale sparse linear systems Ax= b, where A is symmetric positive
> definite
>
How do you know it is positive definite? This is unusual when
On Mon, Sep 26, 2022 at 12:40 PM Duan Junming via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Dear all,
>
> I am using PETSc 3.17.4 on a Linux server, compiling
> with --download-exodus --download-hdf5 --download-openmpi
> --download-triangle --with-fc=0 --with-debugging=0
>
On Mon, Sep 26, 2022 at 8:52 AM fujisan wrote:
> Ok, Thank you.
> I didn't know about MatCreateNormal.
>
> In terms of computer performance, what is best to solve Ax=b with A
> rectangular?
> Is it to keep A rectangular and use KSPLSQR along with PCNONE or
> to convert to normal equations using
Another option are the PCPATCH solvers for multigrid, as shown in this
paper: https://arxiv.org/abs/1912.08516
which I believe solves incompressible elasticity. There is an example in
PETSc for Stokes I believe.
Thanks,
Matt
On Mon, Sep 26, 2022 at 5:20 AM 晓峰 何 wrote:
> Are there other
what you mean?
>
Yes
Matt
> Thanks,
> Feng
>
> --
> *From:* Jose E. Roman
> *Sent:* 21 September 2022 13:07
> *To:* feng wang
> *Cc:* Matthew Knepley ; petsc-users@mcs.anl.gov <
> petsc-users@mcs.anl.gov>
> *Subject:* Re: [petsc-users] Slepc, shell matrix
On Wed, Sep 21, 2022 at 7:41 AM feng wang wrote:
> Hello,
>
> I am using Slepc with a shell matrix. The sequential version seems working
> and now I am trying to make it run in parallel.
>
> The partition of the domain is done, I am not sure how to do the halo
> exchange in the shell matrix in
On Thu, Sep 8, 2022 at 7:44 AM Quentin Chevalier <
quentin.cheval...@polytechnique.edu> wrote:
> Hello PETSc users !
>
> I'm trying to use PETSc routine to read and write complex vectors to
> binary files in parallel.
>
> I'm using dolfinx and petsc4py inside a docker container. My end goal is
>
h the rule that face labels cause any cell in
the support to be labeled, and labeling a cell also labels its faces
2) Partition with a cell overlap
Thanks,
Matt
> [image: image.png]
>
>
> Sincerely
> Nicholas
>
>
>
> On Fri, Sep 2, 2022 at 7:27 AM Matthew Knepley w
sion_debug.cpp:35
> application called MPI_Abort(MPI_COMM_SELF, 77) - process 0
> [unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=77
> :
> system msg for write_line failure : Bad file descriptor
>
>
>
> On Thu, Sep 1, 2022 at 5:32 AM Matthew Knepley wrote:
>
>> On Wed
nt work.
Thanks,
Matt
> Thanks
> Nicholas
>
> On Wed, Aug 31, 2022 at 7:45 AM Matthew Knepley wrote:
>
>> On Tue, Aug 30, 2022 at 6:37 PM Nicholas Arnold-Medabalimi <
>> narno...@umich.edu> wrote:
>>
>>> Hi
>>>
>>> Than
On Wed, Aug 31, 2022 at 9:40 AM wangxq2020--- via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hi
> I run the program on the debuged version of PETSC successfuly. Then I want
> using the optimized version. And configure by ./configure
> --with-debugging=no
On Tue, Aug 30, 2022 at 6:37 PM Nicholas Arnold-Medabalimi <
narno...@umich.edu> wrote:
> Hi
>
> Thanks for all the help. I'm looking to tag several cells based on metrics
> occurring at run time. In addition to the metric-defined cells I also want
> to label the surrounding cells. This is no
On Mon, Aug 29, 2022 at 10:22 AM Lawrence Mitchell wrote:
> On Mon, 29 Aug 2022 at 14:55, Matthew Knepley wrote:
> >
> > On Sun, Aug 28, 2022 at 7:17 PM Matthew Knepley
> wrote:
> >>
> >> On Sun, Aug 28, 2022 at 5:36 PM Mike Michell
> wrote:
> >&g
On Sun, Aug 28, 2022 at 12:44 PM Nicholas Arnold-Medabalimi <
narno...@umich.edu> wrote:
> Hi
>
> I'm working on a finite volume implementation using DMPlex. I was
> wondering what the best-considered way to add cell variables to the DM that
> we don't want to be part of the ksp/snes solve. As
On Sun, Aug 28, 2022 at 5:36 PM Mike Michell wrote:
> Thank you for the reply.
>
> I think it can be more helpful for me if the attached sample code
> (DMInterpolation_Mod.tar) could be checked by you.
>
Okay, you are right. I will run it tomorrow.
Thanks,
Matt
> If you run the
On Sun, Aug 28, 2022 at 12:11 PM Mike Michell wrote:
> Thank you for the reply.
>
> *I cannot quite understand. Are you saying that you have inhomogeneous
> Dirichlet conditions on the boundary, and a 0 guess in the interior, and
> you get all zeros from interpolation? Yes, we have no way of
On Sun, Aug 28, 2022 at 8:51 AM Mike Michell wrote:
> Hi, thank you for the reply.
>
> I was able to manage mapping from cell-center to vertex. Basically, in
> Fortran, it seems DMCreateInterpolation() requires the optional scaling
> vector as a mandatory argument, which is strange.
>
It needs
On Sat, Aug 27, 2022 at 7:17 AM wangxq2020--- via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hi
>
> Submitting a program fails by srun using slurm. This is my command
> srun -n 8 --partition=test --qos=testqos ./test1 -options_file
> options/test1.opts
> and returns
> ./test1: error while
On Thu, Aug 25, 2022 at 7:12 PM Mike Michell wrote:
> Hi, this is a duplication of
> https://lists.mcs.anl.gov/pipermail/petsc-users/2022-August/046746.html
> for in-depth question.
>
> I wrote a short code as attached to test interpolation between two DMPlex
> objects. Goal is to map solution
On Thu, Aug 25, 2022 at 4:25 PM Mike Michell wrote:
> Thank you for the quick note.
> To use PetscFECreateLagrange(), it is unclear about the meaning of "Nc",
> which is explained as 'number of components' on the manual page.
> I looked at an example with Nc=1 in
>
On Thu, Aug 25, 2022 at 1:08 PM Matthew Knepley wrote:
> On Thu, Aug 25, 2022 at 12:38 PM Alexandre Halbach <
> alexandre.halb...@gmail.com> wrote:
>
>> Hello,
>>
>> I am using the metis package provided by the latest petsc version, using
>> --downloa
On Thu, Aug 25, 2022 at 12:38 PM Alexandre Halbach <
alexandre.halb...@gmail.com> wrote:
> Hello,
>
> I am using the metis package provided by the latest petsc version, using
> --download-metis (Linux, 32 bits). This question is arguably more a metis
> issue but since metis is used in petsc and
On Thu, Aug 25, 2022 at 12:52 PM Mike Michell wrote:
> Hi, I am trying to find the function you commented,
> DMPlexCreateInterpolator(), from DMPlex manual page, but I do not think
> this function is visible. Any comment about this?
>
Sorry, it is DMCreateInterpolation(), which in turn calls
otection -fcf-protection -O3
> -fPIC -O2 -flto=auto -ffat-lto-objects -fexceptions -g -grecord-gcc-switches
> -pipe -Wall -Werror=format-security -Wp,-D_FORTIFY_SOURCE=2
> -Wp,-D_GLIBCXX_ASSERTIONS -specs=/usr/lib/rpm/redhat/redhat-hardened-cc1
> -fstack-protector-strong -specs=/usr/
R}/lib/petsc/conf/rules
Thanks,
Matt
> Best regards,
>
> Rafel Amer
>
>
>
>
> El 23/8/22 a les 14:46, Matthew Knepley ha escrit:
>
> Can you run anything in parallel? Say the small sample code that
> calculates pi?
>
>
> https://www.cs.usf
Can you run anything in parallel? Say the small sample code that calculates
pi?
https://www.cs.usfca.edu/~mmalensek/cs220/schedule/code/week09/pi-mpi.c.html
Thanks,
Matt
On Tue, Aug 23, 2022 at 6:34 AM Rafel Amer Ramon wrote:
>
> Hi,
>
> mpicc and mpirun are from the package
On Mon, Aug 22, 2022 at 12:12 PM Patrick Alken
wrote:
> Thank you, I have read that document. I have changed my criteria to:
>
>
> if (j >= first && (j - first) < m) {
>
As Barry says, this should be
if ((j >= first) && (j < last)) {
> /*diagonal*/
>
> else
>
> /*off-diagonal*/
>
>
> It
On Wed, Aug 17, 2022 at 8:31 AM Nicholas Arnold-Medabalimi <
narno...@umich.edu> wrote:
> Good Morning
>
> I have been working with the transitive closure to check to see if a node
> point has a cell that is rooted on a process. However, I think I am having
> a more fundamental issue. Even just
try removing and reconfiguring?
Then I can see everything.
Thanks,
Matt
> Regarding 2.: We need --with-clanguage=cxx since we're also using an
> external C++ library (libMesh). We've been using --with-clanguage=cxx in
> PETSc with libMesh for years...
>
>
>
> On Fri, Aug 12
On Fri, Aug 12, 2022 at 1:14 PM David Knezevic via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> I get a configure error when building PETSc 3.17.4 on Ubuntu 22.04. The
> configure.log is attached.
>
> This error seems to only occur when we use --with-clanguage=cxx. We need
> to use that object
On Fri, Aug 12, 2022 at 9:27 AM David Scott wrote:
> On 12/08/2022 11:26, Matthew Knepley wrote:
>
> This email was sent to you by someone outside the University.
> You should only click on links or attachments if you are certain that the
> email is genuine and the content is safe
On Fri, Aug 12, 2022 at 2:21 AM Quynh Nguyen
wrote:
> Dear petsc developers and users,
>
> I am a new petsc user ^_^.
> I have some questions about using DMDA_BOUNDARY_GHOSTED.
>
> 1. I have 2D Poisson code with domain (nx, ny), and I used
> DMDA_BOUNDARY_NONE, then the size of the matrix is
On Wed, Aug 10, 2022 at 10:00 PM Nicholas Arnold-Medabalimi <
narno...@umich.edu> wrote:
> Hi
>
> I have been working on writing an output of a DMPlex object to interface
> with an existing code that has some strict requirements. I need each
> process to output its partition vertices and cell map
On Tue, Aug 9, 2022 at 11:12 AM Stefano Zampini
wrote:
> PETSc distributes matrices and vectors in parallel. Take a look at
> https://petsc.org/release/docs/manualpages/Vec/VecGetOwnershipRange.html
Longer explanation:
Vectors are distributed by giving a stretch of contiguous rows to each
On Wed, Aug 3, 2022 at 9:17 AM wangxq2020--- via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hello!
> I am a beginner about petsc,and I'm studing the DMDA recently. I hava a
> exercise code below, i run it with 4 processes, why is the size of A matrix
> is 15×15?
>
Are you sure it is not
On Tue, Aug 2, 2022 at 4:00 PM Nicholas Arnold-Medabalimi <
narno...@umich.edu> wrote:
> Hi
>
> Currently, in my distributed DMPlex object, every so often, I want to look
> at total cell counts. To do this in my current implementation, I get the
> cell stratum from the DMplex and search through
Yes, you can use DMPlexCreateSubmesh() (and friends depending on exactly
what kind of submesh you want). This will allow you to create a vector over
only this mesh, and map your volumetric solution to that subvector. Then
you can view the subvector (which pulls in the submesh).
Thanks,
On Mon, Aug 1, 2022 at 11:45 AM E. Ekici wrote:
> Hi,
>
> I would like to preallocate my highly sparse matrices to get speed up
> during assembling.
>
> Up to 7x7, it works fine in parallel.
>
> But when I exceed around 7 rows/columns, I receive following error at
> preallocation
On Thu, Jul 28, 2022 at 6:16 PM Tang, Qi wrote:
> Hi,
>
> Is there some in depth review on different options related to line search
> of snes newtonls? Our solvers are sensitive to the choices there, and I
> could not understand why one works better than the other or how to tune
> those options.
On Mon, Jul 25, 2022 at 10:11 AM Nicholas Arnold-Medabalimi <
narno...@umich.edu> wrote:
> Hi Petsc users,
>
> I have been working on how to read in meshes into a DMPlex object. The
> process of building the cones has been relatively straightforward. The mesh
> files in question have sets of
e in your case) to a natural Vec (before DMPlexDistribute)
>> including constraint dofs.
>> Best,
>> Alexis
>> ---
>> Alexis Marboeuf
>> Postdoctoral fellow, Department of Mathematics & Statistics
>&
On Sun, Jul 24, 2022 at 8:46 AM Duan Junming via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Dear all,
>
>
> I want to add more dofs to the coordinate DM to represent a curved mesh.
>
> I first create a simple box mesh with one cell:
>
>
> PetscCall(DMCreate(comm, ));
>
On Fri, Jul 22, 2022 at 1:32 PM Ahmed Mansur wrote:
> I configured PETSC by downloading the HDF5 library but in my project the
> functions that this library has are not recognized.
> I am trying to use this library to import sparse matrices from MATLAB.
> Thank you very much for the help.
>
in that case.
>
Can you build in scratch and install to user space?
Thanks,
Matt
> On Thu, 21 Jul 2022 at 5:34 PM, Matthew Knepley wrote:
>
>> On Thu, Jul 21, 2022 at 5:16 AM Prateek Gupta
>> wrote:
>>
>>> Hi,
>>> I have a restriction of ins
On Thu, Jul 21, 2022 at 5:16 AM Prateek Gupta
wrote:
> Hi,
> I have a restriction of installing petsc on a community cluster only via
> package manager (APT). But it needs to be configured with lapack and
> parmetis. Is there a way to do this without building from source?
>
Can you just install
On Thu, Jul 21, 2022 at 6:28 AM Emile Soutter
wrote:
> Dear all,
>
> I am struggling with the simple following problem : Having a first matrix
> B1 of size n1xm1, a second matrix B2 of size n2 x m2, build a matrix M of
> size (n1+n2)x(m1+m2) where the blocks B1 and B2 are the "diagonal" of M
>
section includes constraint
>> dofs. I will submit a MR very soon about that. But it will be merged to
>> main.
>>
>> Does this help?
>>
>> Best
>> -----------
>> Alexis Marboeuf
>> Postdoctoral fellow, Department of Mathematics & Statistics
>> Hamilton Hall room 40
number of roots is the number of
local dofs in the global vector. So there would not be 15, which you say is
the global number of unknowns, but rather the local number of unknowns,
which it seems here is 10.
Thanks,
Matt
> Best
>
> On Mon, Jul 18, 2022 at 6:50 PM Matthew Knepley wrot
fining your Section now?
Thanks,
Matt
> Best
>
> On Mon, Jul 18, 2022 at 3:50 PM Matthew Knepley wrote:
>
>> On Mon, Jul 18, 2022 at 12:14 PM Bora Jeong wrote:
>>
>>> Thank you for the corrections. It works.
>>>
>>> However, I s
understand your question. I am not sure what you are using the set
of vertices or cell for.
Thanks,
Matt
> Best
>
> On Sun, Jul 17, 2022 at 9:59 AM Matthew Knepley wrote:
>
>> On Fri, Jul 15, 2022 at 7:05 PM Bora Jeong wrote:
>>
>>> I found that
On Thu, Jul 14, 2022 at 6:33 PM Mike Michell wrote:
> VTK format itself supports pyramids for sure. But I do not think PETSc
> output part through VTK format does not. Is there any possibility to debug
> this issue?
>
VTK output now works with pyramids in
t;
> On Sun, Jul 17, 2022 at 7:59 PM Matthew Knepley wrote:
>
>> On Sun, Jul 17, 2022 at 7:41 AM Prateek Gupta
>> wrote:
>>
>>> Thanks!
>>> On a similar note, is it possible to define data on faces of a parallel
>>> (distributed) dmplex? I have s
unately, Fortran is not very good at reporting declaration errors.
The problem is that you did not include or use the Vec module. I have done
this and your example runs for me. I have included the modified code.
Thanks,
Matt
> Best
>
> On Fri, Jul 15, 2022 at 3:53 PM Matthew Knepley w
erence), where finite volume implementation is the zeroth order special
> case.
>
How is SD different from a spectral element method, with respect to data
layout and communication?
Thanks,
Matt
> Thank you.
> Sincerely,
> Prateek Gupta, PhD
>
>
> On Fri, Jul 15, 202
ion that you get from DMGetGlobalSection(). You can put
any mesh point into the Section and get back an index
using PetscSectionGetOffset().
Maybe give me a more specific example and I can help write some code.
Thanks,
Matt
On Fri, Jul 10, 2015 at 11:14 AM Matthew Knepley wrote:
>
>
ssors? or the number of vertex?
>
Since we need to pass back arrays, you need to pass us in F90 pointers.
Here is an example of doing such a thing:
https://gitlab.com/petsc/petsc/-/blob/main/src/vec/is/sf/tutorials/ex1f.F90#L94
Thanks,
Matt
> Best
>
> On Fri, Jul
On Fri, Jul 15, 2022 at 2:12 PM Randall Mackie
wrote:
>
>
> On Jul 15, 2022, at 11:58 AM, Matthew Knepley wrote:
>
> On Fri, Jul 15, 2022 at 1:46 PM Randall Mackie
> wrote:
>
>> On Jul 15, 2022, at 11:20 AM, Matthew Knepley wrote:
>>
>> On Fri,
On Fri, Jul 15, 2022 at 1:46 PM Randall Mackie
wrote:
> On Jul 15, 2022, at 11:20 AM, Matthew Knepley wrote:
>
> On Fri, Jul 15, 2022 at 11:01 AM Randall Mackie
> wrote:
>
>> I am also interested in converting my DMDA code to DMPlex so that I can
>> use OcTree gri
601 - 700 of 6161 matches
Mail list logo