It seems that the evaluation of shape functions is not thread safe. Is this something that can be fixed easily?