thank you for your fast reply!
With the solution you suggested, I have no more multiple models.
But unfortunately the problem with AKMCS persists.
I have to be clearer with the description of the .m file model which I called
myFunction. In this function also another function is used, which creates the .inp file to be read by the FE software. This because I pass to my FE software a lognormal random field. So
myFunction takes X of size d, where there are (d-1) normal Gaussian random values and another random value which follows a lognormal distribution. Then a function, say
X(1:d-1) and creates the random field. The random field is then passed to a function called, say,
INPGenerator, which creates the .inp file to be read and evaluated from my FE software thanks to
INPGenerator not only sets at the “right place” the values of the random field, but also sets the string
<X0001> at the place where the random variable which follows a lognormal distribution (
X(d)) is going to be inserted. Here how my code now looks like:
function Y = myFunction(X,myUQLinkModel)
RF = RFGenerator(X(1:d-1)); % Creates the random field
INPGenerator(RF); % Creates the file .inp
FEMresult = uq_evalModel(myUQLinkModel,X(d)); % Evaluates the .inp with the FE software
Y = f(FEMresult); % Some calculations to get the output of interest
The fact, that I create the file .inp, could cause the problem with AKMCS?