sandbox/bugs/mixed.c

Boundary functions do not understand u.n, u.t yet. A workaround for this example is to replace u.t with _s in the RHS.

vector u[];
double bi = 2.;

u.t[bottom] = u.t[]*(2. - bi*Δ)/(2. + bi*Δ);
u.t[bottom] = val(_s,0)*(2. - bi*Δ)/(2. + bi*Δ); // this is the workaround

int main() {}