The following implied loop produces a bad performance at compile time
haarWPhi2D = reshape([((((((2**(-j/2))*haarWaveletPhi((2**j)*(m*dx)-k))* & ((2**(-j/2))
