Hi,
I am using OpenTURNS to get Sobol' indices.
I have a function which doesn't return a scalar value but a vector (1D). With
ot.SensitivityAnalysis(), I am able to get the indices by component.
I just iterate over:
output_len = 400
first = [ ]
sa = ot.SensitivityAnalysis(...)
for i in range(output_len):
first.append(sa.getFirstOrderIndices(i))
If I replace the method by ot.FAST I don't get the right output. Indices are
constants.
* I am thinking, maybe ot.FAST doesn't allow to do that?
* Another question, is it possible to set up a model with a function which
has a 2D (or even more) output field?
I would like then to do something like that:
output_len = 400
first = np.empty([output_len, output_len])
sa = ot.SensitivityAnalysis(...)
for i, j in itertools.product(range(output_len), range(output_len)):
first[i, j] = sa.getFirstOrderIndices(i, j)
Thank you for your help :)
Pamphile ROY
UQ Trainee at CERFACS
_______________________________________________
OpenTURNS users mailing list
[email protected]
http://openturns.org/mailman/listinfo/users
