Hello
Using the TSMonitorSet , I am trying to monitor the time-stepping solver. I
use the TSSetConvergedReason() to set a limit on norm error. I found that
even though the error is below the limit at the current step, the solver
continues to solve that step. Then, in the next step, it will be con
Hello
I am trying to solve system of equations using time stepping. I use
TSSetIFunction() to define the DAEs. I tried to add the xdot[] term to all
of variables. As a solution algorithm , I solved with newton line search
with ASM as a preconditioner. In the case, I want to add the xdot[] term to
Yeah! I think! You are completely right! Thanks for your suggestion. I would
definitely look at them!
Thanks
Amir
From: Matthew Knepley [mailto:knep...@gmail.com]
Sent: Thursday, March 22, 2018 2:25 AM
To: Mohammad Hassan Baghaei
Cc: PETSc
Subject: Re: [petsc-users] Global Numbering
: Mohammad Hassan Baghaei
Cc: PETSc
Subject: Re: [petsc-users] Global Numbering
On Wed, Mar 21, 2018 at 1:28 PM, Mohammad Hassan Baghaei
mailto:mhbagh...@mail.sjtu.edu.cn> > wrote:
I am trying to fill a vector based on global coordinate.
You can see how this might be too vague for me to
, March 21, 2018 10:38 PM
To: Mohammad Hassan Baghaei
Cc: PETSc
Subject: Re: [petsc-users] Global Numbering
On Wed, Mar 21, 2018 at 8:18 AM, Mohammad Hassan Baghaei
mailto:mhbagh...@mail.sjtu.edu.cn> > wrote:
Hello
I am trying to create a global vector of coordinates based
Hello
I am trying to create a global vector of coordinates based on the
parallelized dm. I find it really hard to work with the local numbering, how
can I manage to create a kind of numbering on global level so that I could
easily access to the nodes based on the new numbering. Thanks for your tim
Hi
I am trying to fill the global vector . Using single processor, it works
fine. However, switching to several processors, I get the error regarding
the section offset out of range for several points. The short of error
message is below. Can you let me know what causes the offset to become [0,0(
Baghaei
Cc: PETSc
Subject: Re: [petsc-users] Dealing with DMPlexDistribute()
On Thu, Feb 15, 2018 at 5:43 PM, Mohammad Hassan Baghaei
mailto:mhbagh...@mail.sjtu.edu.cn> > wrote:
Hi Matt
In fact, I created a routine for my specific mesh generation. I firstly, create
the DMPlex object at
the distribution in routine,
with this two lines:
DMPlexDistribute(*dm, 0, NULL, &dmDist);
if (dmDist) {DMDestroy(dm); *dm = dmDist;}
Thanks
Amir
From: Matthew Knepley [mailto:knep...@gmail.com]
Sent: Friday, February 16, 2018 6:34 AM
To: Mohammad Hassan Baghaei
Cc: PETSc
Subject
Hi
I am using DMPlex as interface for mesh generation. On single core, I got
around 3 mesh cells. Whenever, I run on multiple core, say 3, in the
output file for DM, VTK, I got 3 times mesh cell and point numbers. Does it
mean that DMPlexDistribute() does not work properly!
Thanks
Amir
, February 14, 2018 5:37 AM
To: Matthew Knepley
Cc: Mohammad Hassan Baghaei ; PETSc
Subject: Re: [petsc-users] Accessing a field values of Staggered grid
On 13 February 2018 at 21:17, Matthew Knepley mailto:knep...@gmail.com> > wrote:
On Tue, Feb 13, 2018 at 3:21 PM, Mohammad Hassan B
Hi
I am filling the local vector from dm , has a section layout. The thing is I
want to know how I can see the field variable values defined on edges, the
staggered grid. In fact, Whenever I output to VTK in preview, I would be
able to see the main grid. But the values which are defined on edges,
Hi
I wanted to use my solve the system of equation using pseudo-transient
continuation. I found that we can refer to this version for driven cavity
example. However, I could not find it within files
(src/snes/examples/tutorials/ex27.c). Would it be possible for you to share
this file?
Thanks
Am
Hello
I want to know whether is it possible that a specific field in the section have
been defined at some time on edges and other times on the vertices. This change
in the dof , I think, may cause problem, especially in the global vector size
of the dm. At times when the field changes to be de
Hello
I am using the DMPlex interface for the solving PDEs. A part of mesh,
considering, is staggered grid, at the location of middle of each edge.
After generation of main grid points, I find it hard to have the staggered
grid at the prescribed location. I want to know how to deal with the
stagg
Hello
I am using DMPlex to construct a fully circular 2D domain for my PDEs. As I
am discretizing the PDEs using finite difference method, I need to know , to
define the layout of PetscSection, how I can find the boundary condition IS.
I know the location of the boundaries in Sieve chart. However,
Hi
I am going to solve the PDEs on a domain of a circle, which I need to use
the polar coordinates. As dm is supposed for Cartesian coordinate system. Is
it possible to define a new coordinate system in that or I need to use
dmplex for that. Thanks for your great time answering me.
17 matches
Mail list logo