Lines Matching full:1
47 nx = nelem + 1
48 nu = nelem * (p - 1) + 1
65 indx[2 * i + 1] = i + 1
66 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
72 indu[p * i + j] = i * (p - 1) + j
73 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
75 strides = np.array([1, q, q], dtype="int32")
76 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
79 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
80 bu = ceed.BasisTensorH1Lagrange(1, 1, p, q, libceed.GAUSS)
86 qf_setup = ceed.QFunction(1, qfs.setup_mass,
88 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
89 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
90 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
92 qf_mass = ceed.QFunction(1, qfs.apply_mass,
94 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
95 qf_mass.add_input("u", 1, libceed.EVAL_INTERP)
96 qf_mass.add_output("v", 1, libceed.EVAL_INTERP)
136 nx = nelem + 1
137 nu = nelem * (p - 1) + 1
154 indx[2 * i + 1] = i + 1
155 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
161 indu[p * i + j] = i * (p - 1) + j
162 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
164 strides = np.array([1, q, q], dtype="int32")
165 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
168 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
169 bu = ceed.BasisTensorH1Lagrange(1, 1, p, q, libceed.GAUSS)
175 qf_setup = ceed.QFunction(1, qfs.setup_mass,
177 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
178 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
179 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
181 qf_mass = ceed.QFunction(1, qfs.apply_mass,
183 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
184 qf_mass.add_input("u", 1, libceed.EVAL_INTERP)
185 qf_mass.add_output("v", 1, libceed.EVAL_INTERP)
204 u.set_value(1.)
226 nx = nelem + 1
227 nu = nelem * (p - 1) + 1
244 indx[2 * i + 1] = i + 1
245 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
251 indu[p * i + j] = 2 * (i * (p - 1) + j)
252 ru = ceed.ElemRestriction(nelem, p, 2, 1, 2 * nu, indu,
254 strides = np.array([1, q, q], dtype="int32")
255 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
258 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
259 bu = ceed.BasisTensorH1Lagrange(1, 2, p, q, libceed.GAUSS)
265 qf_setup = ceed.QFunction(1, qfs.setup_mass,
267 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
268 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
269 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
271 qf_mass = ceed.QFunction(1, qfs.apply_mass_two,
273 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
297 u_array[2 * i] = 1.
298 u_array[2 * i + 1] = 2.
307 total_2 = total_2 + v_array[2 * i + 1]
323 nx = nelem + 1
324 nu = nelem * (p - 1) + 1
341 indx[2 * i + 1] = i + 1
342 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
348 indu[p * i + j] = i * (p - 1) + j
349 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
351 strides = np.array([1, q, q], dtype="int32")
352 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
355 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
356 bu = ceed.BasisTensorH1Lagrange(1, 1, p, q, libceed.GAUSS)
362 qf_setup = ceed.QFunction(1, qfs.setup_mass,
364 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
365 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
366 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
368 qf_mass = ceed.QFunction(1, qfs.apply_mass,
370 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
371 qf_mass.add_input("u", 1, libceed.EVAL_INTERP)
372 qf_mass.add_output("v", 1, libceed.EVAL_INTERP)
391 u.set_value(1)
412 nx = nelem + 1
413 nu = nelem * (p - 1) + 1
422 indx[2 * i + 1] = i + 1
423 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
429 indu[p * i + j] = i * (p - 1) + j
430 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
432 strides = np.array([1, q, q], dtype="int32")
433 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
436 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
437 bu = ceed.BasisTensorH1Lagrange(1, 1, p, q, libceed.GAUSS)
443 qf_setup = ceed.QFunction(1, qfs.setup_mass,
445 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
446 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
447 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
449 qf_mass = ceed.QFunction(1, qfs.apply_mass,
451 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
452 qf_mass.add_input("u", 1, libceed.EVAL_INTERP)
453 qf_mass.add_output("v", 1, libceed.EVAL_INTERP)
487 nx = nelem + 1
488 nu = nelem * (p - 1) + 1
505 indx[2 * i + 1] = i + 1
506 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
512 indu[p * i + j] = i * (p - 1) + j
513 ru = ceed.ElemRestriction(nelem, p, 1, 1, nu, indu,
515 strides = np.array([1, q, q], dtype="int32")
516 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
519 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
520 bu = ceed.BasisTensorH1Lagrange(1, 1, p, q, libceed.GAUSS)
526 qf_setup = ceed.QFunction(1, qfs.setup_mass,
528 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
529 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
530 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
532 qf_mass = ceed.QFunction(1, qfs.apply_mass,
534 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
535 qf_mass.add_input("u", 1, libceed.EVAL_INTERP)
536 qf_mass.add_output("v", 1, libceed.EVAL_INTERP)
555 u.set_value(1.)
567 v.set_value(1.)
590 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
597 x_array[i] = (1. / (nx * 2)) * (i % (nx * 2 + 1))
598 x_array[i + ndofs] = (1. / (ny * 2)) * (i / (nx * 2 + 1))
610 offset = col * 2 + row * (nx * 2 + 1) * 2
613 indx[i * 2 * p + 1] = 9 + offset
615 indx[i * 2 * p + 3] = 1 + offset
629 ru = ceed.ElemRestriction(nelem, p, 1, 1, ndofs, indx,
631 strides = np.array([1, q, q], dtype="int32")
632 rui = ceed.StridedElemRestriction(nelem, q, 1, nqpts, strides)
641 bu = ceed.BasisH1(libceed.TRIANGLE, 1, p, q, interp, grad, qref, qweight)
647 qf_setup = ceed.QFunction(1, qfs.setup_mass_2d,
649 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
651 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
653 qf_mass = ceed.QFunction(1, qfs.apply_mass,
655 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
656 qf_mass.add_input("u", 1, libceed.EVAL_INTERP)
657 qf_mass.add_output("v", 1, libceed.EVAL_INTERP)
697 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
704 x_array[i] = (1. / (nx * 2)) * (i % (nx * 2 + 1))
705 x_array[i + ndofs] = (1. / (ny * 2)) * (i / (nx * 2 + 1))
717 offset = col * 2 + row * (nx * 2 + 1) * 2
720 indx[i * 2 * p + 1] = 9 + offset
722 indx[i * 2 * p + 3] = 1 + offset
736 ru = ceed.ElemRestriction(nelem, p, 1, 1, ndofs, indx,
738 strides = np.array([1, q, q], dtype="int32")
739 rui = ceed.StridedElemRestriction(nelem, q, 1, nqpts, strides)
748 bu = ceed.BasisH1(libceed.TRIANGLE, 1, p, q, interp, grad, qref, qweight)
754 qf_setup = ceed.QFunction(1, qfs.setup_mass_2d,
756 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
758 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
760 qf_mass = ceed.QFunction(1, qfs.apply_mass,
762 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
763 qf_mass.add_input("u", 1, libceed.EVAL_INTERP)
764 qf_mass.add_output("v", 1, libceed.EVAL_INTERP)
783 u.set_value(1.)
805 nx_tet, ny_tet, nx_hex = 3, 1, 3
806 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
812 for i in range(ny * 2 + 1):
813 for j in range(nx * 2 + 1):
814 x_array[i + j * (ny * 2 + 1)] = i / (2 * ny)
815 x_array[i + j * (ny * 2 + 1) + ndofs] = j / (2 * nx)
830 offset = col * 2 + row * (nx * 2 + 1) * 2
833 indx_tet[i * 2 * p_tet + 1] = 9 + offset
835 indx_tet[i * 2 * p_tet + 3] = 1 + offset
849 ru_tet = ceed.ElemRestriction(nelem_tet, p_tet, 1, 1, ndofs, indx_tet,
851 strides = np.array([1, q_tet, q_tet], dtype="int32")
852 rui_tet = ceed.StridedElemRestriction(nelem_tet, q_tet, 1, nqpts_tet,
863 bu_tet = ceed.BasisH1(libceed.TRIANGLE, 1, p_tet, q_hex, interp, grad, qref,
870 qf_setup_tet = ceed.QFunction(1, qfs.setup_mass_2d,
872 qf_setup_tet.add_input("weights", 1, libceed.EVAL_WEIGHT)
874 qf_setup_tet.add_output("rho", 1, libceed.EVAL_NONE)
876 qf_mass_tet = ceed.QFunction(1, qfs.apply_mass,
878 qf_mass_tet.add_input("rho", 1, libceed.EVAL_NONE)
879 qf_mass_tet.add_input("u", 1, libceed.EVAL_INTERP)
880 qf_mass_tet.add_output("v", 1, libceed.EVAL_INTERP)
902 offset = (nx_tet * 2 + 1) * (ny_tet * 2) * (1 + row) + col * 2
907 k * (nx_hex * 2 + 1) + j
912 ru_hex = ceed.ElemRestriction(nelem_hex, p_hex * p_hex, 1, 1, ndofs,
914 strides = np.array([1, q_hex * q_hex, q_hex * q_hex], dtype="int32")
915 rui_hex = ceed.StridedElemRestriction(nelem_hex, q_hex * q_hex, 1,
920 bu_hex = ceed.BasisTensorH1Lagrange(dim, 1, p_hex, q_hex, libceed.GAUSS)
923 qf_setup_hex = ceed.QFunction(1, qfs.setup_mass_2d,
925 qf_setup_hex.add_input("weights", 1, libceed.EVAL_WEIGHT)
927 qf_setup_hex.add_output("rho", 1, libceed.EVAL_NONE)
929 qf_mass_hex = ceed.QFunction(1, qfs.apply_mass,
931 qf_mass_hex.add_input("rho", 1, libceed.EVAL_NONE)
932 qf_mass_hex.add_input("u", 1, libceed.EVAL_INTERP)
933 qf_mass_hex.add_output("v", 1, libceed.EVAL_INTERP)
981 nx_tet, ny_tet, nx_hex = 3, 1, 3
982 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
988 for i in range(ny * 2 + 1):
989 for j in range(nx * 2 + 1):
990 x_array[i + j * (ny * 2 + 1)] = i / (2 * ny)
991 x_array[i + j * (ny * 2 + 1) + ndofs] = j / (2 * nx)
1006 offset = col * 2 + row * (nx * 2 + 1) * 2
1009 indx_tet[i * 2 * p_tet + 1] = 9 + offset
1011 indx_tet[i * 2 * p_tet + 3] = 1 + offset
1025 ru_tet = ceed.ElemRestriction(nelem_tet, p_tet, 1, 1, ndofs, indx_tet,
1027 strides = np.array([1, q_tet, q_tet], dtype="int32")
1028 rui_tet = ceed.StridedElemRestriction(nelem_tet, q_tet, 1, nqpts_tet,
1039 bu_tet = ceed.BasisH1(libceed.TRIANGLE, 1, p_tet, q_hex, interp, grad, qref,
1046 qf_setup_tet = ceed.QFunction(1, qfs.setup_mass_2d,
1048 qf_setup_tet.add_input("weights", 1, libceed.EVAL_WEIGHT)
1050 qf_setup_tet.add_output("rho", 1, libceed.EVAL_NONE)
1052 qf_mass_tet = ceed.QFunction(1, qfs.apply_mass,
1054 qf_mass_tet.add_input("rho", 1, libceed.EVAL_NONE)
1055 qf_mass_tet.add_input("u", 1, libceed.EVAL_INTERP)
1056 qf_mass_tet.add_output("v", 1, libceed.EVAL_INTERP)
1078 offset = (nx_tet * 2 + 1) * (ny_tet * 2) * (1 + row) + col * 2
1083 k * (nx_hex * 2 + 1) + j
1088 ru_hex = ceed.ElemRestriction(nelem_hex, p_hex * p_hex, 1, 1, ndofs,
1090 strides = np.array([1, q_hex * q_hex, q_hex * q_hex], dtype="int32")
1091 rui_hex = ceed.StridedElemRestriction(nelem_hex, q_hex * q_hex, 1,
1096 bu_hex = ceed.BasisTensorH1Lagrange(dim, 1, p_hex, q_hex, libceed.GAUSS)
1099 qf_setup_hex = ceed.QFunction(1, qfs.setup_mass_2d,
1101 qf_setup_hex.add_input("weights", 1, libceed.EVAL_WEIGHT)
1103 qf_setup_hex.add_output("rho", 1, libceed.EVAL_NONE)
1105 qf_mass_hex = ceed.QFunction(1, qfs.apply_mass,
1107 qf_mass_hex.add_input("rho", 1, libceed.EVAL_NONE)
1108 qf_mass_hex.add_input("u", 1, libceed.EVAL_INTERP)
1109 qf_mass_hex.add_output("v", 1, libceed.EVAL_INTERP)
1136 u.set_value(1.)
1158 nx_tet, ny_tet, nx_hex = 3, 1, 3
1159 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
1173 offset = col * 2 + row * (nx * 2 + 1) * 2
1176 indx_tet[i * 2 * p_tet + 1] = 9 + offset
1178 indx_tet[i * 2 * p_tet + 3] = 1 + offset
1192 ru_tet = ceed.ElemRestriction(nelem_tet, p_tet, 1, 1, ndofs, indx_tet,
1194 strides = np.array([1, q_tet, q_tet], dtype="int32")
1195 rui_tet = ceed.StridedElemRestriction(nelem_tet, q_tet, 1, nqpts_tet,
1206 bu_tet = ceed.BasisH1(libceed.TRIANGLE, 1, p_tet, q_hex, interp, grad, qref,
1213 qf_setup_tet = ceed.QFunction(1, qfs.setup_mass_2d,
1215 qf_setup_tet.add_input("weights", 1, libceed.EVAL_WEIGHT)
1217 qf_setup_tet.add_output("rho", 1, libceed.EVAL_NONE)
1219 qf_mass_tet = ceed.QFunction(1, qfs.apply_mass,
1221 qf_mass_tet.add_input("rho", 1, libceed.EVAL_NONE)
1222 qf_mass_tet.add_input("u", 1, libceed.EVAL_INTERP)
1223 qf_mass_tet.add_output("v", 1, libceed.EVAL_INTERP)
1247 offset = (nx_tet * 2 + 1) * (ny_tet * 2) * (1 + row) + col * 2
1252 k * (nx_hex * 2 + 1) + j
1258 ru_hex = ceed.ElemRestriction(nelem_hex, p_hex * p_hex, 1, 1, ndofs,
1260 strides = np.array([1, q_hex * q_hex, q_hex * q_hex], dtype="int32")
1261 rui_hex = ceed.StridedElemRestriction(nelem_hex, q_hex * q_hex, 1,
1266 bu_hex = ceed.BasisTensorH1Lagrange(dim, 1, p_hex, q_hex, libceed.GAUSS)
1269 qf_setup_hex = ceed.QFunction(1, qfs.setup_mass_2d,
1271 qf_setup_hex.add_input("weights", 1, libceed.EVAL_WEIGHT)
1273 qf_setup_hex.add_output("rho", 1, libceed.EVAL_NONE)
1275 qf_mass_hex = ceed.QFunction(1, qfs.apply_mass,
1277 qf_mass_hex.add_input("rho", 1, libceed.EVAL_NONE)
1278 qf_mass_hex.add_input("u", 1, libceed.EVAL_INTERP)
1279 qf_mass_hex.add_output("v", 1, libceed.EVAL_INTERP)
1330 nx_tet, ny_tet, nx_hex = 3, 1, 3
1331 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
1337 for i in range(ny * 2 + 1):
1338 for j in range(nx * 2 + 1):
1339 x_array[i + j * (ny * 2 + 1)] = i / (2 * ny)
1340 x_array[i + j * (ny * 2 + 1) + ndofs] = j / (2 * nx)
1355 offset = col * 2 + row * (nx * 2 + 1) * 2
1358 indx_tet[i * 2 * p_tet + 1] = 9 + offset
1360 indx_tet[i * 2 * p_tet + 3] = 1 + offset
1374 ru_tet = ceed.ElemRestriction(nelem_tet, p_tet, 1, 1, ndofs, indx_tet,
1376 strides = np.array([1, q_tet, q_tet], dtype="int32")
1377 rui_tet = ceed.StridedElemRestriction(nelem_tet, q_tet, 1, nqpts_tet,
1388 bu_tet = ceed.BasisH1(libceed.TRIANGLE, 1, p_tet, q_hex, interp, grad, qref,
1395 qf_setup_tet = ceed.QFunction(1, qfs.setup_mass_2d,
1397 qf_setup_tet.add_input("weights", 1, libceed.EVAL_WEIGHT)
1399 qf_setup_tet.add_output("rho", 1, libceed.EVAL_NONE)
1401 qf_mass_tet = ceed.QFunction(1, qfs.apply_mass,
1403 qf_mass_tet.add_input("rho", 1, libceed.EVAL_NONE)
1404 qf_mass_tet.add_input("u", 1, libceed.EVAL_INTERP)
1405 qf_mass_tet.add_output("v", 1, libceed.EVAL_INTERP)
1427 offset = (nx_tet * 2 + 1) * (ny_tet * 2) * (1 + row) + col * 2
1432 k * (nx_hex * 2 + 1) + j
1438 ru_hex = ceed.ElemRestriction(nelem_hex, p_hex * p_hex, 1, 1, ndofs,
1440 strides = np.array([1, q_hex * q_hex, q_hex * q_hex], dtype="int32")
1441 rui_hex = ceed.StridedElemRestriction(nelem_hex, q_hex * q_hex, 1,
1446 bu_hex = ceed.BasisTensorH1Lagrange(dim, 1, p_hex, q_hex, libceed.GAUSS)
1449 qf_setup_hex = ceed.QFunction(1, qfs.setup_mass_2d,
1451 qf_setup_hex.add_input("weights", 1, libceed.EVAL_WEIGHT)
1453 qf_setup_hex.add_output("rho", 1, libceed.EVAL_NONE)
1455 qf_mass_hex = ceed.QFunction(1, qfs.apply_mass,
1457 qf_mass_hex.add_input("rho", 1, libceed.EVAL_NONE)
1458 qf_mass_hex.add_input("u", 1, libceed.EVAL_INTERP)
1459 qf_mass_hex.add_output("v", 1, libceed.EVAL_INTERP)
1486 u.set_value(1.)
1497 v.set_value(1.)
1521 ndofs = (nx * 2 + 1) * (ny * 2 + 1)
1527 for i in range(nx * 2 + 1):
1528 for j in range(ny * 2 + 1):
1529 x_array[i + j * (nx * 2 + 1) + 0 * ndofs] = i / (2 * nx)
1530 x_array[i + j * (nx * 2 + 1) + 1 * ndofs] = j / (2 * ny)
1542 offset = col * (p - 1) + row * (nx * 2 + 1) * (p - 1)
1545 indx[p * (p * i + k) + j] = offset + k * (nx * 2 + 1) + j
1549 ru = ceed.ElemRestriction(nelem, p * p, 1, 1, ndofs, indx,
1551 strides = np.array([1, q * q, q * q], dtype="int32")
1552 rui = ceed.StridedElemRestriction(nelem, q * q, 1, nqpts, strides)
1556 bu = ceed.BasisTensorH1Lagrange(dim, 1, p, q, libceed.GAUSS)
1575 nx = nelem + 1
1576 nu_coarse = nelem * (p_coarse - 1) + 1
1577 nu_fine = nelem * (p_fine - 1) + 1
1596 indx[2 * i + 1] = i + 1
1597 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
1603 indu_coarse[p_coarse * i + j] = i * (p_coarse - 1) + j
1611 indu_fine[p_fine * i + j] = i * (p_fine - 1) + j
1615 strides = np.array([1, q, q], dtype="int32")
1617 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
1620 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
1621 bu_coarse = ceed.BasisTensorH1Lagrange(1, ncomp, p_coarse, q, libceed.GAUSS)
1622 bu_fine = ceed.BasisTensorH1Lagrange(1, ncomp, p_fine, q, libceed.GAUSS)
1628 qf_setup = ceed.QFunction(1, qfs.setup_mass,
1630 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
1631 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
1632 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
1634 qf_mass = ceed.QFunction(1, qfs.apply_mass_two,
1636 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
1711 nx = nelem + 1
1712 nu_coarse = nelem * (p_coarse - 1) + 1
1713 nu_fine = nelem * (p_fine - 1) + 1
1732 indx[2 * i + 1] = i + 1
1733 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
1739 indu_coarse[p_coarse * i + j] = i * (p_coarse - 1) + j
1747 indu_fine[p_fine * i + j] = i * (p_fine - 1) + j
1751 strides = np.array([1, q, q], dtype="int32")
1753 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
1756 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
1757 bu_coarse = ceed.BasisTensorH1Lagrange(1, ncomp, p_coarse, q, libceed.GAUSS)
1758 bu_fine = ceed.BasisTensorH1Lagrange(1, ncomp, p_fine, q, libceed.GAUSS)
1764 qf_setup = ceed.QFunction(1, qfs.setup_mass,
1766 qf_setup.add_input("weights", 1, libceed.EVAL_WEIGHT)
1767 qf_setup.add_input("dx", 1, libceed.EVAL_GRAD)
1768 qf_setup.add_output("rho", 1, libceed.EVAL_NONE)
1770 qf_mass = ceed.QFunction(1, qfs.apply_mass_two,
1772 qf_mass.add_input("rho", 1, libceed.EVAL_NONE)
1796 1, ncomp, p_coarse, p_fine, libceed.GAUSS_LOBATTO)
1848 ncomp = 1
1849 nx = nelem + 1
1850 nu_coarse = nelem * (p_coarse - 1) + 1
1851 nu_fine = nelem * (p_fine - 1) + 1
1870 indx[2 * i + 1] = i + 1
1871 rx = ceed.ElemRestriction(nelem, 2, 1, 1, nx, indx,
1877 indu_coarse[p_coarse * i + j] = i * (p_coarse - 1) + j
1885 indu_fine[p_fine * i + j] = i * (p_fine - 1) + j
1889 strides = np.array([1, q, q], dtype="int32")
1891 rui = ceed.StridedElemRestriction(nelem, q, 1, q * nelem, strides)
1894 bx = ceed.BasisTensorH1Lagrange(1, 1, 2, q, libceed.GAUSS)
1895 bu_coarse = ceed.BasisTensorH1Lagrange(1, ncomp, p_coarse, q, libceed.GAUSS)
1896 bu_fine = ceed.BasisTensorH1Lagrange(1, ncomp, p_fine, q, libceed.GAUSS)
1925 1, ncomp, p_coarse, p_fine, libceed.GAUSS_LOBATTO)