search

inducer / modepy

Modes and nodes for high-order discretizations
Other
16 stars 7 forks source link

Add new JascowiecSukumar quadrature rules #77

Closed alexfikl closed 5 months ago

alexfikl commented 5 months ago

These are in the same boat as XiaoGimbutasSimplexQuadrature for tets with mostly similar node counts but it goes all the way up to order 20.

>>> [rule["points"].shape[0] for rule in xg_table.values()]
[1, 4, 6, 11, 14, 23, 31, 44, 57, 74, 95, 122, 146, 177, 214]
>>> [rule["points"].T.shape[0] for rule in js_table.values()]
[1, 4, 8, 14, 14, 24, 35, 46, 59, 81, 110, 168, 172, 204, 264, 304, 364, 436, 487, 552]

They're taken from this paper https://onlinelibrary.wiley.com/doi/10.1002/nme.6528. I'm not quite sure what the copyright on these things is though.

For some context: I stumbled upon these and got excited because I remembered Stokes required a lot of oversampling and could use these.. and by the time I realized those are triangles, not tets, it was mostly implemented :(

inducer commented 5 months ago

I'm not quite sure what the copyright on these things is though.

We could simply ask them.

alexfikl commented 5 months ago

We could simply ask them.

Good point! I'll send them an email.

inducer commented 5 months ago

Thanks! And thanks for emailing them regarding the license.