> No, A(t) has five diagonals as indicated by the calls to MatSetValues()
with five entries per row; they correspond to the usual 5 pt stencil in 2d.
This might not be the right place to ask this, but I'm new to PDE's and I
expected to find something like the matrix to be used as the one derived
> On Dec 22, 2018, at 12:29 PM, Matthew Knepley via petsc-users
> wrote:
>
> On Sat, Dec 22, 2018 at 1:18 PM Sajid Ali via petsc-users
> wrote:
> Hi,
>
> I have a few questions about ex13 from TS explaining how to solve the
> following PDE in 2D :
> u_t = uxx + uyy
> 1) Can I get a
Hi,
I have a few questions about ex13 from TS explaining how to solve the
following PDE in 2D :
u_t = uxx + uyy
1) Can I get a reference/explanation for how the RHSJacobian function is
derived ?
2) If the problem is linear, why is TSSetProblemType never set to TS_LINEAR?
3) Since this is just
Wow, this is an old thread.
Sorry if I sound like an old fart talking about the good old days but I
originally did RAP. in Prometheus, in a non work optimal way that might be
of interest. Not hard to implement. I bring this up because we continue to
struggle with this damn thing. I think this
OK, so this thread has drifted, see title :)
On Fri, Dec 21, 2018 at 10:01 PM Fande Kong wrote:
> Sorry, hit the wrong button.
>
>
>
> On Fri, Dec 21, 2018 at 7:56 PM Fande Kong wrote:
>
>>
>>
>> On Fri, Dec 21, 2018 at 9:44 AM Mark Adams wrote:
>>
>>> Also, you mentioned that you are using