Skip to content

autoregressive covariance operator: fix scaling for dim>1 - #5358

Open
JHopeCollins wants to merge 4 commits into
mainfrom
JHopeCollins/covariance-fixes
Open

autoregressive covariance operator: fix scaling for dim>1#5358
JHopeCollins wants to merge 4 commits into
mainfrom
JHopeCollins/covariance-fixes

Conversation

@JHopeCollins

Copy link
Copy Markdown
Member
  1. Corrected diffusion coefficient for autoregressive covariance operator in multiple spatial dimensions.
  2. Fixed the IPDG coefficient for higher degree elements.

@JHopeCollins JHopeCollins self-assigned this Aug 18, 2026

@dham dham left a comment

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Other than the warning, this is good.

L = lengthscale_m(Lar, m, dim)
if dim == 1:
num = (2**(2*m - 1))*factorial(m - 1)**2
den = factorial(2*m - 2)

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

should warn if den is zero.

return 4*math.pi*(m - 1)*(L**2)
elif dim == 3:
num = (2**(2*m - 1))*math.pi*(factorial(m - 2)**2)*(m - 1)
den = factorial(2*m - 4)

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

and here.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants