[Users] boundary conditions with Llama

Miguel Zilhão miguel.zilhao.nogueira at tecnico.ulisboa.pt
Mon Jun 17 17:42:17 CDT 2019

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?


More information about the Users mailing list