-
Notifications
You must be signed in to change notification settings - Fork 19
Second Coupling Integral (times d tilde/dn) #1589
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
base: devel
Are you sure you want to change the base?
Second Coupling Integral (times d tilde/dn) #1589
Conversation
|
It looks like you're making updates unrelated to the coupling integral in this PR, unless I've missed something. Would you mind updating this PR to include just the coupling integral and its test? I'm unsure where the nonconforming fixture is used—please correct me if I'm wrong. Thanks! |
| typename IndexType, typename IntersectionFieldViewType, | ||
| typename IntersectionFactor, typename IntersectionJacobianType, | ||
| typename TransferFunctionDerivativeType, typename CallableType> | ||
| KOKKOS_FUNCTION void coupling_integral_dnshape( |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Precomputing edge-related values and storing them in the assembly can clean up this code, improve readability, and facilitate integration into coupling_integral.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is it alright if we do this in a future PR? I was waiting for this one to go through before updating my work on #1342 in order to put the first coupling_integral in the same test executable.
Precomputation and storage would require changes to nonconforming_interfaces that will bloat this PR even further.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
That sounds good. Can you create another issue for that refactor and link it here?
Also, see my other comment.
lsawade
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I'm not quite understanding why you are rewriting GLL routines here.
See: core/specfem/quadrature/gll/lagrange_poly.cpp
I had not seen that, as I overlooked the |
|
@int-ptr-ptr, how are your notes coming along? |
Description
Please describe the changes/features in this pull request.
Issue Number
If there is an issue created for these changes, link it here
solves #1550
Checklist
Please make sure to check developer documentation on specfem docs.