Method setResponse

SRC/material/section/FiberSectionWarping3d.cpp:1196–1398

}

Response*
FiberSectionWarping3d::setResponse(const char **argv, int argc, OPS_Stream &output)
{
  const ID &type = this->getType();
  int typeSize = this->getOrder();

  Response *theResponse = 0;
  
  static double yLocs[10000];
  static double zLocs[10000];
  static double fiberArea[10000];
  static double omega[10000];
  
  if (sectionIntegr != 0) {
    sectionIntegr->getFiberLocations(numFibers, yLocs, zLocs);
    sectionIntegr->getFiberWeights(numFibers, fiberArea);
    sectionIntegr->getFiberSectorials(numFibers, omega);      
  }  
  else {
    for (int i = 0; i < numFibers; i++) {
      yLocs[i] = matData[4*i];
      zLocs[i] = matData[4*i+1];
      fiberArea[i] = matData[4*i+2];
      omega[i] = matData[4*i+3];	
    }
  }
    
  if (argc > 2 && strcmp(argv[0],"fiber") == 0) {

    int key = numFibers;
    int passarg = 2;
    
      if (argc <= 3)	{  // fiber number was input directly
	
	key = atoi(argv[1]);
	
      } else if (argc > 4) {         // find fiber closest to coord. with mat tag
	int matTag = atoi(argv[3]);
	double yCoord = atof(argv[1]);
	double zCoord = atof(argv[2]);
	double closestDist = 0.0;
	double ySearch, zSearch, dy, dz;
	double distance;
	int j;
	
	// Find first fiber with specified material tag
	for (j = 0; j < numFibers; j++) {
	  if (matTag == theMaterials[j]->getTag()) {
	    ySearch = -matData[4*j];
	    zSearch =  matData[4*j+1];
	    dy = ySearch-yCoord;
	    dz = zSearch-zCoord;
	    closestDist = sqrt(dy*dy + dz*dz);
	    key = j;
	    break;
	  }
	}