Method incrTrialDisp

SRC/domain/node/Node.cpp:831–868

}

int
Node::incrTrialDisp(const Vector &incrDispl)
{
    // check vector arg is of correct size
    if (incrDispl.Size() != numberDOF) {
	opserr << "WARNING Node::incrTrialDisp() - incompatible sizes\n";
	return -2;
    }    

    if (rotation != nullptr && this->getNumberDOF() >= 6)
      rotation[1] = rotation[1]*Versor::from_vector(&disp[3*numberDOF+3]);


    // create a copy if no trial exists and add committed
    if (trialDisp == 0) {
	if (this->createDisp() < 0) {
	    opserr << "FATAL Node::incrTrialDisp() - ran out of memory\n";
	    exit(-1);
	}    
	for (int i = 0; i<numberDOF; i++) {
	  double incrDispI = incrDispl(i);
	  disp[i] = incrDispI;
	  disp[i+2*numberDOF] = incrDispI;
	  disp[i+3*numberDOF] = incrDispI;
	}
	return 0;
    }

    // otherwise set trial = incr + trial
    for (int i = 0; i<numberDOF; i++) {
	  double incrDispI = incrDispl(i);
	  disp[i] += incrDispI;
	  disp[i+2*numberDOF] += incrDispI;
	  disp[i+3*numberDOF] = incrDispI;
    }

    return 0;
}


int