-
Notifications
You must be signed in to change notification settings - Fork 98
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
convert some CTU hydro routines to C++ #731
Conversation
all CPU tests pass |
GPU works for everything except wdmerger. This is a hybrid hydro issue. We get:
|
The issue is the call to |
GPU tests now pass: CPU tests pass too (note: the failures are because of changes on development not merged into this): this is ready to be merged |
Source/hydro/Castro_ctu_hydro.cpp
Outdated
{ | ||
flatn_arr(i,j,k) = flatn_arr(i,j,k) * flatg_arr(i,j,k); | ||
|
||
if (Radiation::flatten_pp_threshold > 0.0) { |
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.
We cannot directly read class member variables in the lambdas on the device. For now we'll need to make a local copy of this variable outside the parallel for.
Source/hydro/Castro_ctu_hydro.cpp
Outdated
if ( q_arr(i-1,j,k,QU) + q_arr(i,j-1,k,QV) + q_arr(i,j,k-1,QW) > | ||
q_arr(i+1,j,k,QU) + q_arr(i,j+1,k,QV) + q_arr(i,j,k+1,QW) ) { | ||
|
||
if (q_arr(i,j,k,QPRES) < Radiation::flatten_pp_threshold * q_arr(i,j,k,QPTOT)) { |
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.
Same comment for this one :-)
PR summary
This converts the CTU consup for hydro (not radiation) , the flattening routines, and src_to_prim to C++
PR checklist
CHANGES
file has been updated