[0, M, 0],
>>
>> [H'(p1), H'(p2), \sigma*M]]
>>
>>
>>
>> where M is a mass matrix, H'(p1) is the jacobian of H(p1, p2) w.r.t. p1 and
>> H'(p2), the jacobian of H(p1, p2) w.r.t. p2. H'(p1) and H'(p2) are
>>
>>
>> where M is a mass matrix, H'(p1) is the jacobian of H(p1, p2) w.r.t. p1
>> and H'(p2), the jacobian of H(p1, p2) w.r.t. p2. H'(p1) and H'(p2) are
>> unnecessary for the solver strategy I want to implement.
>>
>>
>>
>> T
p2), the jacobian of H(p1, p2) w.r.t. p2. H'(p1) and H'(p2) are
> unnecessary for the solver strategy I want to implement.
>
>
>
> Thanks
>
> Miguel
>
>
>
>
>
>
>
> From: Barry Smith mailto:bsm...@petsc.dev>>
> Date: Monda
Ok I will investigate implementing it using SLATE, thanks.
Miguel
From: Matthew Knepley
Date: Tuesday, March 23, 2021 at 12:57 PM
To: "Salazar De Troya, Miguel"
Cc: Barry Smith , "Jorti, Zakariae via petsc-users"
Subject: Re: [petsc-users] Local Discontinuous Galerkin w
H'(p1) is the jacobian of H(p1, p2) w.r.t. p1
> and H'(p2), the jacobian of H(p1, p2) w.r.t. p2. H'(p1) and H'(p2) are
> unnecessary for the solver strategy I want to implement.
>
>
>
> Thanks
>
> Miguel
>
>
>
>
>
>
>
> *From: *Ba
ssary
for the solver strategy I want to implement.
Thanks
Miguel
From: Barry Smith
Date: Monday, March 22, 2021 at 7:42 PM
To: Matthew Knepley
Cc: "Salazar De Troya, Miguel" , "Jorti, Zakariae via
petsc-users"
Subject: Re: [petsc-users] Local Discontinuous Galerkin with
u_t = G(u)
I don't see why you won't just compute any needed u_x from the given u and
then you can use any explicit or implicit TS solver trivially. For implicit
methods it can automatically compute the Jacobian of G for you or you can
provide it directly. Explicit methods will just use
On Mon, Mar 22, 2021 at 7:53 PM Salazar De Troya, Miguel via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hello
>
>
>
> I am interested in implementing the LDG method in “A local discontinuous
> Galerkin method for directly solving Hamilton–Jacobi equations”
> https://www.sciencedirect.com/scie
Hello
I am interested in implementing the LDG method in “A local discontinuous
Galerkin method for directly solving Hamilton–Jacobi equations”
https://www.sciencedirect.com/science/article/pii/S0021999110005255. The
equation is more or less of the form (for 1D case):
p1 = f(u_x)
