From 71559f3e2498051b9205242defdf994e30838c75 Mon Sep 17 00:00:00 2001
From: wyq <yaqiang.wang@gmail.com>
Date: Tue, 8 Apr 2025 09:02:24 +0800
Subject: [PATCH] add toeplitz function

---
 meteoinfo-lab/milconfig.xml                   |  16 ++--
 .../pylib/mipylib/numeric/core/_ndarray.py    |   4 +-
 .../mipylib/numeric/linalg/__init__$py.class  | Bin 2177 -> 2734 bytes
 .../pylib/mipylib/numeric/linalg/__init__.py  |   7 +-
 .../numeric/linalg/_special_matrices.py       |  76 ++++++++++++++++++
 .../org/meteoinfo/ndarray/math/ArrayUtil.java |  25 ++++++
 6 files changed, 117 insertions(+), 11 deletions(-)
 create mode 100644 meteoinfo-lab/pylib/mipylib/numeric/linalg/_special_matrices.py

diff --git a/meteoinfo-lab/milconfig.xml b/meteoinfo-lab/milconfig.xml
index 664636d2..e810f28e 100644
--- a/meteoinfo-lab/milconfig.xml
+++ b/meteoinfo-lab/milconfig.xml
@@ -1,34 +1,32 @@
 <?xml version="1.0" encoding="UTF-8" standalone="no"?>
 <MeteoInfo File="milconfig.xml" Type="configurefile">
-  <Path OpenPath="D:\Working\MIScript\Jython\mis\map\webmap">
-    <RecentFolder Folder="D:\Temp\test"/>
+  <Path OpenPath="D:\Working\MIScript\Jython\mis\array">
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\common_math\stats"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\io\dataconvert"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\io"/>
-    <RecentFolder Folder="D:\Working\MIScript\Jython\mis\array"/>
-    <RecentFolder Folder="D:\Working\MIScript\Jython\mis\common_math"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\common_math\signal"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\test"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\io\netcdf"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\dataframe"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\io\burf"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\io\grib"/>
-    <RecentFolder Folder="D:\Working\MIScript\Jython\mis"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\map"/>
     <RecentFolder Folder="D:\Working\MIScript\Jython\mis\map\webmap"/>
+    <RecentFolder Folder="D:\Working\MIScript\Jython\mis\common_math\linalg"/>
+    <RecentFolder Folder="D:\Working\MIScript\Jython\mis\common_math"/>
+    <RecentFolder Folder="D:\Working\MIScript\Jython\mis"/>
+    <RecentFolder Folder="D:\Working\MIScript\Jython\mis\array"/>
   </Path>
   <File>
     <OpenedFiles>
       <OpenedFile File="D:\Working\MIScript\Jython\mis\io\radar\radar_x_phase_2.py"/>
-      <OpenedFile File="D:\Working\MIScript\Jython\mis\io\grib\grib2_ecmwf_3.py"/>
       <OpenedFile File="D:\Working\MIScript\Jython\mis\map\webmap\new_tilefactory_tianditu-1.py"/>
-      <OpenedFile File="D:\Working\MIScript\Jython\mis\map\webmap\air_path.py"/>
+      <OpenedFile File="D:\Working\MIScript\Jython\mis\common_math\linalg\toeplitz_1.py"/>
     </OpenedFiles>
     <RecentFiles>
       <RecentFile File="D:\Working\MIScript\Jython\mis\io\radar\radar_x_phase_2.py"/>
-      <RecentFile File="D:\Working\MIScript\Jython\mis\io\grib\grib2_ecmwf_3.py"/>
       <RecentFile File="D:\Working\MIScript\Jython\mis\map\webmap\new_tilefactory_tianditu-1.py"/>
-      <RecentFile File="D:\Working\MIScript\Jython\mis\map\webmap\air_path.py"/>
+      <RecentFile File="D:\Working\MIScript\Jython\mis\common_math\linalg\toeplitz_1.py"/>
     </RecentFiles>
   </File>
   <Font>
diff --git a/meteoinfo-lab/pylib/mipylib/numeric/core/_ndarray.py b/meteoinfo-lab/pylib/mipylib/numeric/core/_ndarray.py
index 6f6fcb63..db84cc85 100644
--- a/meteoinfo-lab/pylib/mipylib/numeric/core/_ndarray.py
+++ b/meteoinfo-lab/pylib/mipylib/numeric/core/_ndarray.py
@@ -765,7 +765,7 @@ class NDArray(object):
             else:
                 ArrayMath.setImage(self._array, val)
 
-    def conj(self):
+    def conjugate(self):
         """
         Return the complex conjugate, element-wise.
 
@@ -775,6 +775,8 @@ class NDArray(object):
         """
         return self.array_wrap(ArrayMath.conj(self._array))
 
