fixed various bugs. Main issue was only doing boundary integral on external face with a dirichlet boundary condition.
fixed various bugs. Main issue was only doing boundary integral on external face with a dirichlet boundary condition.

now also solved with inverse of mass matrix

jjo31/ATHAM-Fluidity,jjo31/ATHAM-Fluidity,rjferrier/fluidity,rjferrier/fluidity,jjo31/ATHAM-Fluidity,iakovos-panourgias/fluidity,jrper/fluidity,rjferrier/fluidity,iakovos-panourgias/fluidity,jrper/fluidity,jrper/fluidity,iakovos-panourgias/fluidity,jrper/fluidity,rjferrier/fluidity,rjferrier/fluidity,jjo31/ATHAM-Fluidity,jjo31/ATHAM-Fluidity,iakovos-panourgias/fluidity,iakovos-panourgias/fluidity,jrper/fluidity