Fixing inconsistent results when using seed in examples (reopened)

[sscini]

Fixes # .

When using the reactor_design leave_N_out example in parmest, the results from multiple runs of the example vary a lot, with random seed defined.

Summary/Motivation:

Using a random seed should allow for reproducible results. Currently, the results vary in the leave_N_out reactor design example (pyomo/pyomo/contrib/parmest/examples/reactor_design/leaveNout_example.py) for parmest even with an active random seed. I included a temporary .py file (pyomo/pyomo/contrib/parmest/examples/reactor_design/mwe_bootstrap_error_SSE.py) with the minimum working example showing the results vary with multiple runs.

Changes proposed in this PR:

• Modify leaveNout_bootstrap_test to produce consistent outputs

TODO before converting from draft:

• Inform pyomo team
• Implement solution
• Testing

Legal Acknowledgement

By contributing to this software project, I have read the contribution guide and agree to the following terms and conditions for my contribution:

1. I agree my contributions are submitted under the BSD license.
2. I represent I am authorized to make the contributions and grant the license. If my employer has rights to intellectual property that includes these contributions, I represent that I have received permission to make contributions and grant the required license on behalf of that employer.

[sscini]

Continuation of PR # 3608. Changed the branch name to be more consistent with goal, did not realize that would close PR. Thank you @blnicho for feedback and approval of change. Will make requested changes, make sure tests don't break, and change from draft to full PR.

[sscini]

@adowling2 @djlaky Ready for review, changes made examples reproducible with a set seed. Seeds also added to test_file. Test results:
Ran 44 tests in 14.982s
OK (skipped=2)

[adowling2]

Extra comma on this line?

[adowling2]

Alternate idea: either here or in the examples file (or somewhere else the Pyomo team recommends), define

_RANDOM_SEED_FOR_TESTING = 524

Then in all the files that use this random seed, you can import this value.

Then in this file:

if _RANDOM_SEED_FOR_TESTING:
   np.random.seed(_RANDOM_SEED_FOR_TESTING)

This will mean someone only sets the testing random seed in one spot.

[sscini]

Plan for this week:
Address comments above, confirm tests run, add additional comments, and run black on code for formatting

[sscini]

Notes for Pyomo Dev Meeting 6/17/25:

• Plan to complete plan from previous and hit 'Ready for Review' to convert to normal PR.

Context
Bootstrap, leaveNout, and confidence region examples were not reproducible with a set seed.

Problem
Global seed was set, but not utilized for certain functions, and was not present in graphics.py for any of the scipy commands

Implemented solution
Added global seeds to different functions and in test file, working to optimize and minimize how many were added.

Questions

• Any recommendations for previous comment from Dr. Dowling?
• Would any new documentation or additional tests be needed to make sure this would run properly?