def solve_linear(self, d_outputs, d_residuals, mode): if mode == 'fwd': d_outputs['d'] = lu_solve(self.lu, d_residuals['d'], trans=0) else: d_residuals['d'] = lu_solve(self.lu, d_outputs['d'], trans=1)