Function initSectionCommands

SRC/runtime/commands/modeling/section.cpp:520–644

build the section This function assumes torsion is not NULL when num==3

// build the section
// This function assumes torsion is not NULL when num==3
static int
initSectionCommands(ClientData clientData, Tcl_Interp *interp,
                    int secTag, 
                    UniaxialMaterial *theTorsion, 
                    double Ys, double Zs, 
                    double alpha, 
                    const FiberSectionConfig& options)
{
  assert(clientData != nullptr);
  BasicModelBuilder *builder = static_cast<BasicModelBuilder*>(clientData);

  // Dimension of the structure
  int ndm = builder->getNDM();

  SectionBuilder  *sbuilder = nullptr;
  FrameSection    *section  = nullptr;
  // Create 2d section
  if (ndm == 2) {
    if (options.isND) {
      if (options.isNew) {
        auto sec = new FrameSolidSection3d(secTag, 30);
        sbuilder = new FiberSectionBuilder<2, NDMaterial, FrameSolidSection3d>(*builder, *sec);
        section = sec;
      }
      else if (options.isWarping) {
        auto sec = new NDFiberSectionWarping2d(secTag, 30, alpha);
        sbuilder = new FiberSectionBuilder<2, NDMaterial, NDFiberSectionWarping2d>(*builder, *sec);
        section = sec;
      }
      else {
        auto sec = new NDFiberSection2d(secTag, options.computeCentroid);
        sbuilder = new FiberSectionBuilder<2, NDMaterial, NDFiberSection2d>(*builder, *sec);
        section = sec;
      }
    } else {
      if (options.isThermal) {
        auto sec = new FiberSection2dThermal(secTag, options.computeCentroid);
        sbuilder = new FiberSectionBuilder<2, UniaxialMaterial, FiberSection2dThermal>(*builder, *sec);
        section = sec;
      } else {
        auto sec = new FiberSection2d(secTag, options.computeCentroid);
        sbuilder = new FiberSectionBuilder<2, UniaxialMaterial, FiberSection2d>(*builder, *sec);
        section = sec;
      }
    }
  }

  else if (ndm == 3) {

    if (options.isND) {
      if (options.isNew) {
        auto sec = new FrameSolidSection3d(secTag, 30);
        sbuilder = new FiberSectionBuilder<3, NDMaterial, FrameSolidSection3d>(*builder, *sec);
        section = sec;
      }
      else {
        auto sec = new NDFiberSection3d(secTag,
                                        options.computeCentroid);