def execute(self): # local copy if self._first_exec: self._first_exec = False else: self.config_in.RESTART_SOL = "YES" state = direct(self.config_in) restart2solution(self.config_in, state) for name in _obj_names: setattr(self, name, state.FUNCTIONS[name])
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)
def execute(self): # local copy state = direct(self.config_in) restart2solution(self.config_in, state) for name in _obj_names: setattr(self, name, state.FUNCTIONS[name])