Jie,
Thanks. How does one add a "regular" preconditioner to the deflation
approach?
In the theory of domain decomposition methods one constructs preconditioners by
composing (either additively or multiplicatively) solvers on subspaces of the
entire solution space. For simplicity consid
I guess viewing deflation as a preconditioner may destroy the symmetry of the
system if it is originally symmetric? I wrote a one page technical details
about deflation and it may help clarify my thoughts.
http://www.mcs.anl.gov/~jiechen/tmp/deflate.pdf
Jie
- Original Message -
From:
routine,
> that I am setting
>
> *left->ops->getarray = PETSC_NULL;
>
> Any thoughts on what might be going on?
>
> Thanks,
> Amneet
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/e8737b88/attachment.html>
With the introduction of PetscPowRealInt() should this be moved to
petscmath.h and renamed ?
/* C++ does not promote int64_t to scalar or int32_t for std::pow() */
static PetscScalar Pow(PetscScalar b,PetscInt p)
{
return PetscPowScalar(b,(int)p);
}
st experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which
> their experiments lead.
> -- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/552ad035/attachment.html>
not happening. It is in this routine,
that I am setting
*left->ops->getarray = PETSC_NULL;
Any thoughts on what might be going on?
Thanks,
Amneet
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments
Hi,
Having read that petsc-dev has the possibility of using OpenMP to run PC
and KSP routines in parallel, and that people are getting nice speedup
with this, I wanted to make use of this feature in our code. It is an
incompressible two-phase flow solver programmed in Fortran which used
PETSc
On Feb 28, 2013, at 7:30 AM, John Mousel wrote:
> Let me rephrase. Is a functionality like this something that is worth adding
> to the PETSc source, or should I just write and use my own wrapper.
Just write and use your own wrapper.
Barry
>
>
> On Wed, Feb 27, 2013 at 4:51 PM, Matthe
1) Work with petsc-dev http://www.mcs.anl.gov/petsc/developers/index.html
the threaded stuff is more mature.
2) You should not need petscthreadcomm.h in Fortran. Simply using OpenMP
progmas in your portion of the code.
3) See http://www.mcs.anl.gov/petsc/features/threads.html for a
Hi, -- I'm new here
and I have a somehow seemingly simple task:
We need to get an exact copy of a previously set up SNES including all
function pointers, all options and especially the DM (also in deep copy
-- except the ApplicationContext struct).
is there a showstopper -- or maybe already fu
than any results to which their
experiments lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/3df881fa/attachment.html>
nters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/5a3c5cc8/attachment.html>
g a
> function VecCopyGhosted()?
>
> Thanks,
>
> Matt
>
>
>> John
>>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/b9860df2/attachment.html>
Hi Matt,
Yes. Thank you.
On Thu, Feb 28, 2013 at 6:53 AM, Matthew Knepley wrote:
> On Thu, Feb 28, 2013 at 8:24 AM, Garnet Vaz wrote:
>>
>> Hi,
>>
>> I would like to know what some of the arguments for
>> DMPlexCreateFromCellList are supposed to be for a triangular mesh.
>>
>> The function is:
s of the desired semantics.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/cc9ccd8c/attachment.html>
> -- Norbert Wiener
>
>
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/5650b1dc/attachment.html>
Hi,
I would like to know what some of the arguments for
DMPlexCreateFromCellList are supposed to be for a triangular mesh.
The function is:
DMPlexCreateFromCellList(MPI_Comm comm, PetscInt dim, PetscInt
numCells, PetscInt numVertices, PetscInt numCorners, PetscBool
interpolate, const int cells[],
d...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130228/f5de9578/attachment-0001.html>
18 matches
Mail list logo