.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_reliability_sensitivity/sensitivity_analysis/plot_functional_chaos_sensitivity.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_reliability_sensitivity_sensitivity_analysis_plot_functional_chaos_sensitivity.py: Sobol' sensitivity indices from chaos ===================================== .. GENERATED FROM PYTHON SOURCE LINES 6-23 In this example we are going to compute global sensitivity indices from a functional chaos decomposition. We study the Borehole function that models water flow through a borehole: .. math:: \frac{2 \pi T_u (H_u - H_l)}{\ln{r/r_w}(1+\frac{2 L T_u}{\ln(r/r_w) r^2_w K_w}\frac{T_u}{T_l})} With parameters: - :math:`r_w`: radius of borehole (m) - :math:`r`: radius of influence (m) - :math:`T_u`: transmissivity of upper aquifer (:math:`m^2/yr`) - :math:`H_u`: potentiometric head of upper aquifer (m) - :math:`T_l`: transmissivity of lower aquifer (:math:`m^2/yr`) - :math:`H_l`: potentiometric head of lower aquifer (m) - :math:`L`: length of borehole (m) - :math:`K_w`: hydraulic conductivity of borehole (:math:`m/yr`) .. GENERATED FROM PYTHON SOURCE LINES 25-32 .. code-block:: default from __future__ import print_function import openturns as ot from operator import itemgetter import openturns.viewer as viewer from matplotlib import pylab as plt ot.Log.Show(ot.Log.NONE) .. GENERATED FROM PYTHON SOURCE LINES 33-34 borehole model .. GENERATED FROM PYTHON SOURCE LINES 34-49 .. code-block:: default dimension = 8 input_names = ['rw', 'r', 'Tu', 'Hu', 'Tl', 'Hl', 'L', 'Kw'] model = ot.SymbolicFunction(input_names, ['(2*pi_*Tu*(Hu-Hl))/(ln(r/rw)*(1+(2*L*Tu)/(ln(r/rw)*rw^2*Kw)+Tu/Tl))']) coll = [ot.Normal(0.1, 0.0161812), ot.LogNormal(7.71, 1.0056), ot.Uniform(63070.0, 115600.0), ot.Uniform(990.0, 1110.0), ot.Uniform(63.1, 116.0), ot.Uniform(700.0, 820.0), ot.Uniform(1120.0, 1680.0), ot.Uniform(9855.0, 12045.0)] distribution = ot.ComposedDistribution(coll) distribution.setDescription(input_names) .. GENERATED FROM PYTHON SOURCE LINES 50-51 Freeze r, Tu, Tl from model to go faster .. GENERATED FROM PYTHON SOURCE LINES 51-60 .. code-block:: default selection = [1, 2, 4] complement = ot.Indices(selection).complement(dimension) distribution = distribution.getMarginal(complement) model = ot.ParametricFunction( model, selection, distribution.getMarginal(selection).getMean()) input_names_copy = list(input_names) input_names = itemgetter(*complement)(input_names) dimension = len(complement) .. GENERATED FROM PYTHON SOURCE LINES 61-62 design of experiment .. GENERATED FROM PYTHON SOURCE LINES 62-66 .. code-block:: default size = 1000 X = distribution.getSample(size) Y = model(X) .. GENERATED FROM PYTHON SOURCE LINES 67-68 create a functional chaos model .. GENERATED FROM PYTHON SOURCE LINES 68-74 .. code-block:: default algo = ot.FunctionalChaosAlgorithm(X, Y) algo.run() result = algo.getResult() print(result.getResiduals()) print(result.getRelativeErrors()) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none [0.00472648] [4.13378e-08] .. GENERATED FROM PYTHON SOURCE LINES 75-76 Quick summary of sensitivity analysis .. GENERATED FROM PYTHON SOURCE LINES 76-79 .. code-block:: default sensitivityAnalysis = ot.FunctionalChaosSobolIndices(result) print(sensitivityAnalysis.summary()) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none input dimension: 5 output dimension: 1 basis size: 56 mean: [74.8252] std-dev: [28.86] ------------------------------------------------------------ Index | Multi-indice | Part of variance ------------------------------------------------------------ 1 | [1,0,0,0,0] | 0.667533 3 | [0,0,1,0,0] | 0.0972783 4 | [0,0,0,1,0] | 0.0906816 2 | [0,1,0,0,0] | 0.0798009 5 | [0,0,0,0,1] | 0.0221444 ------------------------------------------------------------ ------------------------------------------------------------ Component | Sobol index | Sobol total index ------------------------------------------------------------ 0 | 0.67653 | 0.705682 1 | 0.0798379 | 0.0891789 2 | 0.0972796 | 0.108678 3 | 0.0916856 | 0.103876 4 | 0.0221444 | 0.0254584 ------------------------------------------------------------ .. GENERATED FROM PYTHON SOURCE LINES 80-81 draw Sobol' indices .. GENERATED FROM PYTHON SOURCE LINES 81-88 .. code-block:: default first_order = [sensitivityAnalysis.getSobolIndex(i) for i in range(dimension)] total_order = [sensitivityAnalysis.getSobolTotalIndex( i) for i in range(dimension)] graph = ot.SobolIndicesAlgorithm.DrawSobolIndices( input_names, first_order, total_order) view = viewer.View(graph) .. image-sg:: /auto_reliability_sensitivity/sensitivity_analysis/images/sphx_glr_plot_functional_chaos_sensitivity_001.png :alt: Sobol' indices :srcset: /auto_reliability_sensitivity/sensitivity_analysis/images/sphx_glr_plot_functional_chaos_sensitivity_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 89-91 We saw that total order indices are close to first order, so the higher order indices must be all quite close to 0 .. GENERATED FROM PYTHON SOURCE LINES 91-97 .. code-block:: default for i in range(dimension): for j in range(i): print(input_names[i] + ' & ' + input_names[j], ":", sensitivityAnalysis.getSobolIndex([i, j])) plt.show() .. rst-class:: sphx-glr-script-out Out: .. code-block:: none Hu & rw : 0.007773952603684 Hl & rw : 0.009536145155773607 Hl & Hu : 5.4681169248072785e-05 L & rw : 0.00917958869007991 L & Hu : 0.0010928197895924347 L & Hl : 0.001315360782351903 Kw & rw : 0.002320368250935267 Kw & Hu : 0.00027291064768849286 Kw & Hl : 0.00031801092178630964 Kw & L : 0.0003073000271333149 .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 2.739 seconds) .. _sphx_glr_download_auto_reliability_sensitivity_sensitivity_analysis_plot_functional_chaos_sensitivity.py: .. only :: html .. container:: sphx-glr-footer :class: sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_functional_chaos_sensitivity.py ` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_functional_chaos_sensitivity.ipynb ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_