swift で一次元、二次元配列を取り扱うクラスを作る
swiftで一次元、二次元配列を操作するクラスを作成しKMeansで有用性を確認しました
ソースはこちら
配列クラス: https://github.com/ikuo0/swift_number/blob/main/numxd.swift
KMeans実行部: https://github.com/ikuo0/swift_number/blob/main/main.swift
コマンド実行する場合は各々同じディレクトリに設置し
swiftc -emit-executable main.swift numxd.swift
とします、コマンドを実行すると main という名前のファイルが出来ますのでそれを実行してください
KMeans実行部
KMeans部分のソースです、70行ほどでKMeansが実行できていますのでだいぶ簡単に配列操作ができていると思います。
funcTest2D(){do{letn1dI=Num1D<Int>()letn1dD=Num1D<Double>()letn2d=Num2D<Double>()lettsvData=tryTSV.Read("seeds_dataset.txt")letdata=tryTSV.ToDouble(tsvData)letxIndexes=n1dI.Arange(0,7)letxData=data.IndexingT(xIndexes)letyData=data.LineT(7)////////////////////////////////////////// scaling////////////////////////////////////////letscaleMean=xData.Mean()letscaleVriance=xData.Variance()letscaleStddev=scaleVriance.Sqrt()letscaledX=(xData-scaleMean)/scaleStddev////////////////////////////////////////// KMeans////////////////////////////////////////letClusters=3letmaxIter=100lettol=1e-5// InitializeletrowIndexes=n1dI.Arange(0,yData.Count)letshuffled=rowIndexes.Shuffle()letinitIndexes=shuffled.Slice(0,Clusters)letinitMeans=scaledX.Indexing(initIndexes)varmeans=initMeansletpredict=n1dI.Create(xData.Row)foriterin0..<maxIter{// EStepformin0..<xData.Row{letdistances=n1dD.Create(Clusters)forclusterin0..<Clusters{letsubtract=scaledX.Line(m)-means.Line(cluster)letpower=subtract.Power(2)lettotal=power.Total()letdistance=sqrt(total)distances[cluster]=distance}predict[m]=distances.ArgMin()}// MStepletnewMeans=n2d.Create(Clusters,scaledX.Col)forclusterin0..<Clusters{letcIndexes=predict.WhereEq(cluster)letcMean=scaledX.Indexing(cIndexes)newMeans[cluster]=cMean.Mean().Value}// thresholdletsubtract=means-newMeansletpower=subtract.Power(2)lettotal=power.TotalX()letmeansDistance=sqrt(total/Double(Clusters))means=newMeansprint("iter=\(iter), means distance=\(meansDistance)")ifmeansDistance<tol{break}}Dump2D(means)Dump1D(predict)}catch{print("exception",error)}}
わざわざクラス化する必要があったかどうか
C/C++と異なりメモリ確保の面倒さが無いためわざわざクラス化する事も無かったんじゃないかという気もしましたが、作ってみたら思いの外短く書けているので有りかなと思いました
ジェネリクス関数の書き込みの仕組み
https://github.com/ikuo0/swift_number/blob/main/numxd.swift#L27
自作の WriteGVariable という関数でジェネリクス関数内の代入処理を行っています
WriteGVariable関数を使用しない場合に型を確定させてからキャストをするコードをその都度書かなければいけません
ifT.selfisDouble.Type{dst[i]=Double(r)as!T}elseifT.selfisInt.Type{dst[i]=Int(r)as!T}// 冗長になってしまう
WriteGVariable関数を使用することで1行で代入処理を書くことができます
WriteGVariable(&dst[i],r)// 1行で済む
swiftが厳格な型付け言語であるためにジェネリクス型を使用した場合に冗長になるか冗長回避のため対策が必要になるのは仕方のないことだと思います。
ジェネリクス型に限らず型違い同士の代入処理等については下記URLにて答えが出ているように思います
https://academy.realm.io/jp/posts/richard-fox-casting-swift-2/
https://qiita.com/ObuchiYuki/items/a945efd14d3a05a19f75
答えが出ているにも関わらず真似しないのはoperator演算子のようなソースを隠蔽してしまう仕様が好みではないのと私のswift理解度が低いからです、理解しきれていない物をコピーするのは気が引けますので。
swiftの勉強のためC++で既に実装した物をswiftで実装し直しました、色々と発見があったりしておもしろかったです。
以上です。