[Users] boundary conditions with Llama
rhaas at illinois.edu
Mon Jun 17 20:45:17 CDT 2019
assuming that your outer boundary is spherical, your should set:
NewRad::z_is_radial = "yes"
> hi all,
> i've been implementing a scalar wave evolution code to work with the Llama multipatch
> infrastructure. i think i have it working except for the boundary conditions. as a first test, i
> tried to call the NewRad infrastructure, in the same way that i usually do for my non-multipatch codes:
> ierr = NewRad_Apply(cctkGH, phi, rhs_phi, phi_0, one, n_phi)
> ierr = NewRad_Apply(cctkGH, Kphi, rhs_Kphi, Kphi_0, one, n_Kphi)
> but this does not seem to work as i get instabilities very early on... am i missing something, or is
> NewRad not designed to work with Llama? if so, what is the "standard" way of imposing BCs for
> multipatch evolutions?
> Users mailing list
> Users at einsteintoolkit.org
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://keys.gnupg.net.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Size: 833 bytes
Desc: OpenPGP digital signature
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20190618/292b5a03/attachment.bin
More information about the Users