tests/xn/curve.rst

   1 >>> import ion.xn.curve as amc
   2 >>> amc._tile ([0,1,2,3], [40, 64, -1, 4], 1, (4,6))
   3 ([0, 1, 2, 3, 4, 5, 6], [40, 64, -1, 4, 40, 64, -1])
   4
   5 >>> amc._tile ([0,1,2,3], [40, 64, -1, 4], -1, (-3, -1))
   6 ([-3, -2, -1, 0, 1, 2, 3], [64, -1, 4, 40, 64, -1, 4])
   7
   8 >>> amc._extend ([0,1,2,3], [40, 64, -1, 4], 1, (4,6))
   9 [0,1,2,3,4,5,6], [40, 64, -1, 4, 44, 68, 3]
  10
  11 >>> amc._extend ([0,1,2,3], [40, 64, -1, 4], -1, (-3, -1))
  12 [-3, -2, -1, 0, 1, 2, 3], [24, -41, -36, 40, 64, -1, 4]
  13
  14
  15 # XXX instead of expanding the spline, wrap the requested X values and add an appropriate offset.
  16 # eg for tile with a spline ranging 0..n values of X, you would just return spl(x % max(spl))
  17 # for min(spl) != 0, normalization must also be performed : spl (((x - min(spl)) % max(spl)) + min(spl))
  18
  19 >>> ms = amc.multiSpline ([0,.25,.50,.75,1.], [[0, 0, 0], [.4, .05, .1], [.5, .25, .3],
  20 ... [1., 0.95, .6], [1., 1., 1.]], k = 3)
  21
  22 >>> import numpy
  23 >>> numpy.round(ms(numpy.arange(0,1.01, 1./8)), 2)
  24 array([[-0.  , -0.  ,  0.  ],
  25        [ 0.31,  0.05,  0.04],
  26        [ 0.4 ,  0.05,  0.1 ],
  27        [ 0.42,  0.09,  0.19],
  28        [ 0.5 ,  0.25,  0.3 ],
  29        [ 0.73,  0.59,  0.44],
  30        [ 1.  ,  0.95,  0.6 ],
  31        [ 1.14,  1.15,  0.79],
  32        [ 1.  ,  1.  ,  1.  ]])
  33
  34 Notice the overflow above: >1 is not legal.
  35 fixing that:
  36
  37 >>> ms = amc.multiSpline ([0,.25,.50,.75,1.], [[0, 0, 0], [.4, .05, .1], [.5, .25, .3],
  38 ... [1., 0.95, .6], [1., 1., 1.]], k = 3, subdiv = True, protected = ( (0,1), (1,2), (2,3),))
  39
  40 >>> numpy.round(ms(numpy.arange(0,1.01, 1./8)), 2)
  41 array([[-0.  , -0.  ,  0.  ],
  42        [ 0.31,  0.05,  0.04],
  43        [ 0.4 ,  0.05,  0.1 ],
  44        [ 0.42,  0.09,  0.19],
  45        [ 0.5 ,  0.25,  0.3 ],
  46        [ 0.73,  0.59,  0.44],
  47        [ 1.  ,  0.95,  0.6 ],
  48        [ 1.  ,  0.97,  0.79],
  49        [ 1.  ,  1.  ,  1.  ]])