'pdetool style' assembling of boundary conditions
[Q,G,H,R,F]=gf_asm_pdetoolbc(mf_u, mf_d, b, e, f_expr)