def METHOD_NAME(self):
    """
    Test linear solve and function evaluations
    """
    # solve
    self.prob.run_model()
    # Test functions values against historical values
    for func_name in self.FUNC_REFS:
        with self.subTest(function=func_name):
            np.testing.assert_allclose(
                self.prob[func_name], self.FUNC_REFS[func_name],
                rtol=self.rtol,
                atol=self.atol,
            )