On Tue, Oct 25, 2016 at 9:20 PM, Barry Smith wrote:
>
> Olivier,
>
> Ok, so I've run the code in the debugger, but I don't not think the
> problem is with the null space. The code is correctly removing the null
> space on all the levels of multigrid.
>
> I think the error comes from cha
Olivier,
Ok, so I've run the code in the debugger, but I don't not think the problem
is with the null space. The code is correctly removing the null space on all
the levels of multigrid.
I think the error comes from changes in the behavior of GAMG. GAMG is
relatively rapidly moving
> On Oct 25, 2016, at 8:28 PM, Matthew Knepley wrote:
>
> On Tue, Oct 25, 2016 at 7:57 PM, Olivier Mesnard
> wrote:
> On 25 October 2016 at 19:59, Matthew Knepley wrote:
> On Tue, Oct 25, 2016 at 6:22 PM, Barry Smith wrote:
>
> > On Oct 25, 2016, at 5:39 PM, Olivier Mesnard
> > wrote:
> >
On Tue, Oct 25, 2016 at 7:57 PM, Olivier Mesnard wrote:
> On 25 October 2016 at 19:59, Matthew Knepley wrote:
>
>> On Tue, Oct 25, 2016 at 6:22 PM, Barry Smith wrote:
>>
>>>
>>> > On Oct 25, 2016, at 5:39 PM, Olivier Mesnard <
>>> olivier.mesna...@gmail.com> wrote:
>>> >
>>> > On 25 October 201
Dear professor:
I ｐａｒｔｉｔｉｏｎｅｄ ｍｙ ２Ｄ ｃａｒｔｅｓｉａｎ
ｇｒｉｄ ｗｉｔｈ ４ｒｏｗｓ＊４ｃｏｌｓ ＣＰＵＳ．
１２ １３ １４ １５
８ ９ １０ １１
４ ５ ６ ７
０ １ ２ ３
Ｎｏｗ ｉ ｎｅｅｄ ｔｏ ｓｃａｔｔｅｒ ｔｈｅ ｖａｌｕｅｓ
ｂｅｌｏｎｇｉｎｇ ｔｏ ｃｐｕ５ ｔｏ ｅｖｅｒｙ ｃｐｕ ａｌｏｎｇ
On 25 October 2016 at 19:59, Matthew Knepley wrote:
> On Tue, Oct 25, 2016 at 6:22 PM, Barry Smith wrote:
>
>>
>> > On Oct 25, 2016, at 5:39 PM, Olivier Mesnard <
>> olivier.mesna...@gmail.com> wrote:
>> >
>> > On 25 October 2016 at 17:51, Barry Smith wrote:
>> >
>> > Olivier,
>> >
>> > I
On Tue, Oct 25, 2016 at 6:22 PM, Barry Smith wrote:
>
> > On Oct 25, 2016, at 5:39 PM, Olivier Mesnard
> wrote:
> >
> > On 25 October 2016 at 17:51, Barry Smith wrote:
> >
> > Olivier,
> >
> > In theory you do not need to change anything else. Are you using a
> different matrix object for
> On Oct 25, 2016, at 5:39 PM, Olivier Mesnard
> wrote:
>
> On 25 October 2016 at 17:51, Barry Smith wrote:
>
> Olivier,
>
> In theory you do not need to change anything else. Are you using a
> different matrix object for the velocity_ksp object than the poisson_ksp
> object?
>
> T
On 25 October 2016 at 17:51, Barry Smith wrote:
>
> Olivier,
>
> In theory you do not need to change anything else. Are you using a
> different matrix object for the velocity_ksp object than the poisson_ksp
> object?
>
> The matrix is different for the velocity_ksp and the poisson_ksp.
Olivier,
In theory you do not need to change anything else. Are you using a
different matrix object for the velocity_ksp object than the poisson_ksp object?
The code change in PETSc is very little but we have a report from another
CFD user who also had problems with the change so the
Hi all,
We develop a CFD code using the PETSc library that solves the Navier-Stokes
equations using the fractional-step method from Perot (1993).
At each time-step, we solve two systems: one for the velocity field, the
other, a Poisson system, for the pressure field.
One of our test-cases is a 2D
>> Stefano Zampini
>> replied to the issue on github and was stating that there is some intention
>> to get MFEM working with PETSc but there is no specific timeframe.
>
> There’s currently an open pull request for PETSc solvers inside the private
> MFEM repo.
> However, I don’t know when the c
On Tue, Oct 25, 2016 at 1:43 PM, Jed Brown wrote:
> Stefano Zampini writes:
> > MATHYPRE could be a shell wrapping Hypre calls for the moment.
> > However, HypreParCSR and MATAIJ are mostly equivalent formats. As far as
> I know, the main (only?) difference resides in the fact that the diagonal
Stefano Zampini writes:
> MATHYPRE could be a shell wrapping Hypre calls for the moment.
> However, HypreParCSR and MATAIJ are mostly equivalent formats. As far as I
> know, the main (only?) difference resides in the fact that the diagonal term
> of the diagonal part is ordered first in the CSR.
> On Oct 25, 2016, at 12:31 PM, Matthew Knepley wrote:
>
> On Tue, Oct 25, 2016 at 12:19 PM, Stefano Zampini
> wrote:
> I have a working conversion from HypreParCSR to PETSc MPIAIJ format.
> I could add this code to PETSc, maybe in the contrib folder. Barry, what do
> you think?
>
> No, no o
On Oct 25, 2016, at 9:10 PM, Julian Andrej wrote:
> We have an implementation which models a real physical application but
> don't have the manpower to implement different preconditioner themes
> (like with fieldsplit) or try out different time solving schemes
> (which is way too easy with TS).
We have an implementation which models a real physical application but
don't have the manpower to implement different preconditioner themes
(like with fieldsplit) or try out different time solving schemes
(which is way too easy with TS). Specifically we create a bunch of
operators (Boundary and Fac
On Oct 25, 2016, at 8:50 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
>> On Tue, Oct 25, 2016 at 12:19 PM, Stefano Zampini >> wrote:
>>
>>> I have a working conversion from HypreParCSR to PETSc MPIAIJ format.
>>> I could add this code to PETSc, maybe in the contrib folder. Barry, what
>>>
On Tue, Oct 25, 2016 at 12:50 PM, Jed Brown wrote:
> Matthew Knepley writes:
>
> > On Tue, Oct 25, 2016 at 12:19 PM, Stefano Zampini <
> stefano.zamp...@gmail.com
> >> wrote:
> >
> >> I have a working conversion from HypreParCSR to PETSc MPIAIJ format.
> >> I could add this code to PETSc, maybe
@Stefano Zampini:
I am planning to do two things. One is to directly use Krylov solvers
and preconditioners available to PETSc to try out some stuff and make
some matrix manipulations like symmetric diagonal scaling or getting a
submatrix. If that works I will implement few things (preconditio
Matthew Knepley writes:
> On Tue, Oct 25, 2016 at 12:19 PM, Stefano Zampini > wrote:
>
>> I have a working conversion from HypreParCSR to PETSc MPIAIJ format.
>> I could add this code to PETSc, maybe in the contrib folder. Barry, what
>> do you think?
>>
>
> No, no one looks there. Add it to src/
On Tue, Oct 25, 2016 at 12:19 PM, Stefano Zampini wrote:
> I have a working conversion from HypreParCSR to PETSc MPIAIJ format.
> I could add this code to PETSc, maybe in the contrib folder. Barry, what
> do you think?
>
No, no one looks there. Add it to src/mat/utils and make an interface
funct
I have a working conversion from HypreParCSR to PETSc MPIAIJ format.
I could add this code to PETSc, maybe in the contrib folder. Barry, what do you
think?
> We tried a similar approach to get MFEM objects to PETSc and the real
> problem is that all other convenient functions like creating
> gri
We tried a similar approach to get MFEM objects to PETSc and the real
problem is that all other convenient functions like creating
gridfunctions and projections, you have to convert them every time
which is basically nothing more than a bad workaround. Stefano Zampini
replied to the issue on github
I will check that. I am preallocating but it may be that I am not
allocating big enough. I still have to figure out nuance differences
between these formats to solve the bugs. I appreciate your answer and
hope Satish knows it.
Thank you,
Abdullah Ali Sivas
On 2016-10-25 12:15 PM, Matthew Kne
Sherry,
We set '-mat_superlu_dist_fact SamePattern' as default in
petsc/superlu_dist on 12/6/15 (see attached email below).
However, Anton must set 'SamePattern_SameRowPerm' to avoid crash in his
code. Checking
http://crd-legacy.lbl.gov/~xiaoye/SuperLU/superlu_dist_code_html/pzgssvx___a_bglobal_
On Tue, Oct 25, 2016 at 10:54 AM, Abdullah Ali Sivas <
abdullahasi...@gmail.com> wrote:
> Hello,
>
> I want to use PETSc with mfem and I know that mfem people will figure out
> a way to do it in few months. But for now as a temporary solution I just
> thought of converting hypre PARCSR matrices (t
Hello,
I want to use PETSc with mfem and I know that mfem people will figure
out a way to do it in few months. But for now as a temporary solution I
just thought of converting hypre PARCSR matrices (that is what mfem uses
as linear solver package) into PETSc MPIAIJ matrices and I have a
semi-
Anton,
I guess, when you reuse matrix and its symbolic factor with updated
numerical values, superlu_dist requires this option. I'm cc'ing Sherry to
confirm it.
I'll check petsc/superlu-dist interface to set this flag for this case.
Hong
On Tue, Oct 25, 2016 at 8:20 AM, Anton Popov wrote:
> Ho
>>>
Looking for ML at git.ml, hg.ml or a directory starting
with petsc-pkg-ml
Could not locate an existing copy of ML:
['ml-6.2']
<<<
So configure was looking for something with 'petsc-pkg-ml' - and it
did not find it. [there was 'ml-6
Always look in configure.log to see the exact error.
No - configure does not do checksums - but it expcects the package to
be in a certain format [and this can change between petsc versions].
So if you are using a url from petsc-3.5 - with 3.7 -- it might not
work..
So for any version of petsc -
Hong,
I get all the problems gone and valgrind-clean output if I specify this:
-mat_superlu_dist_fact SamePattern_SameRowPerm
What does SamePattern_SameRowPerm actually mean?
Row permutations are for large diagonal, column permutations are for
sparsity, right?
Will it skip subsequent matrix pe
On Tue, Oct 25, 2016 at 6:39 AM, Morten Nobel-Jørgensen wrote:
> Dear Matt
>
> Did you (or anyone else) find time to look at our issue?
>
> We are really looking forward to your answer :)
>
Yes, I had a little difficulty understanding what was going on, but now I
think I see. I am
attaching my m
On 10/25/2016 01:58 PM, Anton Popov wrote:
On 10/24/2016 10:32 PM, Barry Smith wrote:
Valgrind doesn't report any problems?
Valgrind hangs and never returns (waited hours for a 5 sec run) after
entering factorization for the second time.
Before it happens it prints this (attached)
On 10/24/2016 10:32 PM, Barry Smith wrote:
Valgrind doesn't report any problems?
Valgrind hangs and never returns (waited hours for a 5 sec run) after
entering factorization for the second time.
On Oct 24, 2016, at 12:09 PM, Anton Popov wrote:
On 10/24/2016 05:47 PM, Hong wrote:
Dear Matt
Did you (or anyone else) find time to look at our issue?
We are really looking forward to your answer :)
Kind regards,
Morten
From: Matthew Knepley [knep...@gmail.com]
Sent: Wednesday, October 12, 2016 3:41 PM
To: Morten Nobel-Jørgensen
Cc: petsc-users@
Here is a small complaint about the error message "unable to
download" that is given when using
--download-PACKAGENAME=/PATH/TO/package.tar.gz with the wrong
tarball. For example, for my previous install, I was using
petsc-3.5.3 with:
--download-ml=/path/to/ml-6.2-win.tar.gz
Using the same file
> El 19 oct 2016, a las 9:54, Jose E. Roman escribió:
>
>>
>> El 19 oct 2016, a las 0:26, Bikash Kanungo escribió:
>>
>> Hi Jose,
>>
>> Thanks for the pointers. Here's what I observed on probing it further:
>>
>> • The ||B - B^H|| norm was 1e-18. So I explicitly made it Hermitian by
>
38 matches
Mail list logo
