Hawaii 37-bus case simulation

[nkoukpaizan]

Description

This is a functional implementation of the Hawaii case
Closes #336

Proposed changes

• Added a workaround to add Jacobian contributions from components to support parallel branches.
• Implemented a different formulation of limits (see Question: behavior of sigmoid and Indicators #320), which removes nan values. Updated the documentation accordingly. Both documentation and implementation could change later if another solution comes up.
• Added monitors to hawaii.json to generate a mon.csv file.
• Updated PhasorDynamics monitoring to remove default stdout outputs in csv format. This add a lot to the outputs on larger systems.

Checklist

• All tests pass.
• Code compiles cleanly with flags -Wall -Wpedantic -Wconversion -Wextra.
• The new code follows GridKit™ style guidelines.
• [N/A] There are unit tests for the new code.
• The new code is documented.
• The feature branch is rebased with respect to the target branch.
• [N/A] I have updated CHANGELOG.md to reflect the changes in this PR. If this is a minor PR that is part of a larger fix already included in the file, state so.

Further comments

[pelesh]

Looks good to me and tests pass.

@abirchfield and @lukelowry: Please review smoothing functions docs and implementation. I would appreciate your approval before we merge this.

[lukelowry]

Great work this is great step forward for cases

[lukelowry]

The indicator_low and indicator_high functions make sense, but I am confused by indicator_zero. Can this be brought into indicator? If Enzyme needs it this way, then that is okay by me.

[nkoukpaizan]

It's not needed by Enzyme, and it could be brought into indicator, I'm not sure what difference that would make. It made sense to me to write it this way, but I'm not particularly attached to the implementation.

[lukelowry]

This is cleaner than my original approach! Can you explain why it is called $\phi_0$?

[nkoukpaizan]

It turns out that it corresponds to f=0, so I called it $\phi_0$. The way I derived it was drawing the full function for each of the f=0, f<0 and f>0 cases. I'm not particularly attached to that notation, though.