Function ops_addPatch

SRC/interpreter/PythonModelBuilder.cpp:617–723

}

PyObject *ops_addPatch(PyObject *self, PyObject *args)
{
    // check model builder
    if(pBuilder.isSet() == false) {
	PyErr_SetString(PyExc_RuntimeError,"ERROR builder has not been set");
	return NULL;
    }
    
    OPS_ResetCommandLine(PyTuple_Size(args), 0, args);

    // check if a section is being processed
    SectionForceDeformation* section = pBuilder.getCurrentSection();
    if(section == 0) {
	PyErr_SetString(PyExc_RuntimeError,"ERROR no current section");
	return NULL;
    }

    // create fiber
    Patch *thePatch =0;
    std::string type = OPS_GetString();
    if(type == "quad" || type == "quadrilateral") {
	thePatch = (Patch*) OPS_QuadPatch();
    } else if(type == "rect" || type == "rectangular") {
	thePatch = (Patch*) OPS_RectPatch();
    } else if(type == "circ" || type == "circular") {
	thePatch = (Patch*) OPS_CircPatch();
    } else {
	PyErr_SetString(PyExc_RuntimeError,"ERROR unknown patch type");
	return NULL;
    }

    if(thePatch == 0) {
    	PyErr_SetString(PyExc_RuntimeError,"ERROR failed to create Patch");
    	return NULL;
    }

    // add fibers to the section
    int numCells = thePatch->getNumCells();
    int matTag = thePatch->getMaterialID();
    Cell** cells = thePatch->getCells();
    if(cells == 0) {
	PyErr_SetString(PyExc_RuntimeError,"ERROR out of run to create fibers");
	delete thePatch;
    	return NULL;
    }
    for(int j=0; j<numCells; j++) {
	// get fiber data
	double area = cells[j]->getArea();
	const Vector& cPos = cells[j]->getCentroidPosition();

	// create fibers
	Fiber *theFiber = 0;
	UniaxialMaterial *material = 0;
	NDMaterial *ndmaterial = 0;
	
	if(section->getClassTag() == SEC_TAG_FiberSection2d) {
	    material = OPS_getUniaxialMaterial(matTag);
	    if(material == 0) {