Contact
CoCalc Logo Icon
StoreFeaturesDocsShareSupport News AboutSign UpSign In
| Download

Interpolation basis

Project: Interpolation
Path: jupyter / src / dbg.py
Views: 887
1
from sympy import diff

2


3
def dbg(show=False,**part):

4
    '''Work like:

5
    dbg(**{"string1":value,...})

6
    '''

7
    if show:

8
        if "lagrangeInterp|part1" in part:

9
            expr=part["lagrangeInterp|part1"]

10
            print("expr=%s"%(expr))

11
            print("diff(expr)=%s"%(diff(expr)))

12
        if "lagrangeInterp|part2" in part:

13
            n,a,d0=part["lagrangeInterp|part2"]

14
            print("len(c[0][:-1])=%s"%(n))

15
            print("a=%s"%(a))

16
            print("d[0]=%s"%(d0))

17
        if "lagrangeInterp|part3" in part:

18
            i,j,cij,aij,deltaAij=part["lagrangeInterp|part3"]

19
            #print('c[%s::%s][:-1]=%s'%(i,j,cij),end=" | ")

20
            #print('a[%s::%s]=%s'%(i,j,aij))

21
            #print('a[i::j][1:]-a[i::j][:-1]=%s'%(deltaAij))

22
        if "lagrangeInterp|part4" in part:

23
            j,dij1,dij2,dij3=part["lagrangeInterp|part4"]

24
            print("d[%s-1][:-%s]=%s"%(j,j,dij1))

25
            print("d[%s-2][:-(%s-1)]=%s"%(j,j,dij2))

26
            print("d[%s-2][:-%s]=%s"%(j,j,dij3))

27
        if "lagrangeInterp|part5" in part:

28
            a,c,d=part["lagrangeInterp|part5"]

29
            print("a=%s"%a)

30
            print("c=%s"%c)

31
            print('d=%s'%d)

32
        if "lagrangeInterp|part6" in part:

33
            A,X=part["lagrangeInterp|part6"]

34
            print("A=%s"%(A))

35
            print("X=%s"%(X))

36
        if "lagrangeInterp|part7" in part:

37
            fa0,exprSum=part["lagrangeInterp|part7"]

38
            print("f(a[0])=%s"%(fa0))

39
            print("exprSum=%s"%(exprSum))

40
    

41
        if "diffSocket|part1" in part:

42
            a,points=part["diffSocket|part1"]

43
            print("a=%s"%(a))

44
            print("points=%s"%(points))

45
        if "diffSocket|part21" in part:

46
            arg0,arg1,arg2,arg3,arg4=part["diffSocket|part21"]

47
            print("b[:-1][%s]=a[:-1][%s]+(x0-points[%s])*b[(%s)]=%s"%(arg0,arg1,arg2,arg3,arg4))

48
        if "diffSocket|part22" in part:

49
            points=part["diffSocket|part22"]

50
            print("points=%s"%(points))

51
        if "diffSocket|part3" in part:

52
            b,deriv=part["diffSocket|part3"]

53
            print("b=%s"%(b))

54
            print("derivations=%s"%deriv)

55
        if "ifEqPoints|part1" in part:

56
            j,ci,di=part["ifEqPoints|part1"]

57
            print('any c[%s-1]=%s'%(j,ci))

58
            print("d[%s-1][:]=%s"%(j,di))

59
        if "ifEqPoints|part21" in part:

60
            j,i,dji=part["ifEqPoints|part21"]

61
            print("if d[%s-1][:-%s][%s]=%s"%(j,j,i,dji))

62
        if "ifEqPoints|part22" in part:    

63
            j,i,dji=part["ifEqPoints|part22"]

64
            print("else d[%s-1][:-%s][%s]=%s"%(j,j,i,dji))

65
        if "leastSquares|part1" in part:    

66
            a,b=part["leastSquares|part1"]

67
            print('a=%s'%a)

68
            print('b=%s'%b)

69
        if "leastSquares|part2" in part:    

70
            HH,b,a,expr=part["leastSquares|part2"]

71
            print("************")

72
            print('HH=%s'%(HH))

73
            print("************")

74
            print("**********")

75
            print('b=%s'%(b))

76
            print('a=%s'%(a))

77
            print(expr)

78
        

79
        if "findParabs|part1" in part:    

80
            a0,a1,b0,b1=part["findParabs|part1"]

81
            print("a0=%s|a1=%s|b0=%s|b1=%s"%(a0,a1,b0,b1))

82


83
        if "findParabs|part2" in part:    

84
            b,invC,p0,p1,c0,c1=part["findParabs|part2"]

85
            print("b=%s"%(b))

86
            print("invC=%s"%(invC))

87
            print("p0,p1= %s,%s"%(p0,p1))

88
            print("c0,c1=%s %s"%(c0,c1))

89


90