Anyone knows why applyBC_3D() is commented out in applyBC(Matrix, ...) and applyBC(Vector, ...)? /Anders