Optimizing volumetric functions is usually about getting the inner loop and the execution overheard to be a low as possible. That often seems to just require lots of experimentation to get right. If you post your shader code for it, it may be possible to optimize a bit. And it also may be easier to control how it performs using a regular material with the custom node. You may need to tweak the raymarching thresholds more aggressively etc.