Issue104 fix triangular + start bsplines deg = 2 on triangles

[Didou09]

Main changes:

• Triangular grid:

  • add_mesh_2d_tri(crop_poly=np.ndarray) implemented !
    • properly handles removing non-needed mesh elements
  • shape-c and shape-k replaced by shape_c and shape_k
  • Fixed testing on triangular grid trifinder (now used point on grid)

• Bsplines of deg = 2 on triangular grid:

  • add_bsplines(deg=2) operational
  • TODO: interpolation (evaluation in class)

TODO: cropbs

Issues:

• Fixes, in devel, issue Fix triangular mesh #104
• Partially fixes, in devel, issue bsplines of degree 2 on triangular meshes #99 => not closing, to be finished on branch 099

Example:

ITER triangular mesh from points cloud density and Delaunay triangulation + cropping

import bsplines2d as bs2

pfe = "C:/Users/dvezinet/Documents/FromOthers/Pedro_CARVHALO/point_cloud.csv"

out = np.loadtxt(pfe, skiprows=1, delimiter=',')

conf = tf.load_config('ITER-V0')

coll = bs2.BSplines2D()

coll.add_mesh_2d_tri(pts_x0=out[:, 0], pts_x1=out[:, 1], crop_poly=conf)

dax = coll.plot_mesh()

And adding bsplines of deg=2:

coll.add_bsplines(deg=2)

coll
 
ref key      size  nb. data  nb. data monot.
-----------  ----  --------  ---------------
m00-nk       1363  2         0              
m00-nc       2615  3         0              
m00-nind     3     1         0              
m00_bs0_nbs  2615  2         0              
m00_bs1_nbs  1363  2         0              
m00_bs2_nbs  2315  2         0              

data         shape      ref                     dim      quant    name  units  bsplines    
-----------  ---------  ----------------------  -------  -------  ----  -----  ------------
m00_k0       (1363,)    ('m00-nk',)                                                        
m00_k1       (1363,)    ('m00-nk',)                                                        
m00_c0       (2615,)    ('m00-nc',)                                                        
m00_c1       (2615,)    ('m00-nc',)                                                        
m00-ind      (2615, 3)  ('m00-nc', 'm00-nind')  indices  indices                           
m00_bs0_ap0  (2615,)    ('m00_bs0_nbs',)                                       ('m00_bs0',)
m00_bs0_ap1  (2615,)    ('m00_bs0_nbs',)                                       ('m00_bs0',)
m00_bs1_ap0  (1363,)    ('m00_bs1_nbs',)                                       ('m00_bs1',)
m00_bs1_ap1  (1363,)    ('m00_bs1_nbs',)                                       ('m00_bs1',)
m00_bs2_ap0  (2315,)    ('m00_bs2_nbs',)                                       ('m00_bs2',)
m00_bs2_ap1  (2315,)    ('m00_bs2_nbs',)                                       ('m00_bs2',)

mesh  nd  type  crop   ntri  knots                 cents                 shape_k  shape_c  ind      subkey  nb bsplines
----  --  ----  -----  ----  --------------------  --------------------  -------  -------  -------  ------  -----------
m00   2d  tri   False  1     ('m00_k0', 'm00_k1')  ('m00_c0', 'm00_c1')  (1363,)  (2615,)  m00-ind  None    3          

bsplines  deg  mesh  ref               ref_bs            apex                            shape    crop 
--------  ---  ----  ----------------  ----------------  ------------------------------  -------  -----
m00_bs0   0    m00   ('m00_bs0_nbs',)  ('m00_bs0_nbs',)  ('m00_bs0_ap0', 'm00_bs0_ap1')  (2615,)  False
m00_bs1   1    m00   ('m00_bs1_nbs',)  ('m00_bs1_nbs',)  ('m00_bs1_ap0', 'm00_bs1_ap1')  (1363,)  False
m00_bs2   2    m00   ('m00_bs2_nbs',)  ('m00_bs2_nbs',)  ('m00_bs2_ap0', 'm00_bs2_ap1')  (2315,)  False

point_cloud.csv