n1 = p1 + n_elements_1 k1 = p1 + 1 verbose = False #verbose = True # ... # ... [u1,w1] = legendre(p1) # ... # ... m1 = n1 + p1 + 1 knots1 = zeros(m1, double) # call to spl knots1 = spl_make_open_knots (n1, p1) # ... # ... TODO fix args of zeros m1 = n_elements_1+1 grid_1 = zeros(m1, double) # call to spl grid_1 = spl_construct_grid_from_knots(p1, n1, n_elements_1, knots1) # ... # ... construct the quadrature points grid points_1 = zeros((k1, n_elements_1), double) weights_1 = zeros((k1, n_elements_1), double) # call to spl
[u1, w1] = legendre(p1) # ... # ... [u2, w2] = legendre(p2) # ... # ... m1 = n1 + p1 + 1 m2 = n2 + p2 + 1 knots1 = zeros(m1, double) knots2 = zeros(m2, double) # call to spl knots1 = spl_make_open_knots(n1, p1) # call to spl knots2 = spl_make_open_knots(n2, p2) # ... # ... TODO fix args of zeros m1 = n_elements_1 + 1 m2 = n_elements_2 + 1 grid_1 = zeros(m1, double) grid_2 = zeros(m2, double) # call to spl grid_1 = spl_construct_grid_from_knots(p1, n1, n_elements_1, knots1)