#! /usr/bin/env python

from openturns import *

TESTPREAMBLE()

try:
    f = NumericalMathFunction("x", "x^2", "y")
    size = 4
    inputSample = NumericalSample(size, 1)
    for i in range(size):
        inputSample[i, 0] = i
    outputSample = f(inputSample)
    print "Is history enabled for f? ", f.isHistoryEnabled()
    print "input history=", f.getHistoryInput()
    print "output history=", f.getHistoryOutput()
    f.enableHistory()
    outputSample = f(inputSample)
    print "Is history enabled for f? ", f.isHistoryEnabled()
    print "input history=", f.getHistoryInput()
    print "output history=", f.getHistoryOutput()
    f.clearHistory()
    print "Is history enabled for f? ", f.isHistoryEnabled()
    print "input history=", f.getHistoryInput()
    print "output history=", f.getHistoryOutput()
    # Perform the computation twice
    outputSample = f(inputSample)
    outputSample = f(inputSample)
    print "input history=", f.getHistoryInput()
    print "output history=", f.getHistoryOutput()

except:
    import sys
    print "t_NumericalMathFunction_history.py", sys.exc_type, sys.exc_value
