Function OPS_sectionLocation

SRC/material/section/SectionForceDeformation.cpp:103–181

}

int OPS_sectionLocation()
{
    // make sure at least one other argument to contain type of system
    if (OPS_GetNumRemainingInputArgs() < 1) {
	opserr << "WARNING want - sectionLocation eleTag? <secNum?> \n";
	return -1;
    }

    //opserr << "sectionLocation: ";
    //for (int i = 0; i < argc; i++)
    //  opserr << argv[i] << ' ' ;
    //opserr << endln;

    int numdata = 1;
    int tag;
    if (OPS_GetIntInput(&numdata, &tag) < 0) {
	opserr << "WARNING sectionLocation eleTag? <secNum?> - could not read int input? \n";
	return -1;
    }

    int secNum = 0;
    if (OPS_GetNumRemainingInputArgs() > 0) {
      if (OPS_GetIntInput(&numdata, &secNum) < 0) {
	opserr << "WARNING sectionLocation eleTag? <secNum?> - could not read int input? \n";
	return -1;
      }
    }
    
    Domain* theDomain = OPS_GetDomain();
    if (theDomain == 0) return -1;

    Element *theElement = theDomain->getElement(tag);
    if (theElement == 0) {
	opserr << "WARNING sectionLocation element with tag " << tag << " not found in domain \n";
	return -1;
    }

    int argcc = 1;
    char a[80] = "integrationPoints";
    const char *argvv[1];
    argvv[0] = a;

    DummyStream dummy;

    Response *theResponse = theElement->setResponse(argvv, argcc, dummy);
    if (theResponse == 0) {
	return 0;
    }

    theResponse->getResponse();
    Information &info = theResponse->getInformation();

    const Vector &theVec = *(info.theVector);
    int Np = theVec.Size();

    if (secNum > 0 && secNum <= Np) { // One IP
      double value = theVec(secNum-1);
      numdata = 1;