In seed.c there is the assignment:
cr->p[6*i+3+2] = SD_vinit * phy->mass[ cr->type[i] ]; cr->dpdt[3*i+2] = SD_vinit;
If SD_vinit is a velocity then the second line assigns it to a force!
Log in to post a comment.