def build_c_fun_args( self, mat, state, ap, vg ):
mat, self.function = choose_scalar_or_in_el( mat, nm.float64,
terms.dw_mass,
NotImplemented )
vec = state()
bf = ap.get_base( 'v', 0, self.integral_name )
return mat, vec, 0, bf, vg, ap.econn
def build_c_fun_args( self, mat, state, ap, vg ):
# terms.dw_mass_rho_in_el is missing
mat, self.function = choose_scalar_or_in_el( mat, nm.float64,
terms.dw_mass,
NotImplemented )
ts, state0 = self.get_args( ['ts', 'parameter'], **kwargs )
dvec = state() - state0()
rhodt = mat / ts.dt
bf = ap.get_base( 'v', 0, self.integral_name )
return rhodt, dvec, 0, bf, vg, ap.econn
