(L, bd_indices, constraints, is_variable)
| 66 | return bd_vertex_indices, bd_vertex_positions; |
| 67 | |
| 68 | def init_harmonic_system(L, bd_indices, constraints, is_variable): |
| 69 | rhs = -L[:,bd_indices].dot(constraints); |
| 70 | |
| 71 | M = L[:, is_variable][is_variable, :]; |
| 72 | rhs = rhs[is_variable, :]; |
| 73 | return M, rhs; |
| 74 | |
| 75 | @pymesh.timethis |
| 76 | def solve(M, rhs, solver_type, tol): |