Changeset 23709 for issm/trunk-jpl/src/py3/classes/qmu/@dakota_method/dakota_method.py

Timestamp:

02/11/19 04:14:06 (6 years ago)

Author:

bdef

Message:

CHG: pyhton 3 migration (not opperational yet)

File:

• issm/trunk-jpl/src/py3/classes/qmu/@dakota_method/dakota_method.py (modified) (13 diffs)

issm/trunk-jpl/src/py3/classes/qmu/@dakota_method/dakota_method.py

@@ -1 +1 @@
 #move this later
 from helpers import *
+
+from MatlabFuncs import *
 import numpy as np
 
@@ -44 +46 @@
 
         def __init__(self,*args):
-                self.method              = ''
-                self.type                        = ''
-                self.variables = []
-                self.lcspec              = []
-                self.responses = []
-                self.ghspec              = []
+                self.method   =''
+                self.type     =''
+                self.variables=[]
+                self.lcspec   =[]
+                self.responses=[]
+                self.ghspec   =[]
                 #properites
-                self.params              = struct()
+                self.params   =struct()
 
         @staticmethod
@@ -75 +77 @@
                         #given argument was a way of constructing a method
                         else:
-                                mlist=[
-                                        'dot_bfgs',
-                                        'dot_frcg',
-                                        'dot_mmfd',
-                                        'dot_slp',
-                                        'dot_sqp',
-                                        'npsol_sqp',
-                                        'conmin_frcg',
-                                        'conmin_mfd',
-                                        'optpp_cg',
-                                        'optpp_q_newton',
-                                        'optpp_fd_newton',
-                                        'optpp_newton',
-                                        'optpp_pds',
-                                        'asynch_pattern_search',
-                                        'coliny_cobyla',
-                                        'coliny_direct',
-                                        'coliny_ea',
-                                        'coliny_pattern_search',
-                                        'coliny_solis_wets',
-                                        'ncsu_direct',
-                                        'surrogate_based_local',
-                                        'surrogate_based_global',
-                                        'moga',
-                                        'soga',
-                                        'nl2sol',
-                                        'nlssol_sqp',
-                                        'optpp_g_newton',
-                                        'nond_sampling',
-                                        'nond_local_reliability',
-                                        'nond_global_reliability',
-                                        'nond_polynomial_chaos',
-                                        'nond_stoch_collocation',
-                                        'nond_evidence',
-                                        'dace',
-                                        'fsu_quasi_mc',
-                                        'fsu_cvt',
-                                        'vector_parameter_study',
-                                        'list_parameter_study',
-                                        'centered_parameter_study',
-                                        'multidim_parameter_study',
-                                        'bayes_calibration']
+                                mlist=['dot_bfgs',
+                                                         'dot_frcg',
+                                                         'dot_mmfd',
+                                                         'dot_slp',
+                                                         'dot_sqp',
+                                                         'npsol_sqp',
+                                                         'conmin_frcg',
+                                                         'conmin_mfd',
+                                                         'optpp_cg',
+                                                         'optpp_q_newton',
+                                                         'optpp_fd_newton',
+                                                         'optpp_newton',
+                                                         'optpp_pds',
+                                                         'asynch_pattern_search',
+                                                         'coliny_cobyla',
+                                                         'coliny_direct',
+                                                         'coliny_ea',
+                                                         'coliny_pattern_search',
+                                                         'coliny_solis_wets',
+                                                         'ncsu_direct',
+                                                         'surrogate_based_local',
+                                                         'surrogate_based_global',
+                                                         'moga',
+                                                         'soga',
+                                                         'nl2sol',
+                                                         'nlssol_sqp',
+                                                         'optpp_g_newton',
+                                                         'nond_sampling',
+                                                         'nond_local_reliability',
+                                                         'nond_global_reliability',
+                                                         'nond_polynomial_chaos',
+                                                         'nond_stoch_collocation',
+                                                         'nond_evidence',
+                                                         'dace',
+                                                         'fsu_quasi_mc',
+                                                         'fsu_cvt',
+                                                         'vector_parameter_study',
+                                                         'list_parameter_study',
+                                                         'centered_parameter_study',
+                                                         'multidim_parameter_study',
+                                                         'bayes_calibration']
 
                                 mlist2=[]
@@ -122 +123 @@
                                         if strncmpi(method,mlist[i],len(method)):
                                                 mlist2.append(mlist[i])
-
-                                #  check for a unique match in the list of methods
-
+                                                #  check for a unique match in the list of methods
                                 l = len(mlist2)
                                 if l == 0:
@@ -134 +133 @@
 
                                 #  assign the default values for the method
+                          # switch dm.method
                                 if dm.method in ['dot_bfgs','dot_frcg']:
                                         dm.type     ='dot'
-                                        dm.variables=['continuous_design','continuous_state']
+                                        dm.variables=['continuous_design',
+                                                                                                'continuous_state']
                                         dm.lcspec   =[]
                                         dm.responses=['objective_function']
@@ -151 +152 @@
                                 elif dm.method in ['dot_mmfd','dot_slp','dot_sqp']:
                                         dm.type     ='dot'
-                                        dm.variables=['continuous_design','continuous_state']
+                                        dm.variables=['continuous_design',
+                                                                                                'continuous_state']
                                         dm.lcspec   =['linear_inequality_constraint',
                                                                                                 'linear_equality_constraint']
@@ -236 +238 @@
                                         dm.params.max_step=1000.
                                         dm.params.gradient_tolerance=1.0e-4
-                                elif dm.method in ['optpp_q_newton','optpp_fd_newton','optpp_newton']:
+
+                                elif dm.method in ['optpp_q_newton',
+                                                                                                         'optpp_fd_newton',
+                                                                                                         'optpp_newton']:
                                         dm.type     ='optpp'
                                         dm.variables=['continuous_design',
@@ -455 +460 @@
                                         dm.params.vol_boxsize_limit=1.0e-8
 
-                                # elif dm.method in ['surrogate_based_local',
-                                #                                                                        'surrogate_based_global']:
+                                        #if dm.method in ['surrogate_based_local',
+                                        #'surrogate_based_global']:
+
                                 elif dm.method == 'moga':
                                         dm.type     ='jega'
@@ -501 +507 @@
                                         dm.params.postprocessor_type=False
                                         dm.params.orthogonal_distance=[0.01]
+
                                 elif dm.method == 'soga':
                                         dm.type     ='jega'
@@ -563 +570 @@
                                         dm.params.covariance=0
                                         dm.params.regression_stressbalances=False
+
                                 elif dm.method == 'nlssol_sqp':
                                         dm.type     ='lsq'
@@ -619 +627 @@
                                         dm.responses=['response_function']
                                         dm.ghspec   =[]
-                                        #not documented, but apparently works
+                                        #                               not documented, but apparently works
                                         dm.params.output=False
                                         dm.params.seed=False
@@ -658 +666 @@
                                         dm.responses=['response_function']
                                         dm.ghspec   =['grad']
-                                        #not documented, but may work
+                                        #                               not documented, but may work
                                         dm.params.output=False
                                         dm.params.x_gaussian_process=False
@@ -860 +868 @@
 
                                 else:
-                                        raise RuntimeError('Unimplemented method: '+str(dm.method)+'.')
+                                        raise RuntimeError('Unimplemented method: {}.'.format(dm.method))
 
                 #  if more than one argument, issue warning