+    conj = conjugate
+
     def min(self, axis=None):
         """
         Get minimum value along an axis.
diff --git a/meteoinfo-lab/pylib/mipylib/numeric/linalg/__init__$py.class b/meteoinfo-lab/pylib/mipylib/numeric/linalg/__init__$py.class
index f273ffda73029b450873e134e6b882eb206457ca..a74e6d2a4d6989d5662b6e6c3b808790d64b4881 100644
GIT binary patch
delta 1003
zcmZuuYflqV5Ixgv>2_HnU4ft|peSGq6h%=HpB0riqM(Sfu4RF0OQBti^@YYy#>Cj%
z_$y2dk60Rs`~VX__$y5O2_~X*X-zb7KkS^DIrrQ%bD#U)musIlKW+dxgy+kei$Z=g
zHD^bnR?LiBPBI#?Q!4y~ibQfcG?#W}60@O5B58+)(tXjC;}x~E%ddCE=bZH5ox66#
zNfD~XGrc>-1eGh*jS|)!7}ZgQT^h=;o8UD~D;6_N6*UAkWjnD%#L9dLWG1{jBy$h;
zYH*{DnbWp2YsGCsOXlnLOixvbj{T^Y-3CHsrmwQWivxsW(`2UQI7!pg5tNJong|{<
zYE4cO+8fygAK$B^1tHmNW$`I9W~H1cfx2}bzs<+8heZd~tL-@ILj{hBdFu3^ASjLb
z2~KM0#3`{s;o@^T&fu&MVVoDbt99E&9;Gv$m|Te2r&M%{Ue|l07kwJb^y?s8kkmm@
ztaO+2(JFLkob#}kF)Yz_5mjpFhIpxrsyA^<!#F0yH|4x~8<y0#Bih_1O$p(Ss56La
zn8sbP>5eO5%!*l04b6$yo?-Rvk2hsgglduTl#8F9SM*T4EE-u1vu(mWTu`NJe(g$W
ziDMB*;yU3mXYZi!4uBUycNVIl1>HH620aoRdF*Apj<T^8>|8_jK+vdd&SD=81XXTu
zmPI2MH4YLtH=9`;99S&Sf~7@&vy6XPe@LzmII@D{OYDR2g#Q)rY+<lZZ26P9bk1um
zI6s6E_72ns)^IwoiY`U@giHNRL8C`e7L^q)5Gnpr<S?*`p=G8@dES0Opp+wV&Aoyl
zc{Sj}6<jTJYgt@&%hnMdGe*{MEsL=n%s>v2yp|7>unTNsi^~>`a1QsD*cyL?7~+y3
zzI*HE6HAmI@F}<2;rhid;`Sm;UKt)iVVsA*<eU7r!|Pb&(=Md(jN@)5Ci#L?HykWT
M(1ZJUAVCW0-xzkfWdHyG

delta 649
zcmZut+fEZv6kU7Hv~xO++D<x33Z;Nfg$^iI#mY?tyc7i$rCv}gA|@tOV$j4lqOU%l
z#IHb7N`whSKfurM0}RhTi2JmCH0NQhv-e(WUpBmt3H#^Y-5&t*SlqM?N?NbV%L}XT
zUK6P3AOSZ9Qw>p@G-p=cED=)bH>K1k@>G`?Or!~-yQk!!6&V}t=uj0gap<IjE}V$L
zM^^2L{-ZqvIb2zJ`F?q6R3cD+#D?34(>9Lbi~~ZqE<LMu&7#NBR<zmtFE}`d^O_B+
zmPm$%RUuN6qquBi1Y_z;WKxb}LVHZAUy(B1P~*`>Ifa`xinyh|Mk~U{9hJ5+G^<9e
zS^1%JeSe?OuEs1+ty`bzi7LpsHJ`^O+=BtBr{+_%=#dzsjKp!mGye5K@B_TS+=g`R
zzzh)&jBiNRkg9TmP;QDeC%s2gn#bX~G@IK&XKo9pg!qo!bZ_AHiZvl>93T>jI70Mq
zA-~Bs+Je__8pIiiV`dF`-3o}|0tTAd4mQ`><R1J-?!_Hk+Qv|bE4dI?8(uM7!}Vsh
zTd(EqTtCEg^+?^0FoWB=KpnTUZ|_#AvCwl)4!}6zlb9BKvNSw}VBVky-1`3+zS{z;
c4b0;OW10)^bA!1eN_e2j6dvNSCaakL2Q_+NzW@LL

diff --git a/meteoinfo-lab/pylib/mipylib/numeric/linalg/__init__.py b/meteoinfo-lab/pylib/mipylib/numeric/linalg/__init__.py
index d677ee8f..414ac93b 100644
--- a/meteoinfo-lab/pylib/mipylib/numeric/linalg/__init__.py
+++ b/meteoinfo-lab/pylib/mipylib/numeric/linalg/__init__.py
@@ -1 +1,6 @@
-from .linalg import *
\ No newline at end of file
+from .linalg import *
+from ._special_matrices import *
+
+__all__ = []
+__all__ += linalg.__all__
+__all__ += _special_matrices.__all__
diff --git a/meteoinfo-lab/pylib/mipylib/numeric/linalg/_special_matrices.py b/meteoinfo-lab/pylib/mipylib/numeric/linalg/_special_matrices.py
new file mode 100644
index 00000000..a1d5a1ad
--- /dev/null
+++ b/meteoinfo-lab/pylib/mipylib/numeric/linalg/_special_matrices.py
@@ -0,0 +1,76 @@
+from org.meteoinfo.ndarray.math import ArrayUtil
+
+from .. import core as np
+
+
+__all__ = ['toeplitz']
+
+
+# -----------------------------------------------------------------------------
+#  matrix construction functions
+# -----------------------------------------------------------------------------
+
+
+def toeplitz(c, r=None):
+    r"""
+    Construct a Toeplitz matrix.
+
+    The Toeplitz matrix has constant diagonals, with c as its first column
+    and r as its first row. If r is not given, ``r == conjugate(c)`` is
+    assumed.
+
+    Parameters
+    ----------
+    c : array_like
+        First column of the matrix.
+    r : array_like, optional
+        First row of the matrix. If None, ``r = conjugate(c)`` is assumed;
+        in this case, if c[0] is real, the result is a Hermitian matrix.
+        r[0] is ignored; the first row of the returned matrix is
+        ``[c[0], r[1:]]``.
+
+        .. warning::
+
+            Beginning in SciPy 1.17, multidimensional input will be treated as a batch,
+            not ``ravel``\ ed. To preserve the existing behavior, ``ravel`` arguments
+            before passing them to `toeplitz`.
+
+    Returns
+    -------
+    A : (len(c), len(r)) ndarray
+        The Toeplitz matrix. Dtype is the same as ``(c[0] + r[0]).dtype``.
+
+    See Also
+    --------
+    circulant : circulant matrix
+    hankel : Hankel matrix
+    solve_toeplitz : Solve a Toeplitz system.
+
+    Notes
+    -----
+    The behavior when `c` or `r` is a scalar, or when `c` is complex and
+    `r` is None, was changed in version 0.8.0. The behavior in previous
+    versions was undocumented and is no longer supported.
+
+    Examples
+    --------
+    >>> from scipy.linalg import toeplitz
+    >>> toeplitz([1,2,3], [1,4,5,6])
+    array([[1, 4, 5, 6],
+           [2, 1, 4, 5],
+           [3, 2, 1, 4]])
+    >>> toeplitz([1.0, 2+3j, 4-1j])
+    array([[ 1.+0.j,  2.-3.j,  4.+1.j],
+           [ 2.+3.j,  1.+0.j,  2.-3.j],
+           [ 4.-1.j,  2.+3.j,  1.+0.j]])
+
+    """
+    c = np.asarray(c)
+    if r is None:
+        r = c.conjugate()
+    else:
+        r = np.asarray(r)
+
+    ta = ArrayUtil.toeplitz(c._array, r._array)
+
+    return np.NDArray(ta)
diff --git a/meteoinfo-ndarray/src/main/java/org/meteoinfo/ndarray/math/ArrayUtil.java b/meteoinfo-ndarray/src/main/java/org/meteoinfo/ndarray/math/ArrayUtil.java
index 0cada629..073fc5a8 100644
--- a/meteoinfo-ndarray/src/main/java/org/meteoinfo/ndarray/math/ArrayUtil.java
+++ b/meteoinfo-ndarray/src/main/java/org/meteoinfo/ndarray/math/ArrayUtil.java
@@ -1426,6 +1426,31 @@ public class ArrayUtil {
         return Array.factory(DataType.INT, new int[]{n, n}, M);
     }
 
+    /**
+     * Returns a Toeplitz matrix in which each descending diagonal
+     * from left to right is constant.
+     *
+     * @param col First column of the matrix. {@code A[i, j] = kl[i - j]} for {@code i >  j}
+     * @param row First row of the matrix. {@code A[i, j] = ku[j - i]} for {@code i <= j}
+     * @return the Toeplitz matrix.
+     */
+    public static Array toeplitz(Array col, Array row) {
+        int nx = (int) row.getSize();
+        int ny = (int) col.getSize();
+        DataType dataType = ArrayMath.commonType(col.getDataType(), row.getDataType());
+        Array r = Array.factory(dataType, new int[]{ny, nx});
+        for (int i = 0; i < ny; i++) {
+            for (int j = 0; j < i; j++) {
+                r.setObject(i * nx + j, col.getObject(i - j));
+            }
+            for (int j = i; j < nx; j++) {
+                r.setObject(i * nx + j, row.getObject(j - i));
+            }
+        }
+
+        return r;
+    }
+
     /**
      * Repeat a value n times
      *