WebSVN - OpenSees - Diff - Rev 4801 and 4803 - /trunk/SRC/reliability/analysis/analysis/ImportanceSamplingAnalysis.cpp

                                                                                 ProbabilityTransformation *passedProbabilityTransformation,
                                                                                 FunctionEvaluator *passedGFunEvaluator,
                                                                                 RandomNumberGenerator *passedRandomNumberGenerator,
                                                         bool passedStartAtOrigin,
-                                                                                Tcl_Interp *passedInterp,
+                            Tcl_Interp *passedInterp,
-                                                                                int passedNumberOfSimulations,
+                                                        long int passedNumberOfSimulations,
-                                                                                double passedTargetCOV,
+                            double passedTargetCOV,
-                                                                                double passedSamplingStdv,
+                                                        double passedSamplingStdv,
-                                                                                int passedPrintFlag,
+                                                        int passedPrintFlag,
-                                                                                TCL_Char *passedFileName,
+                                                        TCL_Char *passedFileName,
-                                                                                int passedAnalysisTypeTag)
+                                                        int passedAnalysisTypeTag)
 int
 ImportanceSamplingAnalysis::analyze(void)
+{
         // Alert the user that the simulation analysis has started
-        opserr << "Sampling Analysis is running ... " << endln;
+        opserr << "ImportanceSampling Analysis is running ... " << endln;
         // Declaration of some of the data used in the algorithm
         double gFunctionValue;
-        int result, I, k = 1, seed = 1;
+        int result, I, seed = 1;
+    long int k = 1;
         double det_covariance, phi, h, q;
         int numRV = theReliabilityDomain->getNumberOfRandomVariables();
         int numLsf = theReliabilityDomain->getNumberOfLimitStateFunctions();
         Vector x(numRV);
         Vector z(numRV);
                         inputFile.close();
+                }
+        }
-        // Transform start point into standard normal space,
-        // unless it is the origin that is to be sampled around
+    // get starting x values from parameter directly
-        Vector startPointY(numRV);
-        if (startAtOrigin)
-          startPointY.Zero();
-        else {
-          //theReliabilityDomain->getStartPoint(x);
-          for (int j = 0; j < numRV; j++) {
+    for (int j = 0; j < numRV; j++) {
-            RandomVariable *theParam = theReliabilityDomain->getRandomVariablePtrFromIndex(j);
+        RandomVariable *theParam = theReliabilityDomain->getRandomVariablePtrFromIndex(j);
-            x(j) = theParam->getStartValue();
+        double rvVal = theParam->getStartValue();
+        x(j) = rvVal;
+        // now we should update the parameter value
+        // KRM should this be a parameter reference now?
+        theParam->setCurrentValue(rvVal);
+    }
+        // Transform start point into standard normal space
+        Vector startPointY(numRV);
-          result = theProbabilityTransformation->transform_x_to_u(startPointY);
+        result = theProbabilityTransformation->transform_x_to_u(startPointY);
-          if (result < 0) {
+        if (result < 0) {
             opserr << "ImportanceSamplingAnalysis::analyze() - could not " << endln
                    << " transform x to u. " << endln;
             return -1;
-          }
+    }
         // Initial declarations
         Vector cov_of_q_bar(numLsf);
         Vector q_bar(numLsf);
         bool isFirstSimulation = true;
         while( (k<=numberOfSimulations && govCov>targetCOV || k<=2) ) {
                 // Keep the user posted
                 if (printFlag == 1 || printFlag == 2) {
+            sprintf(myString,"%li",k);
-                        opserr << "Sample #" << k << ":" << endln;
+                        opserr << "Sample #" << myString << ":" << endln;
+                }
                 // Create array of standard normal random numbers
                 if (isFirstSimulation) {
                 x = theProbabilityTransformation->get_x();
                 */
                 result = theProbabilityTransformation->transform_u_to_x(u, x);
                 if (result < 0) {
                   opserr << "ImportanceSamplingAnalysis::analyze() - could not transform u to x. " << endln;
+                  return -1;
+                }
+        // set values in the variable namespace
+        if (theGFunEvaluator->setVariables(x) < 0) {
+            opserr << "ImportanceSamplingAnalysis::analyze() - " << endln
+                << " could not set variables in namespace. " << endln;
-                  return -1;
+            return -1;
+        }
                 // Evaluate limit-state function
                 FEconvergence = true;
-                result = theGFunEvaluator->runAnalysis(x);
+                if (theGFunEvaluator->runAnalysis(x) < 0) {
-                if (result < 0) {
                         // In this case a failure happened during the analysis
                         // Hence, register this as failure
                         FEconvergence = false;
+                }
                         theReliabilityDomain->setTagOfActiveLimitStateFunction(lsfTag);
                         // set namespace variable for tcl procedures
                         Tcl_SetVar2Ex(interp,"RELIABILITY_lsf",NULL,Tcl_NewIntObj(lsfTag),TCL_NAMESPACE_ONLY);
-                        // Get value of limit-state function
+            // set and evaluate LSF
+            const char *lsfExpression = theLimitStateFunction->getExpression();
+            theGFunEvaluator->setExpression(lsfExpression);
-                        gFunctionValue = theGFunEvaluator->evaluateExpression();
+            gFunctionValue = theGFunEvaluator->evaluateExpression();
-                        if (!FEconvergence) {
+            if (!FEconvergence) {
                                 gFunctionValue = -1.0;
+                        }
                   int lsfTag = theLimitStateFunction->getTag();
                   int lsfIndex = theReliabilityDomain->getLimitStateFunctionIndex(lsfTag);
                         if ( q_bar(lsfIndex) == 0.0 ) {
                                 resultsOutputFile << "#######################################################################" << endln;
-                                resultsOutputFile << "#  SAMPLING ANALYSIS RESULTS, LIMIT-STATEFUNCDTION NUMBER   "
+                                resultsOutputFile << "#  SAMPLING ANALYSIS RESULTS, LIMIT-STATE FUNCTION NUMBER   "
                                         <<setiosflags(ios::left)<<setprecision(1)<<setw(4)<<lsfTag <<"      #" << endln;
                                 resultsOutputFile << "#                                                                     #" << endln;
                                 resultsOutputFile << "#  Failure did not occur, or zero response!                           #" << endln;
                                 resultsOutputFile << "#                                                                     #" << endln;
                                 resultsOutputFile << "#######################################################################" << endln << endln << endln;
         // Print summary of results to screen
-        opserr << "Simulation Analysis completed." << endln;
+        opserr << "ImportanceSampling Analysis completed." << endln;
         // Clean up
         resultsOutputFile.close();
         return 0;

Subversion Repositories OpenSees

(root)/trunk/SRC/reliability/analysis/analysis/ImportanceSamplingAnalysis.cpp @ 4803 - Rev 4801 → 4803