@@ -204,8 +204,7 @@ def time_update(Sol,flux,mpv,t,tout,ud,elem,node,steps,th,bld=None,writer=None,d
204204 up = Sol_half_new .rhou / Sol_half_new .rho
205205 vp = Sol_half_new .rhov / Sol_half_new .rho
206206 Yp = Sol_half_new .rhoY / Sol_half_new .rho - mpv .HydroState .Y0 .reshape (1 ,- 1 )
207- # vp = 0.0
208- # Yp = 0.0
207+
209208 pi = mpv_half_new .p2_nodes
210209
211210 boundary .rayleigh_damping (Sol , mpv , ud , [up , vp , Yp , pi , t + 0.5 * dt ])
@@ -218,10 +217,6 @@ def time_update(Sol,flux,mpv,t,tout,ud,elem,node,steps,th,bld=None,writer=None,d
218217 ud .rf_bot .eigenfunction ((t + 0.5 * dt ), 1 , grid = 'n' )
219218 _ , _ , _ , pi_n = ud .rf_bot .dehatter (th , grid = 'n' )
220219
221- # kernel = np.ones([2]*mpv.p2_nodes.ndim)
222- # pi_n = np.zeros_like(mpv.p2_nodes)
223- # pi_n[node.i1] = signal.fftconvolve(pi, kernel, mode='valid') / kernel.sum()
224-
225220 boundary .rayleigh_damping (Sol , mpv , ud , [up , vp , Yp , pi_n , t + 0.5 * dt ])
226221 boundary .set_explicit_boundary_data (Sol , elem , ud , th , mpv )
227222
0 commit comments