Multidimensional Interpolation with knots condition?

General Usage
,
1

Hi Everyone!
I am trying to make up an interpolated multidimansional function with 5 dimensions and 1 variable. In this case:

X1 X2 X3 X4 X5 Density

X1=Tempertaure
X2=Preassure
X3 to X5: % of composition of 3 elements
Density.

So, in my data table I have something like:

   #X1      X2     X3    X4     X5    Density1 Density2 Density3 
 200.0    0.5    0.0   70.0   30.0  6.2762  3.2896  2.6654
 400.0    0.2    0.0   70.0   30.0  6.2085  3.2603  2.643
 1000.0  0.6    0.0   90.0   10.0  6.3092  3.273   2.6636

FULL DATAFILE at: SWF/Data.txt at main · marianoarnaiz/SWF · GitHub

So X1 and X2 are sampled


X1=[200 400 600 800 1000]
X2=[0.0 0.1 0.2 0.3 0.4 0.5 0.6]

But X3 to X5 have a condition, they must sum = 100.

So:
[215 0.05 75 12.5 12.5] #Should exist from the interpolation

But:
[215 0.05 100 10 10] #Should NOT exist from the interpolation.

When I am building my interpolation function:


knots = ([t for t in X1], [p for p in X2], [qz for qz in X3], [cl for cl in X4], [cy for cy in X5]);
Density = interpolate(knots, D_M, Gridded(Linear()));

I need to create a D_M (Density multidimensional array). I have tried filling it with missing, NaN, Big values… It all gives me problems because I can only evaluate in the nodes with data and Interpolation in a none knot gives me a problem.

Any ideas?