Hello Miguel, assuming that your outer boundary is spherical, your should set:
NewRad::z_is_radial = "yes" Yours, Roland > 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? > > thanks, > Miguel > _______________________________________________ > Users mailing list > [email protected] > http://lists.einsteintoolkit.org/mailman/listinfo/users -- 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.
pgp6nu1wa59P2.pgp
Description: OpenPGP digital signature
_______________________________________________ Users mailing list [email protected] http://lists.einsteintoolkit.org/mailman/listinfo/users
