Skip to content

Commit

Permalink
Minor fix for SORM
Browse files Browse the repository at this point in the history
  • Loading branch information
dimtsap committed Mar 18, 2021
1 parent bb401b4 commit 0e0c076
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions src/UQpy/Reliability/TaylorSeries/SORM.py
Original file line number Diff line number Diff line change
Expand Up @@ -108,10 +108,10 @@ def normalize(v):
matrix_b = np.dot(np.dot(r1, hessian_g), r1.T) / np.linalg.norm(dg_u_record[-1])
kappa = np.linalg.eig(matrix_b[:self.dimension-1, :self.dimension-1])
if self.call is None:
self.Pf_sorm = [stats.norm.cdf(-self.beta_form) * np.prod(1 / (1 + self.beta_form * kappa[0]) ** 0.5)]
self.Pf_sorm = [stats.norm.cdf(-1*self.beta_form) * np.prod(1 / (1 + self.beta_form * kappa[0]) ** 0.5)]
self.beta_sorm = [-stats.norm.ppf(self.Pf_sorm)]
else:
self.Pf_sorm = self.Pf_sorm + [stats.norm.cdf(-self.beta_form) * np.prod(1 / (1 + self.beta_form *
self.Pf_sorm = self.Pf_sorm + [stats.norm.cdf(-1*self.beta_form) * np.prod(1 / (1 + self.beta_form *
kappa[0]) ** 0.5)]
self.beta_sorm = self.beta_sorm + [-stats.norm.ppf(self.Pf_sorm)]

Expand Down

0 comments on commit 0e0c076

Please sign in to comment.