On Sun, Sep 10, 2017 at 5:51 PM, David Gross wrote:
> Hello,
> I was wondering if there was a matrix equivalent to the vecDot function
> (Frobenius inner product)? As far as I can tell the closest thing is
> MatNorm with NORM_FROBENIUS, but obviously this is acting on only
Hello,
I was wondering if there was a matrix equivalent to the vecDot function
(Frobenius inner product)? As far as I can tell the closest thing is
MatNorm with NORM_FROBENIUS, but obviously this is acting on only one
matrix.
If there is not a built in function, what is the best way to compute
Matthew Knepley writes:
> On Sun, Sep 10, 2017 at 9:15 AM, Praveen C wrote:
>
>>
>>
>>>
>>> You can use SNESSetJacobian, but DMSNESSetJacobianLocal is more
>>> convenient in my opinion.
>>>
>>> > I still have some doubt. If
>>> >
>>> > da0 --> for residual
On Sun, Sep 10, 2017 at 9:15 AM, Praveen C wrote:
>
>
>>
>> You can use SNESSetJacobian, but DMSNESSetJacobianLocal is more
>> convenient in my opinion.
>>
>> > I still have some doubt. If
>> >
>> > da0 --> for residual
>> > da1 --> for jacobian
>> >
>> > I still need to call
>
>
> You can use SNESSetJacobian, but DMSNESSetJacobianLocal is more
> convenient in my opinion.
>
> > I still have some doubt. If
> >
> > da0 --> for residual
> > da1 --> for jacobian
> >
> > I still need to call SNESetDM. So here I pass da0 ? But then how will
> snes
> > know about da1 ?
>
>
Praveen C writes:
> Hello Matt
>
>
>> No. You want 2 DAs, 1 for evaluating residuals with s = 3, and one for
>> making a Jacobian with s = 1. You give these
>> two DAs to
>>
>> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/
>> DMDASNESSetFunctionLocal.html
Praveen C writes:
> Dear all
>
> I am solving 3 coupled first order PDE in 2-D with SNES. The stencil is
> (3+1+3) in both directions due to weno scheme. I want to use matrix-free
> gmres. For preconditioner, I will use first order scheme with stencil
> (1+1+1) in both
Hello Matt
> No. You want 2 DAs, 1 for evaluating residuals with s = 3, and one for
> making a Jacobian with s = 1. You give these
> two DAs to
>
> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/
> DMDASNESSetFunctionLocal.html
>
On Sun, Sep 10, 2017 at 2:39 AM, Praveen C wrote:
> Dear all
>
> I am solving 3 coupled first order PDE in 2-D with SNES. The stencil is
> (3+1+3) in both directions due to weno scheme. I want to use matrix-free
> gmres. For preconditioner, I will use first order scheme with
Dear all
I am solving 3 coupled first order PDE in 2-D with SNES. The stencil is
(3+1+3) in both directions due to weno scheme. I want to use matrix-free
gmres. For preconditioner, I will use first order scheme with stencil
(1+1+1) in both directions. I am using DMDA.
Two things that I am
10 matches
Mail list logo