Function OPS_transformUtoX

SRC/interpreter/OpenSeesReliabilityCommands.cpp:2790–2859

}

int OPS_transformUtoX() {
    ProbabilityTransformation *theTransf =
        cmds->getProbabilityTransformation();
    if (theTransf == 0) {
        opserr << "ERROR: probability transformation has not been set"
               << endln;
        return -1;
    }

    ReliabilityDomain *theReliabilityDomain = cmds->getDomain();
    int nrv = theReliabilityDomain->getNumberOfRandomVariables();

    if (OPS_GetNumRemainingInputArgs() < nrv) {
        opserr << "ERROR: transformUtoX insufficient # args" << endln;
        return -1;
    }
    if (OPS_GetNumRemainingInputArgs() > nrv &&
        OPS_GetNumRemainingInputArgs() < 2 * nrv) {
        opserr << "ERROR: transformUtoX insufficient # rv tags" << endln;
        return -1;
    }

    int numData = 1;
    double val;
    Vector u(nrv);
    int loc = 0;
    // Read in u-values
    while (loc < nrv && OPS_GetNumRemainingInputArgs() > 0) {
        if (OPS_GetDoubleInput(&numData, &val) < 0) {
            OPS_ResetCurrentInputArg(-1);
            break;
        }
        u(loc) = val;
        loc++;
    }

    ID rvIndex(nrv);
    // Initialize the index to be sequential (default)
    for (int i = 0; i < nrv; i++) rvIndex(i) = i;

    int rvTag;
    loc = 0;
    // Now read in rv tags and get their indices
    while (loc < nrv && OPS_GetNumRemainingInputArgs() > 0) {
        if (OPS_GetIntInput(&numData, &rvTag) < 0) {
            OPS_ResetCurrentInputArg(-1);
            break;
        }
        rvIndex(loc) = theReliabilityDomain->getRandomVariableIndex(rvTag);
        loc++;
    }

    // Map in
    Vector uSorted(nrv);
    for (int i = 0; i < nrv; i++) uSorted(rvIndex(i)) = u(i);

    Vector x(nrv);
    theTransf->transform_u_to_x(uSorted, x);