OpenSees » SRC » reliability » domain » distributions

ParetoRV.cpp

Go to the documentation of this file.
00001 /* ****************************************************************** **
00002 **    OpenSees - Open System for Earthquake Engineering Simulation    **
00003 **          Pacific Earthquake Engineering Research Center            **
00004 **                                                                    **
00005 **                                                                    **
00006 ** (C) Copyright 2001, The Regents of the University of California    **
00007 ** All Rights Reserved.                                               **
00008 **                                                                    **
00009 ** Commercial use of this program without express permission of the   **
00010 ** University of California, Berkeley, is strictly prohibited.  See   **
00011 ** file 'COPYRIGHT'  in main directory for information on usage and   **
00012 ** redistribution,  and for a DISCLAIMER OF ALL WARRANTIES.           **
00013 **                                                                    **
00014 ** Developed by:                                                      **
00015 **   Frank McKenna (fmckenna@ce.berkeley.edu)                         **
00016 **   Gregory L. Fenves (fenves@ce.berkeley.edu)                       **
00017 **   Filip C. Filippou (filippou@ce.berkeley.edu)                     **
00018 **                                                                    **
00019 ** Reliability module developed by:                                   **
00020 **   Terje Haukaas (haukaas@ce.berkeley.edu)                          **
00021 **   Armen Der Kiureghian (adk@ce.berkeley.edu)                       **
00022 **                                                                    **
00023 ** ****************************************************************** */
00024                                                                         
00025 // $Revision: 1.6 $
00026 // $Date: 2003/03/04 00:44:34 $
00027 // $Source: /usr/local/cvs/OpenSees/SRC/reliability/domain/distributions/ParetoRV.cpp,v $
00028 
00029 
00030 //
00031 // Written by Terje Haukaas (haukaas@ce.berkeley.edu) 
00032 //
00033 
00034 #include <ParetoRV.h>
00035 #include <math.h>
00036 #include <string.h>
00037 #include <classTags.h>
00038 #include <OPS_Globals.h>
00039 
00040 ParetoRV::ParetoRV(int passedTag, 
00041                  double passedParameter1,
00042                  double passedParameter2,
00043                  double passedParameter3,
00044                  double passedParameter4,
00045                  double passedStartValue)
00046 :RandomVariable(passedTag, RANDOM_VARIABLE_pareto)
00047 {
00048         tag = passedTag ;
00049         k = passedParameter1;
00050         u = passedParameter2;
00051         startValue = passedStartValue;
00052 }
00053 ParetoRV::ParetoRV(int passedTag, 
00054                  double passedParameter1,
00055                  double passedParameter2,
00056                  double passedParameter3,
00057                  double passedParameter4)
00058 :RandomVariable(passedTag, RANDOM_VARIABLE_pareto)
00059 {
00060         tag = passedTag ;
00061         k = passedParameter1;
00062         u = passedParameter2;
00063         startValue = getMean();
00064 }
00065 
00066 
00067 ParetoRV::~ParetoRV()
00068 {
00069 }
00070 
00071 
00072 void
00073 ParetoRV::Print(OPS_Stream &s, int flag)
00074 {
00075 }
00076 
00077 
00078 double
00079 ParetoRV::getPDFvalue(double rvValue)
00080 {
00081         double result;
00082         if ( u <= rvValue ) {
00083                 result = k/u * pow(u/rvValue,k+1);
00084         }
00085         else {
00086                 result = 0.0;
00087         }
00088         return result;
00089 }
00090 
00091 
00092 double
00093 ParetoRV::getCDFvalue(double rvValue)
00094 {
00095         double result;
00096         if ( u <= rvValue ) {
00097                 result = 1-pow(u/rvValue,k);
00098         }
00099         else {
00100                 result = 0.0;
00101         }
00102         return result;
00103 }
00104 
00105 
00106 double
00107 ParetoRV::getInverseCDFvalue(double rvValue)
00108 {
00109         return 0.0;
00110 }
00111 
00112 
00113 const char *
00114 ParetoRV::getType()
00115 {
00116         return "PARETO";
00117 }
00118 
00119 
00120 double 
00121 ParetoRV::getMean()
00122 {
00123         return k*u/(k-1);
00124 }
00125 
00126 
00127 
00128 double 
00129 ParetoRV::getStdv()
00130 {
00131         return sqrt(k/(k-2))*(u/(k-1));
00132 }
00133 
00134 
00135 double 
00136 ParetoRV::getStartValue()
00137 {
00138         return startValue;
00139 }
00140 
00141 double ParetoRV::getParameter1()  {return k;}
00142 double ParetoRV::getParameter2()  {return u;}
00143 double ParetoRV::getParameter3()  {opserr<<"No such parameter in r.v. #"<<tag<<endln; return 0.0;}
00144 double ParetoRV::getParameter4()  {opserr<<"No such parameter in r.v. #"<<tag<<endln; return 0.0;}
00145

Generated on Mon Oct 23 15:05:26 2006 for OpenSees by doxygen 1.5.0