[petsc-users] Finite Element Method of 3d wave equation

[Jed Brown]

Please always use "reply-all" so that your messages go to the list.
This is standard mailing list etiquette.  It is important to preserve
threading for people who find this discussion later and so that we do
not waste our time re-answering the same questions that have already
been answered in private side-conversations.  You'll likely get an
answer faster that way too.

吕超 <luchao at mail.iggcas.ac.cn> writes:

> Hi Jed Brown:
>
>        Thank you for your reply. My Phd's thesis is about 3d wave
>        equation solved by spectral element method(high order finite
>        element method). So I want to kown whether you have some docs
>        or Petsc's programs about my thesis direction. I am so sorry to
>        cause you so much trouble. Thanks!

PETSc does not have special code for spectral element methods, though
there is at least one example that uses spectral elements
(src/ksp/ksp/examples/tutorials/ex59.c) and PetscFE (new in petsc-dev)
has support for arbitrary order.

The formulation in src/ts/examples/tutorials/ex12.c is a reasonable
formulation of the time-dependent problem, though it is common to ask
for a symplectic time integrator.  Gauss methods are one way to get
there (not necessarily the most efficient, but very high order and can
be vectorized--new research topic), though there is a large literature
on symplectic methods.  Additional time integration methods are easy to
add.

> -----原始邮件-----
> 发件人: "吕超" <luchao at mail.iggcas.ac.cn>
> 发送时间: 2013年12月26日 星期四
> 收件人: "Jed Brown" <jedbrown at mcs.anl.gov>
> 抄送:
> 主题: Re: Re: 3d acoustic wave equation
>
>
>
>> -----原始邮件-----
>> 发件人: "Jed Brown" <jedbrown at mcs.anl.gov>
>> 发送时间: 2013年12月26日 星期四
>> 收件人: "吕超" <luchao at mail.iggcas.ac.cn>
>> 抄送: petsc-users <petsc-users at mcs.anl.gov>
>> 主题: Re: 3d acoustic wave equation
>> 
>> 吕超 <luchao at mail.iggcas.ac.cn> writes:
>> 
>> > hi, Jed Brown:
>> >
>> >        Thank you for your reply.
>> >
>> >        I had met you in Changchun ,China, at about 2013/08. After
>> >        looking your operation of PETSC in your pc on that meeting, I
>> >        started to konw and study the Petsc. Do you have some Petsc's
>> >        program about the 3d acoustic wave equation(u_tt = ( c2 ) *
>> >        LAPLACIAN( u )). In Petsc's ts's tutorials, I just find the
>> >        heat equation(u_t = u_xx).
>> 
>> Look at src/ts/examples/tutorials/ex12.c
>> 
>> It could also be easily added as a first-order system tojj
>> src/ts/examples/tutorials/ex11.c
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20131226/97442269/attachment.pgp>