Numerical differentiation with python

Question

I have such a task:

You need to write a function with the following signature: def derForward(f, I, h) This function should get as input function and segment as Python list of two elements --- ends of the segment. In the function you are asked to divide the segment into small segments of length ℎ , thus getting a grid . Your function should return dy --- forward differences for each point x (except the border since the formula asks for the next value). You should return both x and dy arrays of the same length.

My answer looks this way:

xx = [] #list of x values
frw = [] #list of forw. differencies of x

def derForward(f, I, h):
    x = np.arange(I[0], I[1], h)
    f = np.vectorize(f)
    for x_ in x[:-1]:
        dxdy = (f(x_+h)-f(x_))/ h
        xx.append(x_)
        frw.append(dxdy)
    x = np.array(xx)
    dy = np.array(frw)
    return x, dy

This function has to pass 3 tests for errors and accuracy. But now it is failing with tests 2 and 3. And I don't understand why I can pass only 1 test. The checker looks like this:

import numpy as np
from math import *


from time import time

def findif_check(derForward):
    count=0
    I=[0.001, 2*np.pi]
    f=lambda x: sin(x)
    h=0.01
    st=time()
    x, dy=derForward(f, I, h)
    dur=time()-st

    if x.shape[0]!=dy.shape[0]:
        print('FAILED: x and dy shape mismatch!')
    else:
        df=lambda x: cos(x)    
        err=np.max(np.abs(dy-np.vectorize(df)(x)))
        print('Test 1 |::|  err=', np.max(np.abs(dy-np.cos(x))), '  |::|  time: ', dur, 's')
        if err<2*0.0075:
            count+=1
            print('Test 1 |::|  accuracy OK') 
        else:
            print('Test 1 |::|  accuracy FAIL') 
        f=lambda x: x**x
        I=[0.001, 1]

        st=time()
        x, dy=derForward(f, I, h)
        dur=time()-st
        
        df=lambda x: x**x*(log(x)+1)
        err=np.max(np.abs(dy-np.vectorize(df)(x)))
        print('Test 2 |::|  err=', err, '  |::|  time: ', dur, 's')
        if err<2:
            count+=1
            print('Test 2 |::|  accuracy OK')
        else:
            print('Test 2 |::|  accuracy FAIL')

        f=lambda x: e**(-x**2)
        I=[0.001, 1]

        st=time()
        x, dy=derForward(f, I, h)
        dur=time()-st
        
        df=lambda x: -2*x*e**(-x**2)
        err=np.max(np.abs(dy-np.vectorize(df)(x)))
        print('Test 3 |::|  err=', err, '  |::|  time: ', dur, 's')
        if err<2*0.01:
            count+=1
            print('Test 3 |::|  accuracy OK')
        else:
            print('Test 3 |::|  accuracy FAIL')
    print('Passed: ', count, '/ 3')

cg-zhou · Answer 1 · 2021-12-29T04:41:00.093

Because when the test 2 and 3 is running, variable xx and frw contains data of test 1.

So move the definition of xx and frw into def derForward(f, I, h): as below:

def derForward(f, I, h):
    xx = [] #list of x values
    frw = [] #list of forw. differencies of x

    x = np.arange(I[0], I[1], h)
    f = np.vectorize(f)
    // ... ... ... ...

Then all 3 tests passed:

Test 1 |::|  err= 0.004999976659897233   |::|  time:  0.04198646545410156 s
Test 1 |::|  accuracy OK
Test 2 |::|  err= 1.715675908387758   |::|  time:  0.007221221923828125 s
Test 2 |::|  accuracy OK
Test 3 |::|  err= 0.009999270029526776   |::|  time:  0.00669407844543457 s
Test 3 |::|  accuracy OK
Passed:  3 / 3

Numerical differentiation with python

1 Answers1