Method convectParticle

SRC/element/PFEMElement/BackgroundMesh.cpp:3531–3641

}

int BackgroundMesh::convectParticle(Particle* pt, VInt index,
                                    int nums) {
    Domain* domain = OPS_GetDomain();
    if (domain == 0) return 0;

    int ndm = OPS_GetNDM();

    // check dt
    double dt = pt->getDt();
    if (dt <= 0) {
        return 0;
    }

    // get corners
    VVInt indices, ndtags;
    VVDouble crds;
    std::vector<BackgroundType> types;
    VVDouble vels, dvns;
    VDouble pns, dpns;

    // convect in a cell
    while (pt->getDt() > 0) {
        // get subdt
        double subdt = dt / nums;
        if (pt->getDt() < subdt) {
            subdt = pt->getDt();
        }

        // particle crds
        const VDouble& pcrds = pt->getCrds();

        // new index
        VInt newIndex;
        lowerIndex(pcrds, newIndex);

        // update corners
        if (indices.empty() || newIndex != index) {
            index = newIndex;

            VVInt temp;
            getCorners(index, 1, temp);
            indices = temp;
            indices[2] = temp[3];
            indices[3] = temp[2];
            if (ndm == 3) {
                indices[6] = temp[7];
                indices[7] = temp[6];
            }

            // get corner coordinates, types, and velocities
            // pressures for structural damping
            ndtags.assign(indices.size(), VInt());
            crds.assign(indices.size(), VDouble());
            types.assign(indices.size(), BACKGROUND_FIXED);
            vels.assign(indices.size(), VDouble());
            dvns.assign(indices.size(), VDouble());
            pns.assign(indices.size(), 0.0);
            dpns.assign(indices.size(), 0.0);