Update cecj script and values

Author: PhilipFackler

tests/unit/opflow/constraint_jacobian/equality/README.md

@@ -132,10 +132,10 @@ With the parameters of the example network the matrix is:
 <td>Q1</td> <td>-1.6</td> <td>-2.0</td> <td>1.6</td> <td>-0.4</td> <td>0</td> <td>0</td> <td>0</td> <td>0</td> <td>0</td> <td>0</td> <td>0</td> <td>0</td>
 </tr>
 <tr>
-<td>P2</td> <td>-0.8</td> <td>-0.8</td> <td>0.8</td> <td>2.8</td> <td>0.8</td> <td>-0.2</td> <td>0</td> <td>0</td> <td>-0.8</td> <td>-0.8</td> <td>0</td> <td>0</td>
+<td>P2</td> <td>-0.8</td> <td>-0.8</td> <td>0.8</td> <td>1.81</td> <td>0.8</td> <td>-0.2</td> <td>0</td> <td>0</td> <td>-0.8</td> <td>-0.8</td> <td>0</td> <td>0</td>
 </tr>
 <tr>
-<td>Q2</td> <td>1.6</td> <td>-0.4</td> <td>-3.6</td> <td>-3.8</td> <td>0.4</td> <td>0.4</td> <td>0</td> <td>0</td> <td>1.6</td> <td>-0.4</td> <td>0</td> <td>0</td>
+<td>Q2</td> <td>1.6</td> <td>-0.4</td> <td>-3.6</td> <td>-3.998</td> <td>0.4</td> <td>0.4</td> <td>0</td> <td>0</td> <td>1.6</td> <td>-0.4</td> <td>0</td> <td>0</td>
 </tr>
 <tr>
 <td>P3</td> <td>0</td> <td>0</td> <td>-0.8</td> <td>0.2</td> <td>0.8</td> <td>1.8</td> <td>-1</td> <td>0</td> <td>0</td> <td>0</td> <td>0</td> <td>0</td>

tests/unit/opflow/constraint_jacobian/equality/Untitled.m

@@ -11,6 +11,7 @@
 theta4=deg2rad(0);
 theta5=deg2rad(0);
 
+mvabase = 100;
 
 R=2;
 X=1;
@@ -28,8 +29,8 @@
 Pd=-3.4;
 Qd=8.8;
 
-Gl=0.25;
-Bl=-0.05;
+Gl=0.25 / mvabase;
+Bl=-0.05 / mvabase;
 
 
 theta12=theta1-theta2;

tests/unit/opflow/constraint_jacobian/equality/cecj.csv

@@ -10,15 +10,15 @@
 3,1,-0.8
 3,2,-0.8
 3,3,0.8
-3,4,2.8
+3,4,1.81
 3,5,0.8
 3,6,-0.2
 3,9,-0.8
 3,10,-0.8
 4,1,1.6
 4,2,-0.4
 4,3,-3.6
-4,4,-3.8
+4,4,-3.998
 4,5,0.4
 4,6,0.4
 4,9,1.6

tests/unit/opflow/constraint_jacobian/equality/jac_eq_acopf.cpp

@@ -89,7 +89,7 @@ int main(int argc, char **argv) {
   PetscCall(VecView(X, PETSC_VIEWER_STDOUT_SELF));
   PetscCall(MatView(J_eq, PETSC_VIEWER_STDOUT_SELF));
 
-  PetscCall(MatAXPY(J_eq, -1.0, J_eq_ref, SAME_NONZERO_PATTERN));
+  PetscCall(MatAXPY(J_eq, -1.0, J_eq_ref, UNKNOWN_NONZERO_PATTERN));
   PetscReal norm = 0.0;
   PetscCall(MatNorm(J_eq, NORM_INFINITY, &norm));
   std::cout << "Error norm: " << norm << std::endl;