fix terminal id

data/dsntp/ieeeds/case13/dev/meas_def_ieee13.csv

 Index,ID,Point ID,Dev ID,Terminal ID,Desc
 Index,ID,Point ID,Dev ID,Terminal ID,Desc
-1,1,1000001,1,650,TOTAL
-2,2,1000002,5,634,A
-3,3,1000003,5,634,A
-4,4,1000004,5,634,B
-5,5,1000005,5,634,B
-6,6,1000006,5,634,C
-7,7,1000007,5,634,C
-8,8,1000008,3,645,B
-9,9,1000009,3,645,B
-10,10,1000010,6,646,B
-11,11,1000011,6,646,B
-12,12,1000012,8,652,A
-13,13,1000013,8,652,A
-14,14,1000014,9,671,A
-15,15,1000015,9,671,A
-16,16,1000016,9,671,B
-17,17,1000017,9,671,B
-18,18,1000018,9,671,C
-19,19,1000019,9,671,C
-20,20,1000020,14,675,A
-21,21,1000021,14,675,A
-22,22,1000022,14,675,B
-23,23,1000023,14,675,B
-24,24,1000024,14,675,C
-25,25,1000025,14,675,C
-26,26,1000026,12,692,C
-27,27,1000027,12,692,C
-28,28,1000028,13,611,C
-29,29,1000029,13,611,C
+1,1,1000001,1,1,TOTAL
+2,2,1000002,15,26,A
+3,3,1000003,15,26,A
+4,4,1000004,15,26,B
+5,5,1000005,15,26,B
+6,6,1000006,15,26,C
+7,7,1000007,15,26,C
+8,8,1000008,16,27,B
+9,9,1000009,16,27,B
+10,10,1000010,17,28,B
+11,11,1000011,17,28,B
+12,12,1000012,18,29,A
+13,13,1000013,18,29,A
+14,14,1000014,19,30,A
+15,15,1000015,19,30,A
+16,16,1000016,19,30,B
+17,17,1000017,19,30,B
+18,18,1000018,19,30,C
+19,19,1000019,19,30,C
+20,20,1000020,20,31,A
+21,21,1000021,20,31,A
+22,22,1000022,20,31,B
+23,23,1000023,20,31,B
+24,24,1000024,20,31,C
+25,25,1000025,20,31,C
+26,26,1000026,21,32,C
+27,27,1000027,21,32,C
+28,28,1000028,22,33,C
+29,29,1000029,22,33,C
\ No newline at end of file

mems/examples/ds-3phase-pf/src/nlp.rs

@@ -113,9 +113,37 @@ fn get_pq_of_acline(r_x: Array<Complex64, Ix2>) -> Option<(String, String)> {
                 r_x[[1, 0]], r_x[[1, 1]], r_x[[1, 2]],
                 r_x[[1, 0]], r_x[[1, 1]], r_x[[1, 2]],
                 r_x[[2, 0]], r_x[[2, 1]], r_x[[2, 2]]);
                 r_x[[2, 0]], r_x[[2, 1]], r_x[[2, 2]]);
             let gb = rx.try_inverse().unwrap();
             let gb = rx.try_inverse().unwrap();
-            format!("{}*x1-{}*x2-{}*x3", gb.m11, gb.m12, gb.m13)
+            format!("{:.4}*x1-{:.4}*x2-{:.4}*x3", gb.m11.re, gb.m12.re, gb.m13.re)
         }
         }
         _ => { return None; }
         _ => { return None; }
     };
     };
     Some((result.clone(), result.clone()))
     Some((result.clone(), result.clone()))
+}
+
+// test
+#[cfg(test)]
+mod test {
+    use ndarray::array;
+    use super::*;
+
+    #[test]
+    fn test_get_pq_of_acline() {
+        // 原矩阵：
+        // 0.3465+1.0179j  0.1560+0.5017j  0.1580+0.4236j
+        // 0.1560+0.5017j  0.3375+1.0478j  0.1535+0.3849j
+        // 0.1580+0.4236j  0.1535+0.3849j  0.3414+1.0348j
+        // 求逆的结果：
+        //   0.4338 - 1.2502i  -0.1840 + 0.4622i  -0.1008 + 0.3455i
+        //   -0.1840 + 0.4622i   0.3798 - 1.1847i  -0.0478 + 0.2639i
+        //   -0.1008 + 0.3455i  -0.0478 + 0.2639i   0.3359 - 1.1176i
+        let arr = array![  [Complex64::new(0.3465,1.0179), Complex64::new(0.1560,0.5017), Complex64::new(0.1580,0.4236)],
+                                            [Complex64::new(0.1560,0.5017), Complex64::new(0.3375,1.0478), Complex64::new(0.1535,0.3849)],
+                                            [Complex64::new(0.1580,0.4236), Complex64::new(0.1535,0.3849), Complex64::new(0.3414,1.0348)]];
+        let (p, q) = get_pq_of_acline(arr).unwrap();
+        assert_eq!(p, "0.4338*x1--0.1840*x2--0.1008*x3");
+        let arr = array![  [Complex64::new(0.0,0.0), Complex64::new(0.0,0.0), Complex64::new(0.0,0.0)],
+                                            [Complex64::new(0.0,0.0), Complex64::new(0.0,0.0), Complex64::new(0.0,0.0)],
+                                            [Complex64::new(0.0,0.0), Complex64::new(0.0,0.0), Complex64::new(0.3414,1.0348)]];
+        let (p, q) = get_pq_of_acline(arr).unwrap();
+    }
 }
 }
\ No newline at end of file