def solve_nonlinear(self, inputs, outputs): force_vector = np.concatenate([self.metadata['force_vector'], np.zeros(2)]) self.lu = lu_factor(inputs['K']) outputs['d'] = lu_solve(self.lu, force_vector)