Add test for #5007 regression

Co-authored-by: kc611 <[email protected]>

pymc/tests/test_logprob.py

@@ -190,3 +190,16 @@ def test_logcdf_helper():
 
     x_logp = logcdf(x, [0, 1], sum=False)
     np.testing.assert_almost_equal(x_logp.eval(), sp.norm(0, 1).logcdf([0, 1]))
+
+
+def test_model_unchanged_logprob_access_():
+    # Issue #5007
+    with Model() as model:
+        a = Normal("a")
+        c = Uniform("c", lower=a - 1, upper=1)
+
+    original_inputs = set(aesara.graph.graph_inputs([c]))
+    # Extract model.logpt
+    model.logpt
+    new_inputs = set(aesara.graph.graph_inputs([c]))
+    assert original_inputs == new_inputs