Updates to CMake and code style.

Author: pelesh

examples/CMakeLists.txt

@@ -13,6 +13,7 @@ if(TARGET SUNDIALS::idas)
   add_subdirectory(RLCircuit)
   add_subdirectory(Microgrid)
   add_subdirectory(ScaleMicrogrid)
+  add_subdirectory(PhasorDynamics)
   if(GRIDKIT_ENABLE_IPOPT)
     add_subdirectory(DynamicConstrainedOpt)
     add_subdirectory(GenConstLoad)
@@ -24,5 +25,3 @@ endif()
 if(GRIDKIT_ENABLE_ENZYME)
   add_subdirectory(Enzyme)
 endif()
-
-add_subdirectory(PhasorDynamics)

src/Model/PhasorDynamics/BusFault/BusFault.hpp

@@ -64,10 +64,11 @@ namespace GridKit
       {
         if (status_)
         {
-          double B  = -X_ / (X_ * X_ + R_ * R_);
-          double G  = R_ / (X_ * X_ + R_ * R_);
-          Ir()     += -Vr() * G + Vi() * B;
-          Ii()     += -Vr() * B - Vi() * G;
+          double B = -X_ / (X_ * X_ + R_ * R_);
+          double G = R_ / (X_ * X_ + R_ * R_);
+
+          Ir() += -Vr() * G + Vi() * B;
+          Ii() += -Vr() * B - Vi() * G;
         }
         return 0;
       }

src/Model/PhasorDynamics/SynchronousMachine/GENROUwS/README.md

@@ -93,9 +93,14 @@ then
 Sat(\psi'')=Sat(\vert V_{r}+jV_{i} \vert)
 ```
 
-It is important to point out that finding the initial value of $`\delta`$ for the model without saturation direct method can be used. In case when saturation is considered some "claver" math is needed. Key insight for determining initial $`\delta`$ is that the magnitude of the saturation depends upon the magnitude of $`\psi''`$, which is independent of $`\delta`$.
+It is important to point out that finding the initial value of $`\delta`$ for
+the model without saturation direct method can be used. In case when saturation
+is considered some "claver" math is needed. Key insight for determining initial
+$`\delta`$ is that the magnitude of the saturation depends upon the magnitude
+of $`\psi''`$, which is independent of $`\delta`$.
 ```math
-\delta=tan^{-1}(\dfrac{K_{sat}V_{iterm}+K_{sat}R_{a}I_{i}+(K_{sat}X''_{d}+X_{q}-X''_{q})I_{r}}{K_{sat}V_{rterm}+K_{sat}R_{a}I_{r}-(K_{sat}X''_{d}+X_{q}-X''_{q})I_{i}})
+\delta=\tan^{-1}\left(\dfrac{K_{sat}V_{iterm}+K_{sat}R_{a}I_{i}+(K_{sat}X''_{d}+X_{q}-X''_{q})I_{r}}
+                            {K_{sat}V_{rterm}+K_{sat}R_{a}I_{r}-(K_{sat}X''_{d}+X_{q}-X''_{q})I_{i}} \right)
 ```
 where
 ```math