Method tryFix

SRC/element/PFEMElement/BackgroundFixData.cpp:60–139

}

int
BackgroundFixData::tryFix(int ndtag, Domain& domain)
{
    // get node
    Node* node = domain.getNode(ndtag);
    if (node == 0) {
	opserr<<"WARNING: node "<<ndtag<<" not exist\n";
	return -1;
    }
    int ndf = node->getNumberDOF();

    // crds
    const Vector& crds = node->getCrds();

    // check each info
    int index = -1;
    for (int i=0; i<(int)min.size(); i++) {
	if (crds.Size() != min[i].Size()) {
	    opserr << "WARNING: ndm for the nodes and fix range are not compatible\n";
	    return -1;
	}
	if (crds.Size() != max[i].Size()) {
	    opserr << "WARNING: ndm for the nodes and fix range are not compatible\n";
	    return -1;
	}

	// if in the fix area
	bool fixed = true;
	for (int j=0; j<crds.Size(); j++) {
	    if (crds(j)<min[i](j) || crds(j)>max[i](j)) {
		fixed = false;
		break;
	    }
	}
	if (fixed) {
	    index = i;
	    break;
	}
    }

    // if not fixed
    if (index == -1) return 0;

    // set nodal states
    Vector disp = node->getTrialDisp();
    Vector vel = node->getTrialVel();
    Vector accel = node->getTrialVel();
    for (int i=0; i<ndf; i++) {
	if (i>=fix[index].Size()) break;
	if (fix[index](i) != 0) {
	    disp(i) = 0.0;
	    vel(i) = 0.0;
	    accel(i) = 0.0;
	}
    }
    node->setTrialDisp(disp);
    node->setTrialVel(vel);
    node->setTrialAccel(accel);