Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Mesh check in make_same_mesh_connection too strict? #358

Open
majosm opened this issue Nov 1, 2022 · 2 comments · May be fixed by #359
Open

Mesh check in make_same_mesh_connection too strict? #358

majosm opened this issue Nov 1, 2022 · 2 comments · May be fixed by #359
Assignees

Comments

@majosm
Copy link
Collaborator

majosm commented Nov 1, 2022

@anderson2981 and @tulioricci have both run into errors in make_same_mesh_connection when restarting with different orders. Specifically, this check is failing:

    if from_discr.mesh is not to_discr.mesh:
        raise ValueError("from_discr and to_discr must be based on "
                "the same mesh")

I suspect it's related to the mesh boundary tag relabeling that grudge does when constructing DiscretizationCollections. Strangely, though, changing the check condition to

   if from_discr.mesh.groups is not to_discr.mesh.groups:
       ...

doesn't appear to help, so maybe something else is going on? (Edit: I don't think this would be correct anyway, since it wouldn't detect transformed meshes.)

This can be worked around for now by just commenting out the check, but I should take a closer look.

@majosm majosm self-assigned this Nov 1, 2022
@alexfikl
Copy link
Collaborator

alexfikl commented Nov 1, 2022

I had something in this direction here if it helps: alexfikl@6b35e06

It mostly just checks that the structure of the discretizations is the same, but it doesn't have any special handling for transformed meshes. If you think it would be useful, I can open a PR.

@majosm majosm linked a pull request Nov 1, 2022 that will close this issue
@majosm
Copy link
Collaborator Author

majosm commented Nov 1, 2022

I had something in this direction here if it helps: alexfikl@6b35e06

It mostly just checks that the structure of the discretizations is the same, but it doesn't have any special handling for transformed meshes. If you think it would be useful, I can open a PR.

#359 might be good enough for this particular issue; we'll see.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging a pull request may close this issue.

2 participants