def create_high_fidelity_models(self):
rawData = mt.read_tabulated_data(self.dataPath,False,1)
self.xNorm = self.xNormalization.normalize(rawData[:,:9])
self.highCL = mt.RbfMod(self.xNorm,rawData[:,9])
self.highCD = mt.RbfMod(self.xNorm,rawData[:,10])
self.highCM = mt.RbfMod(self.xNorm,rawData[:,11])
return rawData
def create_scaling_models(self):
self.rawData = mt.read_tabulated_data(self.dataPath,False,1)
self.gammaCL = self.rawData[:,9] - self.rawData[:,12]
self.gammaCD = self.rawData[:,10] - self.rawData[:,13]
self.gammaCM = self.rawData[:,11] - self.rawData[:,14]
self.xNorm = self.xNormalization.normalize(self.rawData[:,:9])
self.scalingCL = mt.RbfMod(self.xNorm,self.gammaCL)
self.scalingCD = mt.RbfMod(self.xNorm,self.gammaCD)
self.scalingCM = mt.RbfMod(self.xNorm,self.gammaCM)
def run_test1():
data = mt.read_tabulated_data('samples.txt',False,1)
x0 = np.array([0.06343078, 0.1400427, 0.1070076, 0.13289899, -0.01762751,
-0.08041535, 0.02156679])
xl = np.array([0.04343078, 0.0900427, 0.0570076, 0.08289899, -0.06762751,
-0.13041535, -0.02843321])
xu = np.array([0.08343078, 0.1900427, 0.1570076, 0.18289899, 0.03237249,
-0.03041535, 0.07156679])
xNormalization = mt.Normalization(xl,xu)
def create_scaling_models2(self):
# after AIAA conference - 100DOE for low-fi
hifiPath = self.dataPath
lowfiPath = 'LHS9_100_CFD.txt'
rdHi = mt.read_tabulated_data(hifiPath,False,1)
rdLo = mt.read_tabulated_data(lowfiPath,False,1)
xNormLo = self.xNormalization.normalize(rdLo[:,:9])
xNormHi = self.xNormalization.normalize(rdHi[:,:9])
# --- low-fidelity model ---
self.lowCL = mt.RbfMod(xNormLo, rdLo[:,9])
self.lowCD = mt.RbfMod(xNormLo, rdLo[:,10])
self.lowCM = mt.RbfMod(xNormLo, rdLo[:,11])
# --- scaling ---
nh = len(xNormHi)
gammaCL = np.zeros(nh)
gammaCD = np.zeros(nh)
gammaCM = np.zeros(nh)
for i,xx in enumerate(xNormHi):
gammaCL[i] = rdHi[i,9] - self.lowCL(xx)
gammaCD[i] = rdHi[i,10] - self.lowCD(xx)
gammaCM[i] = rdHi[i,11] - self.lowCM(xx)
self.scalingCL = mt.RbfMod(xNormHi, gammaCL)
self.scalingCD = mt.RbfMod(xNormHi, gammaCD)
self.scalingCM = mt.RbfMod(xNormHi, gammaCM)
