整理思路：http://www.360doc.com/content/21/0531/20/21578602_979870548.shtml
- 更新时间：2022.09.12
- version：v1.0.0
一、下载与安装说明
Numahk库
提取码：
mono复制
解压码：无
建议使用pia安装：（参考链接）
pia install print
pia install ahktype
pia install xa
pia install numahk源码：
; Author: Mono
; Time: 2022.09.08
; Version: 1.0.0
#Include <ahktype\ahktype>
#Include <xa\xa>
; DefProp
DefProp(Numahk.NDArray().Base, "T", {Get: NumahkTranspose})
NumahkAdd(this)
{
Return this.Add
}
NumahkTranspose(this)
{
Return Numahk.Transpose(this)
}
Class Numahk
{
Class NDArray Extends Array
{
; Basic
__New(Obj*)
{
if !Obj.Length
Return
Obj.Standardization()
if Obj.NDim == 1
this.Extend(Obj)
else
{
Obj := Obj[1]
For i in Obj
{
if i.NDim == 0
this.Extend(i)
else
this.Push(Numahk.NDArray(i))
}
}
}
__Item[Pos*]
{
Get => this.GetMethod(Pos*)
Set => this.SetMethod(Value, True, Pos*)
}
GetMethod(Lst_Pos*)
{
if Lst_Pos.NDim == 1
Return this.GetSingleMethod(Lst_Pos*)
Tmp := Lst_Pos.Shape
Tmp.Pop()
Lst_Pos := Numahk.NDArray(Lst_Pos)
New_Pos := []
For i in Tmp
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
TmpNDArray := Numahk.Zeros(Tmp)
For i in New_Pos
TmpNDArray.SetMethod(this.GetSingleMethod(Lst_Pos[i*]*), False, i*)
Return TmpNDArray
}
GetSingleMethod(Pos*)
{
Tmp := [this]
Pos := Pos.Clone()
While Pos.Length
{
TmpPos := Pos.RemoveAt(1)
Tmpthis := Tmp.Pop()
For i in Tmpthis
{
if TmpPos == A_Index
Tmp.Push(i)
}
}
Return Tmp[-1]
}
SetMethod(Value, Flag, Pos*)
{
Value := (Value.NDim) ? Numahk.NDArray(Value) : Value
Pos := Pos.Clone()
LastPos := Pos.Pop()
Tmp := this.GetSingleMethod(Pos*)
if ListCmp(Tmp[LastPos].Shape, Value.Shape) && Flag
Throw ValueError("SetMethod Error!`nTarget shape is " Print.ToString(Tmp[LastPos].Shape) ".`nBut got an object with shape " Print.ToString(Value.Shape) ".")
Tmp.RemoveAt(LastPos)
Tmp.InsertAt(LastPos, Value)
}
; Operator
Add(NDArray := 0)
{
Return Numahk.Add(this, NDArray)
}
Div(NDArray := 1)
{
Return Numahk.Divide(this, NDArray)
}
Mul(NDArray := 1)
{
Return Numahk.Multiply(this, NDArray)
}
Sub(NDArray := 0)
{
Return Numahk.Subtract(this, NDArray)
}
; Function
All()
{
Return Numahk.All()
}
AMax(Axis := "axis=None")
{
Return Numahk.Max(this, Axis)
}
AMin(Axis := "axis=None")
{
Return Numahk.Min(this, Axis)
}
Any()
{
Return Numahk.Any()
}
ArgMax(Axis := "axis=None")
{
Return Numahk.ArgMax(this, Axis)
}
ArgMin(Axis := "axis=None")
{
Return Numahk.ArgMin(this, Axis)
}
CumProd(Axis := "axis=None")
{
Return Numahk.CumProd(this, Axis)
}
CumSum(Axis := "axis=None")
{
Return Numahk.CumSum(this, Axis)
}
Dot(NDArray*)
{
Return Numahk.Dot(this, NDArray*)
}
Max(Axis := "axis=None")
{
Return Numahk.Max(this, Axis)
}
Mean(Axis := "axis=None")
{
Return Numahk.Mean(this, Axis)
}
Min(Axis := "axis=None")
{
Return Numahk.Min(this, Axis)
}
Reshape(Shape*)
{
Return Numahk.NDArray(Super.Reshape(Shape*))
}
Sort(Axis := "axis=None")
{
Return Numahk.Sort(this, Axis)
}
Std(Axis := "axis=None")
{
Return Numahk.Std(this, Axis)
}
Sum(Axis := "axis=None")
{
Return Numahk.Sum(this, Axis)
}
SwapAxes(Axis1, Axis2)
{
Tmp := Range(0, this.NDim - 1).Array()
Tmp.Swap(Axis1 + 1, Axis2 + 1)
Return this.Transpose(Tmp)
}
Transpose(Shape*)
{
if !Shape.Length
Return this.T
if Shape.NDim == 2
Shape := Shape[1]
New_Shape := Shape.Clone()
Tmp := this.Shape
For i in Shape
New_Shape[A_Index] := Tmp[Shape.Index(A_Index - 1)]
New_NDArray := Numahk.Zeros(New_Shape)
Pos_Old := []
Pos_New := []
Loop Tmp.Length
Pos_Old.Push(Range(1, Tmp[A_Index]).Array())
Pos_Old := Dynloop_Loop(Pos_Old)
For i in Pos_Old
{
k := i.Clone()
Loop k.Length
k[A_Index] := i[Shape.Index(A_Index - 1)]
Pos_New.Push(k)
}
Loop Pos_Old.Length
New_NDArray[Pos_New[A_Index]*] := this[Pos_Old[A_Index]*]
Return New_NDArray
}
Var(Axis := "axis=None")
{
Return Numahk.Var(this, Axis)
}
}
; NoUse
Static Transpose(NDArray)
{
New_Shape := NDArray.Shape.Clone().Reverse()
New_NDArray := Numahk.Zeros(New_Shape)
Pos_Old := []
Pos_New := []
Loop NDArray.Shape.Length
Pos_Old.Push(Range(1, NDArray.Shape[A_Index]).Array())
Pos_Old := Dynloop_Loop(Pos_Old)
For i in Pos_Old
{
k := i.Clone()
Pos_New.Push(k.Reverse())
}
Loop Pos_Old.Length
New_NDArray[Pos_New[A_Index]*] := NDArray[Pos_Old[A_Index]*]
Return New_NDArray
}
; Create
Static ARange(Args*)
{
Return Numahk.NDArray(Range(Args*).Array())
}
Static Array(Obj*)
{
Return Numahk.NDArray(Obj*)
}
Static AsArray(NDArray)
{
Tmp := NDArray.Copy()
Return Numahk.NDArray(Tmp)
}
Static Empty(Shape)
{
if Type(Shape).in(["Integer", "String"])
Shape := [Shape]
Res_Lst := [""]
Loop Shape.Length
{
Res_Lst.mul(Shape[-A_Index])
if A_Index !== Shape.Length
Res_Lst := [Res_Lst]
}
Res_Lst := Numahk.NDArray(Res_Lst)
Return Res_Lst
}
Static Empty_Like(NDArray)
{
Shape := NDArray.Shape
Res_Lst := [""]
Loop Shape.Length
{
Res_Lst.mul(Shape[-A_Index])
if A_Index !== Shape.Length
Res_Lst := [Res_Lst]
}
Res_Lst := Numahk.NDArray(Res_Lst)
Return Res_Lst
}
Static Eye(N, M := 0, k := 0)
{
if M
NDArray := Numahk.Zeros([N, M])
else
NDArray := Numahk.Zeros([N, N])
Loop N
{
Try
NDArray[A_Index, A_Index + k] := 1
}
Return NDArray
}
Static Identity(N)
{
NDArray := Numahk.Zeros([N, N])
Loop N
NDArray[A_Index, A_Index] := 1
Return NDArray
}
Static In1d(NDArray1, NDArray2, Invert := False)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
NDArray := Numahk.Zeros_Like(NDArray1)
New_Pos := []
For i in NDArray1.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For j in New_Pos
{
if NDArray2.NDim
{
flag := Invert
For k in NDArray2
{
flag := Invert ? (NDArray1[j*] != k) : (NDArray1[j*] = k)
if flag != Invert
{
NDArray[j*] := flag
Continue 2
}
}
NDArray[j*] := flag
}
else
NDArray[j*] := Invert ? (NDArray1[j*] != NDArray2) : (NDArray1[j*] = NDArray2)
}
Return NDArray
}
Static MeshGrid(NDArray*)
{
}
Static Ones(Shape)
{
if Type(Shape).in(["Integer", "String"])
Shape := [Shape]
Res_Lst := [1]
Loop Shape.Length
{
Res_Lst.mul(Shape[-A_Index])
if A_Index !== Shape.Length
Res_Lst := [Res_Lst]
}
Res_Lst := Numahk.NDArray(Res_Lst)
Return Res_Lst
}
Static Ones_Like(NDArray)
{
Shape := NDArray.Shape
Res_Lst := [1]
Loop Shape.Length
{
Res_Lst.mul(Shape[-A_Index])
if A_Index !== Shape.Length
Res_Lst := [Res_Lst]
}
Res_Lst := Numahk.NDArray(Res_Lst)
Return Res_Lst
}
Static Where(Cond, NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
For i in NDArray1.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For j in New_Pos
NDArray1[j*] := ((Cond is Func) ? Cond(NDArray1[j*]) : Cond[j*]) ? NDArray1[j*] : ((NDArray2.NDim) ? NDArray2[j*] : NDArray2)
Return NDArray1
}
Static Zeros(Shape)
{
if Type(Shape).in(["Integer", "String"])
Shape := [Shape]
Res_Lst := [0]
Loop Shape.Length
{
Res_Lst.mul(Shape[-A_Index])
if A_Index !== Shape.Length
Res_Lst := [Res_Lst]
}
Res_Lst := Numahk.NDArray(Res_Lst)
Return Res_Lst
}
Static Zeros_Like(NDArray)
{
Shape := NDArray.Shape
Res_Lst := [0]
Loop Shape.Length
{
Res_Lst.mul(Shape[-A_Index])
if A_Index !== Shape.Length
Res_Lst := [Res_Lst]
}
Res_Lst := Numahk.NDArray(Res_Lst)
Return Res_Lst
}
; Matrix
Static Diag(NDArray, k := 0)
{
if NDArray.NDim == 1
{
Tmp := Numahk.Zeros([Abs(k) + NDArray.Length, Abs(k) + NDArray.Length])
Loop NDArray.Length
{
Try
Tmp[A_Index, A_Index + k] := NDArray[A_Index]
}
Return Tmp
}
Tmp := []
Loop NDArray.Length
{
Try
Tmp.Push(NDArray[A_Index, A_Index + k])
}
Return Numahk.NDArray(Tmp)
}
Static Dot(NDArray*)
{
Temp := NDArray.RemoveAt(1)
if !NDArray.Length
Return Temp
if NDArray.Length == 1
{
Temp2 := NDArray.RemoveAt(1)
Temp3 := Numahk.Zeros([Temp.Shape[1], Temp2.Shape[2]])
For i in Range(Temp.Shape[1])
{
For j in Range(Temp2.Shape[2])
{
Sum := 0
For k in Range(Temp.Shape[2])
Sum += Temp[i, k] * Temp2[k, j]
Temp3[i, j] := Round(Sum, 10)
}
}
Return Temp3
}
else
{
Temp2 := Numahk.Dot(Temp, NDArray.RemoveAt(1))
Return Numahk.Dot(Temp2, NDArray*)
}
}
Static Trace(NDArray, offset := 0)
{
New_Pos := []
Tmp := NDArray.Shape
if !(Tmp.Length - 2)
Return Numahk.Diag(NDArray, offset).Sum
Loop Tmp.Length - 2
New_Pos.InsertAt(1, Range(1, Tmp[-A_Index]).Array())
New_Pos := Dynloop_Loop(New_Pos)
Tmp2 := []
For i in New_Pos
{
Tmp3 := []
Loop Tmp.Length
{
Try
Tmp3.Push(NDArray[A_Index, A_Index + offset, i*])
}
Tmp2.Push(Tmp3.Sum)
}
Tmp.RemoveAt(1, 2)
Return Numahk.NDArray(Tmp2).Reshape(Tmp)
}
; Sort
Static Intersect1d(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
New_NDArray := []
For i in NDArray1.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Value := NDArray1[i*]
if Value.in(NDArray2) && !Value.in(New_NDArray)
New_NDArray.Push(Value)
}
Return Numahk.NDArray(New_NDArray.Sort())
}
Static Setdiff1d(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
New_NDArray := []
For i in NDArray1.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Value := NDArray1[i*]
if !Value.in(NDArray2) && !Value.in(New_NDArray)
New_NDArray.Push(Value)
}
Return Numahk.NDArray(New_NDArray.Sort())
}
Static Setxor1d(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
New_NDArray := []
For i in NDArray1.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Value := NDArray1[i*]
if !Value.in(NDArray2) && !Value.in(New_NDArray)
New_NDArray.Push(Value)
}
For i in New_Pos
{
Value := NDArray2[i*]
if !Value.in(NDArray1) && !Value.in(New_NDArray)
New_NDArray.Push(Value)
}
Return Numahk.NDArray(New_NDArray.Sort())
}
Static Sort(NDArray, Axis := "axis=-1")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
Return Numahk.NDArray(NDArray.Ravel().Sort())
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
New_NDArray := []
Tmp := NDArray.Shape[Axis + 1]
New_Pos := []
New_Shape := NDArray.Shape
New_Shape.RemoveAt(Axis + 1)
Loop New_Shape.Length
New_Pos.Push(Range(1, New_Shape[A_Index]).Array())
New_Pos := Dynloop_Loop(New_Pos)
New_PosTmp := []
For i in New_Pos
{
Loop Tmp
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
New_PosTmp.Push(iTmp)
}
}
Loop NDArray.Shape.Product // Tmp
{
Index := A_Index - 1
Tmp2 := []
Loop Tmp
Tmp2.Push(NDArray[New_PosTmp[A_Index + Index * Tmp]*])
Tmp2.Sort()
Loop Tmp
NDArray[New_PosTmp[A_Index + Index * Tmp]*] := Tmp2[A_Index]
}
Return NDArray
}
Static Union1d(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
New_NDArray := []
For i in NDArray1.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
New_NDArray.Push(NDArray1[i*])
For i in New_Pos
{
Value := NDArray2[i*]
if !Value.in(New_NDArray)
New_NDArray.Push(Value)
}
Return Numahk.NDArray(New_NDArray.Sort())
}
Static Unique(NDArray, Return_Index := False, Return_Counts := False)
{
Ret := []
New_NDArray := []
New_IndexArray := []
For i in NDArray
{
if !i.in(New_NDArray)
{
New_NDArray.Push(i)
New_IndexArray.Push(A_Index)
}
}
if !Return_Index && !Return_Counts
Return Numahk.NDArray(New_NDArray)
Ret.Push(New_NDArray)
if Return_Index
Ret.Push(New_IndexArray)
if Return_Counts
{
Tmp := []
For i in New_NDArray
Tmp.Push(NDArray.Count(i))
Ret.Push(Tmp)
}
Return Numahk.NDArray(Ret)
}
; 1-D Calculate
Static Abs(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Abs(NDArray[i*])
Return NDArray
}
Static All(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if !NDArray[i*]
Return False
}
Return True
}
Static Any(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray[i*]
Return True
}
Return False
}
Static ArcCos(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := ACos(NDArray[i*])
Return NDArray
}
Static ArcCosh(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := Ln(Tmp + Sqrt(Tmp ** 2 - 1))
}
Return NDArray
}
Static ArcSin(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := ASin(NDArray[i*])
Return NDArray
}
Static ArcSinh(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := Ln(Tmp + Sqrt(Tmp ** 2 + 1))
}
Return NDArray
}
Static ArcTan(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := ATan(NDArray[i*])
Return NDArray
}
Static ArcTanh(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := 0.5 * Ln((1 + Tmp) / (1 - Tmp))
}
Return NDArray
}
Static ArgMax(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Return Tmp.Index(Max(Tmp*))
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
Return Tmp.Index(Max(NDArray*))
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.ArgMax(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
Static ArgMin(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Return Tmp.Index(Min(Tmp*))
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
Return Tmp.Index(Min(NDArray*))
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.ArgMin(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
Static Ceil(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Ceil(NDArray[i*])
Return NDArray
}
Static Cos(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Cos(NDArray[i*])
Return NDArray
}
Static Cosh(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := (Exp(Tmp) + Exp(-Tmp)) / 2
}
Return NDArray
}
Static CumProd(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Loop Tmp.Length - 1
Tmp[A_Index + 1] *= Tmp[A_Index]
Return Tmp
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
{
Loop NDArray.Length - 1
NDArray[A_Index + 1] *= NDArray[A_Index]
Return NDArray
}
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Tmp2)
}
Loop Lst_Res.Length
{
Index := A_Index
Loop Lst_Res[Index].Length - 1
Lst_Res[Index][A_Index + 1] *= Lst_Res[Index][A_Index]
}
Lst_Res.Ravel()
Index := 1
For i in Lst_Pos
{
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
NDArray[iTmp*] := Lst_Res[Index++]
}
}
Return NDArray
}
}
Static CumSum(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Loop Tmp.Length - 1
Tmp[A_Index + 1] += Tmp[A_Index]
Return Tmp
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
{
Loop NDArray.Length - 1
NDArray[A_Index + 1] += NDArray[A_Index]
Return NDArray
}
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Tmp2)
}
Loop Lst_Res.Length
{
Index := A_Index
Loop Lst_Res[Index].Length - 1
Lst_Res[Index][A_Index + 1] += Lst_Res[Index][A_Index]
}
Lst_Res.Ravel()
Index := 1
For i in Lst_Pos
{
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
NDArray[iTmp*] := Lst_Res[Index++]
}
}
Return NDArray
}
}
Static Exp(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Exp(NDArray[i*])
Return NDArray
}
Static FAbs(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Abs(NDArray[i*])
Return NDArray
}
Static Floor(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Floor(NDArray[i*])
Return NDArray
}
Static Log(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Ln(NDArray[i*])
Return NDArray
}
Static Log10(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Log(NDArray[i*])
Return NDArray
}
Static Log1p(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Ln(NDArray[i*] + 1)
Return NDArray
}
Static Log2(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Ln(NDArray[i*]) / Ln(10)
Return NDArray
}
Static Logical_Not(NDArray, Function := 0)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := Function ? Function(Tmp) : !Tmp
}
Return NDArray
}
Static Max(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Return Max(Tmp*)
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
Return Max(NDArray*)
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.Max(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
Static Mean(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Return Tmp.Sum / Tmp.Length
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
Return NDArray.Sum / NDArray.Length
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.Mean(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
Static Min(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Return Min(Tmp*)
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
Return Min(NDArray*)
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.Min(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
Static Modf(NDArray)
{
New_Pos := []
New_NDArray := [Numahk.Zeros_Like(NDArray), Numahk.Zeros_Like(NDArray)]
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
New_NDArray[1][i*] := Integer(Tmp - Mod(Tmp, 1))
New_NDArray[2][i*] := Mod(Tmp, 1)
}
Return New_NDArray
}
Static Rint(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
if Tmp < 0
NDArray[i*] := (Abs(Tmp - Floor(Tmp)) <= 0.5) ? Floor(Tmp) : Ceil(Tmp)
else
NDArray[i*] := (Abs(Tmp - Floor(Tmp)) < 0.5) ? Floor(Tmp) : Ceil(Tmp)
}
Return NDArray
}
Static Sign(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := (Tmp > 0) ? 1 : ((Tmp < 0) ? -1 : 0)
}
Return NDArray
}
Static Sin(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Sin(NDArray[i*])
Return NDArray
}
Static Sinh(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := (Exp(Tmp) - Exp(Tmp)) / 2
}
Return NDArray
}
Static Sqrt(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Sqrt(NDArray[i*])
Return NDArray
}
Static Square(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := NDArray[i*] ** 2
Return NDArray
}
Static Std(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
avg := Tmp.Sum / Tmp.Length
sum := 0
For i in Tmp
sum += (i - avg) ** 2
Return sum / Tmp.Length
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
{
avg := NDArray.Sum / NDArray.Length
sum := 0
For i in NDArray
sum += (i - avg) ** 2
Return sum / NDArray.Length
}
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.Std(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
Static Sum(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
Return Tmp.Sum
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
Return NDArray.Sum
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.Sum(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
Static Tan(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
NDArray[i*] := Tan(NDArray[i*])
Return NDArray
}
Static Tanh(NDArray)
{
New_Pos := []
For i in NDArray.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
Tmp := NDArray[i*]
NDArray[i*] := (Exp(Tmp) - Exp(-Tmp)) / (Exp(Tmp) + Exp(-Tmp))
}
Return NDArray
}
Static Var(NDArray, Axis := "axis=None")
{
if InStr(Axis, "=")
Axis := Trim(StrSplit(Axis, "=")[2])
if Axis = "None"
{
Tmp := NDArray.Copy().Ravel()
avg := Tmp.Sum / Tmp.Length
sum := 0
For i in Tmp
sum += (i - avg) ** 2
Return Sqrt(sum / Tmp.Length)
}
Axis := (Axis == -1) ? NDArray.NDim - 1 : Axis
if NDArray.NDim == 1
{
avg := NDArray.Sum / NDArray.Length
sum := 0
For i in NDArray
sum += (i - avg) ** 2
Return Sqrt(sum / NDArray.Length)
}
else
{
Lst_Res := []
Lst_Pos := []
Tmp := NDArray.Shape
Loop NDArray.NDim
{
if Axis + 1 !== A_Index
Lst_Pos.Push(Range(1, Tmp[A_Index]).Array())
}
Lst_Pos := Dynloop_Loop(Lst_Pos)
For i in Lst_Pos
{
Tmp2 := []
Loop Tmp[Axis + 1]
{
iTmp := i.Clone()
iTmp.InsertAt(Axis + 1, A_Index)
Tmp2.Push(NDArray[iTmp*])
}
Lst_Res.Push(Numahk.Var(Tmp2))
}
Lst_Res := Numahk.NDArray(Lst_Res)
Tmp.RemoveAt(Axis + 1)
Lst_Res.Reshape(Tmp)
Return Lst_Res
}
}
; N-D Calculate
Static Add(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 + NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] + NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] + NDArray2
else
New_NDArray[i*] := NDArray1 + NDArray2[i*]
}
Return New_NDArray
}
Static Cmp(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return ListCmp(NDArray1, NDArray2)
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := ListCmp(NDArray1[i*], NDArray2[i*])
else if NDArray1.NDim
New_NDArray[i*] := ListCmp(NDArray1[i*], NDArray2)
else
New_NDArray[i*] := ListCmp(NDArray1, NDArray2[i*])
}
Return New_NDArray
}
Static CopySign(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
New_Pos := []
For i in NDArray1.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray2.NDim
{
Tmp := NDArray1[i*]
NDArray1[i*] := (Tmp.Sign == NDArray2[i*].Sign) ? Tmp : -Tmp
}
else
{
Tmp := NDArray1[i*]
NDArray1[i*] := (Tmp.Sign == NDArray2.Sign) ? Tmp : -Tmp
}
}
Return NDArray1
}
Static Divide(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 / NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] / NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] / NDArray2
else
New_NDArray[i*] := NDArray1 / NDArray2[i*]
}
Return New_NDArray
}
Static Equal(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 = NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] = NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] = NDArray2
else
New_NDArray[i*] := NDArray1 = NDArray2[i*]
}
Return New_NDArray
}
Static Floor_Divide(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return Floor(NDArray1 / NDArray2)
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := Floor(NDArray1[i*] / NDArray2[i*])
else if NDArray1.NDim
New_NDArray[i*] := Floor(NDArray1[i*] / NDArray2)
else
New_NDArray[i*] := Floor(NDArray1 / NDArray2[i*])
}
Return New_NDArray
}
Static FMax(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return Max(NDArray1, NDArray2)
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := Max(NDArray1[i*], NDArray2[i*])
else if NDArray1.NDim
New_NDArray[i*] := Max(NDArray1[i*], NDArray2)
else
New_NDArray[i*] := Max(NDArray1, NDArray2[i*])
}
Return New_NDArray
}
Static FMin(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return Min(NDArray1, NDArray2)
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := Min(NDArray1[i*], NDArray2[i*])
else if NDArray1.NDim
New_NDArray[i*] := Min(NDArray1[i*], NDArray2)
else
New_NDArray[i*] := Min(NDArray1, NDArray2[i*])
}
Return New_NDArray
}
Static Greater(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 > NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] > NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] > NDArray2
else
New_NDArray[i*] := NDArray1 > NDArray2[i*]
}
Return New_NDArray
}
Static Greater_Equal(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 >= NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] >= NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] >= NDArray2
else
New_NDArray[i*] := NDArray1 >= NDArray2[i*]
}
Return New_NDArray
}
Static Ix_(NDArray1, NDArray2)
{
Return Dynloop_Loop([NDArray1, NDArray2]).Reshape([NDArray1.Length, NDArray2.Length, 2])
}
Static Less(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 < NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] < NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] < NDArray2
else
New_NDArray[i*] := NDArray1 < NDArray2[i*]
}
Return New_NDArray
}
Static Less_Equal(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 <= NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] <= NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] <= NDArray2
else
New_NDArray[i*] := NDArray1 <= NDArray2[i*]
}
Return New_NDArray
}
Static Logical_And(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 & NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] & NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] & NDArray2
else
New_NDArray[i*] := NDArray1 & NDArray2[i*]
}
Return New_NDArray
}
Static Logical_Or(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 | NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] | NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] | NDArray2
else
New_NDArray[i*] := NDArray1 | NDArray2[i*]
}
Return New_NDArray
}
Static Logical_XOr(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 ^ NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] ^ NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] ^ NDArray2
else
New_NDArray[i*] := NDArray1 ^ NDArray2[i*]
}
Return New_NDArray
}
Static Maximum(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return Max(NDArray1, NDArray2)
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := Max(NDArray1[i*], NDArray2[i*])
else if NDArray1.NDim
New_NDArray[i*] := Max(NDArray1[i*], NDArray2)
else
New_NDArray[i*] := Max(NDArray1, NDArray2[i*])
}
Return New_NDArray
}
Static Minimum(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return Min(NDArray1, NDArray2)
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := Min(NDArray1[i*], NDArray2[i*])
else if NDArray1.NDim
New_NDArray[i*] := Min(NDArray1[i*], NDArray2)
else
New_NDArray[i*] := Min(NDArray1, NDArray2[i*])
}
Return New_NDArray
}
Static Mod(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return Mod(NDArray1, NDArray2)
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := Mod(NDArray1[i*], NDArray2[i*])
else if NDArray1.NDim
New_NDArray[i*] := Mod(NDArray1[i*], NDArray2)
else
New_NDArray[i*] := Mod(NDArray1, NDArray2[i*])
}
Return New_NDArray
}
Static Multiply(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 * NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] * NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] * NDArray2
else
New_NDArray[i*] := NDArray1 * NDArray2[i*]
}
Return New_NDArray
}
Static Not_Equal(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 != NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] != NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] != NDArray2
else
New_NDArray[i*] := NDArray1 != NDArray2[i*]
}
Return New_NDArray
}
Static Subtract(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 - NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] - NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] - NDArray2
else
New_NDArray[i*] := NDArray1 - NDArray2[i*]
}
Return New_NDArray
}
Static Power(NDArray1, NDArray2)
{
NDArray1 := (NDArray1.NDim) ? Numahk.NDArray(NDArray1) : NDArray1
NDArray2 := (NDArray2.NDim) ? Numahk.NDArray(NDArray2) : NDArray2
New_Pos := []
if !NDArray1.NDim && !NDArray2.NDim
Return NDArray1 ** NDArray2
NewShape := NDArray1.NDim ? NDArray1 : NDArray2
New_NDArray := Numahk.Zeros_Like(NewShape)
For i in NewShape.Shape
New_Pos.Push(Range(1, i).Array())
New_Pos := Dynloop_Loop(New_Pos)
For i in New_Pos
{
if NDArray1.NDim && NDArray2.NDim
New_NDArray[i*] := NDArray1[i*] ** NDArray2[i*]
else if NDArray1.NDim
New_NDArray[i*] := NDArray1[i*] ** NDArray2
else
New_NDArray[i*] := NDArray1 ** NDArray2[i*]
}
Return New_NDArray
}
; File
Static Load(String)
{
Return Numahk.NDArray(XA.XA_Load(String ".nahk"))
}
Static Loadz(String)
{
NDArrayz := []
TmpArray := XA.XA_Load(String ".nahk")
For i in TmpArray
NDArrayz.Push(Numahk.NDArray(i))
Return NDArrayz
}
Static Save(String, NDArray)
{
NDArrayFile := FileOpen(String ".nahk", "w")
NDArrayFile.Write(XA.XA_Save(NDArray))
NDArrayFile.Close()
}
Static Savez(String, NDArray*)
{
NDArrayFile := FileOpen(String ".nahk", "w")
NDArrayFile.Write(XA.XA_Save(NDArray))
NDArrayFile.Close()
}
; Random
Class Random
{
}
; Linalg
Class Linalg
{
}
}
Class Range
{
__New(Start, Stop := "", Step := 1)
{
if !Step
Throw ValueError("Step cannot be 0.")
if Stop == ""
{
Tmp := Start
Start := 1
Stop := Tmp
}
this.Length := Max(0, (Stop - Start) // Step + 1)
this.Start := Start
this.Step := Step
}
__Enum(_)
{
this.LoopTimes := this.Length
this.Index := 0
Return Fn
Fn(&Idx := 0, &Value := 0)
{
if !this.LoopTimes
Return False
Idx := (IsSet(Value)) ? this.Start + this.Index * this.Step : this.Index + 1
Value := this.Start + this.Index * this.Step
this.Index++
Return this.LoopTimes--
}
}
Array()
{
Ret := []
Loop this.Length
Ret.Push(this.Start + (A_Index - 1) * this.Step)
Return Ret
}
}
Dynloop_Loop(Data)
{
if Data.Length == 1
{
Lst_rst := []
For i in Data[1]
Lst_rst.Push([i])
Return Lst_rst
}
Max_y_idx := Data.Length
Row_max_idx := 1
Arr_len := []
Lst_row := []
Lst_rst := []
Arr_idx := [0]
Arr_idx.mul(Max_y_idx)
For item in Data
{
_n := item.Length
Arr_len.Push(_n)
Lst_row.Extend(item)
Row_max_idx *= _n
}
For row_idx in Range(0, Row_max_idx - 1)
{
For y_idx in Range(0, Max_y_idx - 1)
{
_pdt := 1
For n in Range(y_idx + 1, Arr_len.Length - 1)
_pdt *= Arr_len[n + 1]
_offset := 0
For n in Range(0, y_idx - 1)
_offset += Arr_len[n + 1]
Arr_idx[y_idx + 1] := Mod((row_idx // _pdt), Arr_len[y_idx + 1]) + _offset
}
_lst_tmp := []
For idx in Arr_idx
_lst_tmp.Push(Lst_row[idx + 1])
Lst_rst.Push(_lst_tmp)
}
Return Lst_rst
}二、重置原因与使用说明
我重置Numahk主要是因为挂载一个自定义的类型List，有两点原因。一个是不便于初始化，对普通用户不够友好；另一个是索引混乱不符合AHK自有逻辑。
所以我进行了重置，这里说明Numahk的三个特点：
- 存取访问速度较前版本快了3-5倍。
- 下标索引统一从1开始，符合AHK用户逻辑。
- 优化所有函数以及源码呈现形式，在源码中按功能进行分块，便于快速定位源码。
三、函数说明
见参考文档，后续会整理。注意，目前Random和Linalg模块暂未开放，敬请期待。
???