Restructured test code, added LU solve for Helmholtz matrix
Bug fix: incorrect access descriptor for function u in LU solve

The function is modified in-place, so it should be RW instead of READ