Function OPS_runFORMAnalysis

SRC/interpreter/OpenSeesReliabilityCommands.cpp:2918–3003

}

int OPS_runFORMAnalysis() {
    if (OPS_GetNumRemainingInputArgs() < 1) {
        opserr << "WARNING: Wrong number of input parameter to FORM "
                  "analysis\n";
        return -1;
    }

    // get file name
    const char *filename = OPS_GetString();

    ReliabilityDomain *theReliabilityDomain = cmds->getDomain();
    if (theReliabilityDomain == 0) {
        opserr << "FORMAnalysis -- ReliabilityDomain is not defined\n";
        return -1;
    }
    
    // Check for essential ingredients
    FunctionEvaluator *theFunctionEvaluator = cmds->getFunctionEvaluator();
    if (theFunctionEvaluator == 0) {
        opserr << "Need theGFunEvaluator before a FORMAnalysis can "
                  "be created\n";
        return -1;
    }

    FindDesignPointAlgorithm *theFindDesignPointAlgorithm =
        cmds->getFindDesignPointAlgorithm();
    if (theFindDesignPointAlgorithm == 0) {
        opserr << "Need theFindDesignPointAlgorithm before a "
                  "FORMAnalysis "
                  "can be created\n";
        return -1;
    }

    ProbabilityTransformation *theProbabilityTransformation =
        cmds->getProbabilityTransformation();
    if (theProbabilityTransformation == 0) {
      opserr << "FORMAnalysis - probability transformation not defined - ";
      opserr << "assuming all independent transformation" << endln;
	theProbabilityTransformation =
            new AllIndependentTransformation(theReliabilityDomain, 0);
	cmds->setProbabilityTransformation(theProbabilityTransformation);      
    }

    Domain *theStructuralDomain = cmds->getStructuralDomain();
    if (theStructuralDomain == 0) {
        opserr << "Structural Domain is not defined\n";
        return -1;
    }

    // Read input parameter(s)
    int relSensTag = 0;
    if (OPS_GetNumRemainingInputArgs() > 1) {
        const char *type = OPS_GetString();
        if (strcmp(type, "-relSens") == 0) {
            int numdata = 1;
            if (OPS_GetIntInput(&numdata, &relSensTag) < 0) {
                opserr << "ERROR: invalid input: relSensTag \n";
                return -1;