Skip to content

Commit

Permalink
Update t_FunctionalChaos_conditionalExpectation.py
Browse files Browse the repository at this point in the history
  • Loading branch information
mbaudin47 committed Sep 20, 2024
1 parent d81e4e7 commit e235c15
Showing 1 changed file with 4 additions and 3 deletions.
7 changes: 4 additions & 3 deletions python/test/t_FunctionalChaos_conditionalExpectation.py
Original file line number Diff line number Diff line change
@@ -1,6 +1,7 @@
#! /usr/bin/env python

import openturns as ot
import openturns.experimental as otexp
import math
from openturns.usecases import ishigami_function

Expand Down Expand Up @@ -41,7 +42,7 @@
print("absoluteTolerance = ", absoluteTolerance)
marginalDistribution = im.distributionX.getMarginal(0)
experiment = ot.MonteCarloExperiment(marginalDistribution, sampleSizeTest)
integration = ot.ExperimentIntegration(experiment)
integration = otexp.ExperimentIntegration(experiment)
error = integration.computeL2Norm(conditionalPCEFunction - functionEgivenX1Exact)
print("PCE | X1, error = ", error)
assert error[0] < absoluteTolerance
Expand All @@ -53,7 +54,7 @@
functionEgivenX2Exact = ot.ParametricFunction(parametricEgivenX2, [1, 2], [im.a, im.b])
marginalDistribution = im.distributionX.getMarginal(1)
experiment = ot.MonteCarloExperiment(marginalDistribution, sampleSizeTest)
integration = ot.ExperimentIntegration(experiment)
integration = otexp.ExperimentIntegration(experiment)
error = integration.computeL2Norm(conditionalPCEFunction - functionEgivenX2Exact)
print("PCE | X2, error = ", error)
assert error[0] < absoluteTolerance
Expand All @@ -66,7 +67,7 @@
sampleSizeTest = 10000
marginalDistribution = im.distributionX.getMarginal(2)
experiment = ot.MonteCarloExperiment(marginalDistribution, sampleSizeTest)
integration = ot.ExperimentIntegration(experiment)
integration = otexp.ExperimentIntegration(experiment)
error = integration.computeL2Norm(conditionalPCEFunction - functionEgivenX3Exact)
print("PCE | X3, error = ", error)
assert error[0] < absoluteTolerance

0 comments on commit e235c15

Please sign in to comment.