def linearize(self): """ Calculate adjoint and read it in """ for name in _obj_names: config_in.ADJ_OBJ_FUNC = name state = adjoint(self.config_in) csvname = self.config_in.SURFACE_ADJ_FILENAME+'.csv' suffix = SU2.io.tools.get_adjointSuffix(name) csvname = SU2.io.tools.add_suffix(csvname, suffix)
def linearize(self): """ Create jacobian from adjoint results.""" self.J = None for i, name in enumerate(_obj_names): self.config_in.ADJ_OBJFUNC = name state = adjoint(self.config_in) restart2solution(self.config_in, state) csvname = self.config_in.SURFACE_ADJ_FILENAME+'.csv' col = get_sensitivities(csvname) if self.J is None: self.J = np.zeros((len(col),len(_obj_names))) self.J[:,i] = np.array(col)
def provideJ(self): """ Create jacobian from adjoint results.""" if self._first_lin: self._first_lin = False else: self.config_in.RESTART_SOL = "YES" self.J = None for i, name in enumerate(_obj_names): self.config_in.ADJ_OBJFUNC = name state = adjoint(self.config_in) restart2solution(self.config_in, state) csvname = self.config_in.SURFACE_ADJ_FILENAME + '.csv' col = get_sensitivities(csvname) if self.J is None: self.J = np.zeros((len(col), len(_obj_names))) self.J[:, i] = np.array(col)
def provideJ(self): """ Create jacobian from adjoint results.""" if self._first_lin: self._first_lin = False else: self.config_in.RESTART_SOL = "YES" self.J = None for i, name in enumerate(_obj_names): self.config_in.ADJ_OBJFUNC = name state = adjoint(self.config_in) restart2solution(self.config_in, state) csvname = self.config_in.SURFACE_ADJ_FILENAME+'.csv' col = get_sensitivities(csvname) if self.J is None: self.J = np.zeros((len(col),len(_obj_names))) self.J[:,i] = np.array(col)