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