#include "udf.h"
#include "dpm.h"
DEFINE_DPM_SCALAR_UPDATE(ceshi, c, t, initialize, p)
{
Injection *I;
Injection *Ilist = Get_dpm_injections();
loop(I, Ilist)
{
Particle *p_new;
p_new = new_particle(I, FALSE);
p_new->part_id = get_next_part_id();
P_INJECTION(p_new) = I;
p_new->I_id = 0;
((p_new)->stream_index) = 1;
P_POS(p_new)[0] = 0;
P_POS(p_new)[1] = 0;
P_POS(p_new)[2] = 5;
P_VEL(p_new)[0] = 0;
P_VEL(p_new)[1] = 0;
P_VEL(p_new)[2] = 0;
P_DIAM(p_new) = 0.001;
p_new->n_steps = 0;
p_new->next_time_step = solver_par.flow_time_step;
P_TIME(p_new) = solver_par.flow_time;
memcpy((char *)&(p_new->init_state),
(char *)&(p_new->state), sizeof(particle_state_t));
P_FLOW_RATE(p_new) = 0.001;
p_new->time_of_birth = 0.0;
alloc_pvars(p_new);
append_particle_to_list(p_new, &(I->p), I);/* 问题函数 */
I->n_particles++;
}
}
append_particle_to_list(p_new, &(I->p), I); 这个函数总是出问题，有人知道怎么用么这个函数
