SciCompMod / GMGPolar

High Order Geometric Multigrid for planes in curvilinear coordinates
Apache License 2.0
14 stars 0 forks source link

Anisotropy in theta needs adaptation of define_line_splitting() #42

Open mknaranja opened 1 year ago

mknaranja commented 1 year ago

If we allow anisotropy in theta (which basically should be possible) but was never really used due to our use cases, we have to adapt:

https://github.com/mknaranja/GMGPolar/blob/main/src/smoother.cpp#L79

Further, when we calculate "delete circles" in level::define_line_splitting() (smoother.cpp line 73), we use that the grid spacing in theta is constant, dont we ? If anisotrophy is chosen in theta I think one should update that function as well!

_Originally posted by @CodingAllan in https://github.com/mknaranja/GMGPolar/pull/37#discussion_r1282896862_