Entropy filter - point set combination

Cases
,
1

Hello,

I was trying to use the entropy filtering but I don’t understand why this point set combination is not valid


[solver-interfaces-quad]
flux-pts = gauss-legendre-lobatto

[solver-interfaces-tri]
flux-pts = alpha-opt

[solver-elements-pri]
flux-pts = alpha-opt
soln-pts = alpha-opt~gauss-legendre-lobatto

and it causes the error:

ValueError: Entropy filter requires flux points to be a subset of solution points or a convex combination thereof.

Any insight on that?

Best

2

Hi,

I think this should work. It is possible that there are still some tolerance issues with how the points are tabulated for the tri and quad vs the prism. What order does it fail at? And is the mesh a pure prism mesh?

Also, FYI, I don’t think specifying flux points for the elements is required.

Peter

3

Can you confirm what version of PyFR you are using?

4

Hi, thanks for the answer.

And is the mesh a pure prism mesh?

Yes I confirm that the mesh is pure prism.

Can you confirm what version of PyFR you are using?

I’m using the latest 2.0.2

And is the mesh a pure prism mesh?

It works with order 1, 2, 3 but fails with order 4

Best

5

Can you run the following script and let me know if it gives any errors:

import numpy as np

from pyfr.inifile import Inifile
from pyfr.shapes import PriShape

cfg = """
[solver]
order = 4

[solver-interfaces-quad]
flux-pts = gauss-legendre-lobatto

[solver-interfaces-tri]
flux-pts = alpha-opt

[solver-elements-pri]
soln-pts = alpha-opt~gauss-legendre-lobatto
"""
cfg = Inifile(cfg)

p = PriShape(None, cfg)
m0 = p.m0
mrowsum = np.max(np.abs(np.sum(m0, axis=1) - 1.0))
if np.min(m0) < -1e-8 or mrowsum > 1e-8:
    raise ValueError('Entropy filter requires flux points to be a '
                     'subset of solution points or a convex '
                     'combination thereof.')

Regards, Freddie.

6

I get the same error

7

That script works fine for me under a fresh install of 2.0.2.

Regards, Freddie.

8

That script also works for me with latest release. Could you double check the version you are running?