ples/tutorials/swarm_ex3.c
> This demonstrates how you push particles from one sub-domain to another.
>
> * src/dm/examples/tutorials/swarm_ex1.c
> This demonstrates how to define a collection operation to gather particles
> from neighbour cells (cells being defined via DMDA)
>
>
ment support, in other words we still have a lot to do.
>
>
> Barry
>
>
>
>
>
> > On Oct 14, 2016, at 9:54 PM, Sang pham van wrote:
> >
> > Hi Barry,
> >
> > Thank your for your answer. I am writing a parallel code for
> smoothed-particle hydrod
use this information we may have
> suggestions on how to proceed.
>
> Barry
>
>
>
> > On Oct 14, 2016, at 9:23 PM, Sang pham van wrote:
> >
> > Hi Barry,
> >
> > In 2 processes case, the problem is simple, as I know all ghost cells of
> partition
the
partition 0? and can I get a list of ghost cells managing by a neighboring
partition?
Please let me know if my question is still not clear.
Many thanks.
On Sat, Oct 15, 2016 at 8:59 AM, Barry Smith wrote:
>
> > On Oct 14, 2016, at 8:50 PM, Sang pham van wrote:
> >
>
Hi,
I am using DM Vec for a FV code, for some reasons, I want to know partition
of all ghost cells of a specific partition. is there a way do that?
Many thanks.
Best,
(DM da,Mat J)*dm/impls/da/fdda.c:#define
> __FUNCT__ "DMCreateMatrix_DA_3d_MPIAIJ_Fill"
> dm/impls/da/fdda.c:PetscErrorCode DMCreateMatrix_DA_3d_MPIAIJ_Fill(DM
> da,Mat J)
>
>
> On 11 February 2016 at 04:08, Matthew Knepley wrote:
>
>> On Wed, Feb 10, 2016 at
The irregular rows is quite many. The matrix really needs to be
preallocated.
Could you show me how to use DMCreateMatrix_DA_3d_MPIAIJ() directly?
Pham
On Feb 11, 2016 9:52 AM, "Matthew Knepley" wrote:
> On Wed, Feb 10, 2016 at 8:44 PM, Sang pham van
> wrote:
>
>> That
That is because my matrix has some rows which need more entries than usual.
Where can i find source code of DMCreateMatrix()?
Pham.
On Feb 11, 2016 8:35 AM, "Matthew Knepley" wrote:
> On Wed, Feb 10, 2016 at 6:14 PM, Sang pham van
> wrote:
>
>> Hi,
>>
>> I
Hi,
I am trying to create a DM matrix with DMCreateMatrix_DA_3d_MPIAIJ()
instead of using DMCreateMatrix().
Which header file should I include to use that routine? also, what is the
source file containing the DMCreateMatrix() routine?
Many thanks in advance.
Pham
is not.
I am using MatSetValuesStencil() to put values in to the matrix.
I think my problem is fitting well with the Box Stencil (3 dof) of PETSc,
am I right?
Many thanks,
Pham
On Wed, Oct 28, 2015 at 9:05 AM, Barry Smith wrote:
>
> > On Oct 27, 2015, at 8:55 PM, Sang pham v
Pham
On Wed, Oct 28, 2015 at 8:48 AM, Matthew Knepley wrote:
> On Tue, Oct 27, 2015 at 8:18 PM, Sang pham van
> wrote:
>
>> Hi Barry,
>>
>> I made my function for preallocating the DM matrix. I am using immersed
>> boundary method to solve a problem of solid mecha
Many thanks, Dave!
Best,
Pham
On Sun, Oct 25, 2015 at 12:46 AM, Barry Smith wrote:
>
> > On Oct 24, 2015, at 1:43 AM, Sang pham van wrote:
> >
> > Thank you, Dave,
> >
> > Can I just call My_DMCreateMatrix_DA_3d_MPIAIJ(DM da, Mat J, void*
> my_context)
cs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/PetscContainerCreate.html#PetscContainerCreate
>
> You would then set your user context pointer inside the container and then
> use PetscObjectCompose() to attach the container to the DM
>
> Thanks,
> Dave
>
>
> On 24 O
matrix for you.
>
>Barry
>
>
>
>
> > On Sep 3, 2015, at 11:28 AM, Sang pham van wrote:
> >
> > Hi,
> >
> > I am using DMCreateMatrix to create matrix from a existed DM object and
> defined stencil.
> > In my code, boundary nodes need to invo
Hi,
I have the below code section, every time it runs through this section,
used memory increases very much and not decreases then. Could you please
let me know what object did I forgot to release in the for-loop?
Many thanks.
PVS
int mu_x[20]={0,1,0,0,2,0,0,1,0,1,3,0,0,2,2,1,0,1,0,1};
int
Hi,
I am using DMCreateMatrix to create matrix from a existed DM object and
defined stencil.
In my code, boundary nodes need to involve many inner nodes, thus matrix
rows corresponding to boundary nodes are almost dense. How can I tell
petsc that those rows need to be preallocated with more entri
Dear everyone,
I want to purchase a high performance workstation which should have from 32
to 64 nodes, price is not so important. The workstation will be used for
cfd simulation purpose.
Would any one help me to chose one? I would very much appreciate your
suggestions.
Thank you.
Dear everyone,
I want to purchase a high performance workstation which should have from 32
to 64 nodes, price is not so important. The workstation will be used for
cfd simulation purpose.
Would any one help me to chose one? I would very much appreciate your
suggestions.
Thank you.
Sang
On Thu,
Hi Barry,
When using -snes_grid_sequence , how does SNES refines the mesh? Does it
detect where should the mesh be refine locally like an adaptive refinement?
Thank you.
S.
On Tue, Jun 17, 2014 at 4:09 PM, Barry Smith wrote:
>
> On Jun 17, 2014, at 3:03 PM, Sang pham van wrote:
&g
convenient for me when some of my routines need such info
from very first moment. I wish PETSc could allow us to pass in the finest
mesh for grid sequencing.
S,
On Tue, Jun 17, 2014 at 3:46 PM, Barry Smith wrote:
>
> On Jun 17, 2014, at 2:30 PM, Sang pham van wrote:
>
> > Thanks
what options should I use to have pure Newton method?
S.
On Tue, Jun 17, 2014 at 3:20 PM, Barry Smith wrote:
>
> On Jun 17, 2014, at 2:12 PM, Sang pham van wrote:
>
> > Hi,
> >
> > I am using DM structure and SNES to solve a 3D problem. In the problem I
> have 3
Hi,
I am using DM structure and SNES to solve a 3D problem. In the problem I
have 3 variables.
I got SNES converged with a grid. Obtain result are physically right.
However when I refine the grid, SNES does not always converge, the reason
of divergence is line search fail or linear solver failed.
I just modify coordinates
attached with the DM, or I need more implementation?
Many thanks.
Sang
On Wed, Apr 16, 2014 at 6:23 PM, Matthew Knepley wrote:
> On Wed, Apr 16, 2014 at 4:48 PM, Sang pham van wrote:
>
>> Hi Jed,
>>
>> I modified the ex43 code to enforce no-sl
Hi Dave,
I guess you are the one contributed the ex42 in KSP's examples. I want to
modify the example to solve for stokes flow driven by volume force in 3D
duct. Please help me to understand the code by answering the following
questions:
1. Firstly, just for confirmation, the equations you're solv
Hi Matt,
Is it true that the option -snes_linesearch_type basic will make SNES run
with the pure Newton method (without line search nor true region)?
Thank you.
Sang
On Wed, Jun 5, 2013 at 1:04 PM, Matthew Knepley wrote:
> On Wed, Jun 5, 2013 at 12:33 PM, Gaurish Telang wrote:
>
>> I am tryin
rectory", "None",
lu->id.ooc_tmpdir, lu->id.ooc_tmpdir, 256, PETSC_NULL);
?in the function:
PetscErrorCode PetscSetMUMPSOptions(Mat F, Mat A)
in the file: mumps.c
So that the users can enter the temporary directory of their choice.
Kind regards,
Pham Van Ha
---
I think the sentence:
The matrix has row and column ownership ranges.
is very important. I did not get it from the manual (maybe my fault).
Kind regards,
Pham Van Ha
From: "petsc-users-request at mcs.anl.gov"
To: petsc-users at mcs.anl.gov
Se
does a great job.
Nevertheless it would be great if you guy (developers) could consider putting
some comments on creation of rectangular matrix.
Thank you for your very quick response.
Kind regards,
Pham Van Ha
From: Jed Brown
To: Pham Van ; PETSc users list
ith single process and
painfully slow when 2 processes employed.
Thank you in advance.
Kind regards,
Pham Van Ha
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120428/9d782e4d/attachment.htm>
Dear Team,
I used to compile petsc with external package mumps in version 3.1
However, today I tried 3.2 and get many of such error massage:
External package SCALAPACK does not work on Microsoft Windows
It happens also with blacs. How can I build mumps without those pakages?
Kind regards,
PVH
30 matches
Mail list logo