Method line

SRC/element/PFEMElement/ParticleGroup.cpp:414–440

}

int ParticleGroup::line(const VDouble &p1, const VDouble &p2, int num,
                        const VDouble &vel0, double p0) {
    if (num <= 0)
        return 0;

    if (p1.size() != p2.size())
        return -1;

    VDouble p1p2 = p2;
    p1p2 -= p1;
    p1p2 /= num;

    VDouble crds(p1);
    VDouble vel(crds.size());
    for (int i = 0; i < (int) vel.size(); i++) {
        if (i < (int) vel0.size()) {
            vel[i] = vel0[i];
        }
    }

    for (int i = 0; i <= num; i++) {
        this->addParticle(crds, vel, p0);
        crds += p1p2;
    }

    return 0;
}

int ParticleGroup::pointlist(VDouble &pointdata) {
    int ndm = OPS_GetNDM();