def METHOD_NAME(self, inputs, outputs, residuals):
    num_rhs = self.options['num_rhs']
    self.K = self.assemble_CSC_K(inputs)
    for j in range(num_rhs):
        force_vector = np.concatenate([self.options['force_vector'][:, j], np.zeros(2)])
        residuals['d_%d' % j] = self.K.dot(outputs['d_%d' % j]) - force_vector