(N, x0=0, y0=0, dx=1, dy=1)
| 29 | return [(x0+i*dx/N,y0+j*dy/N) for j in range(N+1) for i in range(N+1-j)] |
| 30 | |
| 31 | def _make_quad(N, x0=0, y0=0, dx=1, dy=1): |
| 32 | return [(x0+i*dx/N,y0+j*dy/N) for j in range(N+1) for i in range(N+1-j)] + [(x0+dx-i*dx/N,1-(y0+j*dy/N)) for j in range(N+1) for i in range(N+1-j)] |
| 33 | |
| 34 | _intrules = {} |
| 35 | def get_intrules(dim:int, order: int): |