-
-
Couldn't load subscription status.
- Fork 72
[WIP] (issue #237) Allow distinct type for grid vs. function value in BC constructor #260
base: master
Are you sure you want to change the base?
Conversation
Test Robin BC, Dirichlet, Neumann with complex function values
|
|
||
| denoml = αl[2] .+ αl[3:end] ⋅ s0_l | ||
| denomr = αr[2] .+ αr[3:end] ⋅ s0_r | ||
| denoml = αl[2] .+ αl[3:end]' ⋅ s0_l |
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.
why adjoint?
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.
⋅(a,b) complex conjugates a. I believe this should not be complex conjugated.
Should probably have been conj. instead of adjoint (same result due to the dot).
| #RobinBC(al::T, bl::T, cl::T, dx_l::T, ar::T, br::T, cr::T, dx_r::T, order = 1) where T = RobinBC([al,bl, cl], [ar, br, cr], dx_l, order) | ||
|
|
||
| Base.:*(Q::AffineBC, u::AbstractVector) = BoundaryPaddedVector(Q.a_l ⋅ u[1:length(Q.a_l)] + Q.b_l, Q.a_r ⋅ u[(end-length(Q.a_r)+1):end] + Q.b_r, u) | ||
| Base.:*(Q::AffineBC, u::AbstractVector) = BoundaryPaddedVector( Q.a_l'⋅ u[1:length(Q.a_l)] + Q.b_l, Q.a_r' ⋅ u[(end-length(Q.a_r)+1):end] + Q.b_r, u ) |
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.
The idea that the dot-product should not complex conjugate came from this line. I figured it should be the same for GeneralBC
Hi there,
This should allow RobinBC and GeneralBC to take different types between function values and grid spacing (#237).
Suggestions welcome! I'm still new to the Github workflow, so please be gentle ;)