Hi Jan
If building your own trilinear interpolation engine is a feasible option (the math is computationally inexpensive), you can find a brief introduction via Paul Bourke's web presence.
http://astronomy.swin.edu.au/~pbourke/analysis/trilinear/
Depending on the type of colour transformation you wish to employ in your 3D LUT (the difference between adjacent response vectors) & flexibility/precision you have with that data (ie. sampled empirically or generated procedurally) you may find it more appropriate to reduce the sampling frequency in each dimension (less memory) & employ a higher order interpolation algorithm (eg. tricubic).
What type of "complex colour transformations" are you intending to use the plugin for?
Dominic Glynn