spherepts Numpy > np.concatenate()と自前のarray結合の処理時間比較 > np.concatenate()が3割ほど速いperformancesphereptsdifferencenumpynumpy spherepts > Delaunay > ScipyのDelaunayに置き換えても同じ結果になる (適用する問題によるのだろう)Delaunaysphereptsscipydifferencescipy MATLAB > API > Delaunay | Voronoi 試してみたVoronoigeometrysphereptsmatlabDelaunaymatlab MATLAB > freeBoundary() > 四面体結合物への適用 > 隠面以外の三角形の情報を取得するgeometrysphereptsmatlabapimatlab MATLAB > 四面体メッシュプロット > tetramesh() | view(30,30) | サブセット表示 tetramesh(tet(2:3,:),X)sphereptsmatlabapimatlab onlineTool > MATLAB > Execute MATLAB/Octave Online | ideoneonlineToolmatlabsphereptsmatlab MATLAB > load tetmesh > tetとX > nodeとelementmatlabsphereptskeywordmatlab Ubuntu 16.04 > MATLABセットアップ > Home Editionmatlabsetupsphereptsmatlab MATLAB > tetmesh > tet: mesh, X: element > X:座標情報, tet: 四角形の座標インデックス?matlabsphereptsmatlab MATLAB > freeBoundary() > 四角形をfreeBoundaryしてみる > 2つでなく4つの三角形のインデックスが得られるapimatlabsphereptsmatlab geometry + Python > 四面体結合体 > 隠面の「頂点インデックス組合せ」を見つける > v0.1geometrymatlabsphereptsmatlab MATLAB > array > 行をswapするmatlabsphereptslistOperationmatlab MATLAB > .mファイルを作って実行する > 関数の定義と利用apimatlabsphereptsmatlab MATLAB > random > 整数の乱数を得る > X = round(rand()*100)matlabsphereptsRandommatlab MATLAB > Error: Too many output arguments. > 関数の戻り値を指定することで対処するerrorapimatlabsphereptsmatlab MATLAB > NxM array > 比較 > y1 == y2 | 行列なのかarrayなのかmatlabsphereptsPrimermatlab MATLAB > NxM array > 比較 > y1とy2が異なるかを確認する > unique(y1 == y2) にて 0がある(異なる) > isequal() , ifでtrueが返される条件, all()matlabsphereptsmatlab MATLAB > COMMAND WINDOWの出力をファイル保存 > diarymatlabcommandsphereptsredirectmatlab MATLAB > 変数の内容をファイルに保存する > save('pi.txt', 'X', '-ascii')apimatlabsphereptsfileIOmatlab MATLAB > IDE > ショートカットのキー追加 > コピーのCtrl+c割当 | 貼り付けのCtrl+vは割当できないmatlabsetupshortcutsphereptsmatlab Numpy > fileIO > MATLAB saveコマンドによるasciiファイルを読込む > genfromtxt()使用numpymatlabsphereptsfileIOmatlab MATLAB > Delaunay > T = delaunay(x, y, {'Qt','Qbb', 'Qc', 'Qz'}) > Error: DELAUNAY では、Qhull 固有のオプションはサポートされなくなったか、不要になりました。errormatlabsphereptsDelaunaymatlab MATLAB > fileIO > load() > save()の結果と同じ書式でなくても読込むことができるmatlabsphereptsfileIOmatlab MATLAB > ./ > rdivideを使うと良い > ソースリーディングの観点から | リンク: 省略表記 | ドキュメント和訳matlabsphereptsdifficultToSearchmatlab MATLAB > rdivide > right array divisionとは | Right-array division > rdivide(A,B) > Bによる除算 | 二つのAPIを用意する欠点と利点apimatlabsphereptskeywordmatlab MATLAB > rdivide > arrayと一変数での右除算 | arrayと列ベクトルでの右除算 | sum(A.^2, 2)apimatlabsphereptsmatlab MATLAB: repmat(inpx, [1 3]) > Numpy: np.tile(inpx, (1,3))apimatlabsphereptsmatlab MATLAB: acos(44) > 0.0000 + 4.4772i | Numpy: np.arccos(44) > nannumpymatlabsphereptsmatlab MATLAB: dot(x1, x2, 2) > dot product > Numpy: 実装 v0.1, v0.2 > 一次元と二次元は確認numpymatlabsphereptsmatlab MATLAB > unique(vs, 'rows')を使ってみた > [res, ~, idx] = unique(vs, 'rows')matlabsphereptsmatlab MATLAB: unique(vs, 'rows') > Numpy: 実装 > v0.1,..v0.3 | np.where(unqs == alist)は符号に関係なくインデックスを返すようだnumpymatlabsphereptstrapmatlab MATLAB:fix() ゼロ方向への丸め > Numpy:実装 > 正: floor(), 負: ceil()numpymatlabsphereptsmatlab MATLAB: fliplr() > Numpy:実装 > 1次元とN(>1)次元で別処理が必要 | np.flipud(), np.fliplr(), [::-1], ascontiguousarraynumpymatlabsphereptsmatlab Numpy > MATLABのsph2cart()の実装 > sph2cart_180204.py > 直交座標から球座標への変換 v0.1..v0.3sphereptsnumpygeometrymatlabpySphereptsmatlab
Numpy > np.concatenate()と自前のarray結合の処理時間比較 > np.concatenate()が3割ほど速いperformancesphereptsdifferencenumpynumpy spherepts > Delaunay > ScipyのDelaunayに置き換えても同じ結果になる (適用する問題によるのだろう)Delaunaysphereptsscipydifferencescipy MATLAB > API > Delaunay | Voronoi 試してみたVoronoigeometrysphereptsmatlabDelaunaymatlab MATLAB > freeBoundary() > 四面体結合物への適用 > 隠面以外の三角形の情報を取得するgeometrysphereptsmatlabapimatlab MATLAB > 四面体メッシュプロット > tetramesh() | view(30,30) | サブセット表示 tetramesh(tet(2:3,:),X)sphereptsmatlabapimatlab onlineTool > MATLAB > Execute MATLAB/Octave Online | ideoneonlineToolmatlabsphereptsmatlab MATLAB > load tetmesh > tetとX > nodeとelementmatlabsphereptskeywordmatlab Ubuntu 16.04 > MATLABセットアップ > Home Editionmatlabsetupsphereptsmatlab MATLAB > tetmesh > tet: mesh, X: element > X:座標情報, tet: 四角形の座標インデックス?matlabsphereptsmatlab MATLAB > freeBoundary() > 四角形をfreeBoundaryしてみる > 2つでなく4つの三角形のインデックスが得られるapimatlabsphereptsmatlab geometry + Python > 四面体結合体 > 隠面の「頂点インデックス組合せ」を見つける > v0.1geometrymatlabsphereptsmatlab MATLAB > array > 行をswapするmatlabsphereptslistOperationmatlab MATLAB > .mファイルを作って実行する > 関数の定義と利用apimatlabsphereptsmatlab MATLAB > random > 整数の乱数を得る > X = round(rand()*100)matlabsphereptsRandommatlab MATLAB > Error: Too many output arguments. > 関数の戻り値を指定することで対処するerrorapimatlabsphereptsmatlab MATLAB > NxM array > 比較 > y1 == y2 | 行列なのかarrayなのかmatlabsphereptsPrimermatlab MATLAB > NxM array > 比較 > y1とy2が異なるかを確認する > unique(y1 == y2) にて 0がある(異なる) > isequal() , ifでtrueが返される条件, all()matlabsphereptsmatlab MATLAB > COMMAND WINDOWの出力をファイル保存 > diarymatlabcommandsphereptsredirectmatlab MATLAB > 変数の内容をファイルに保存する > save('pi.txt', 'X', '-ascii')apimatlabsphereptsfileIOmatlab MATLAB > IDE > ショートカットのキー追加 > コピーのCtrl+c割当 | 貼り付けのCtrl+vは割当できないmatlabsetupshortcutsphereptsmatlab Numpy > fileIO > MATLAB saveコマンドによるasciiファイルを読込む > genfromtxt()使用numpymatlabsphereptsfileIOmatlab MATLAB > Delaunay > T = delaunay(x, y, {'Qt','Qbb', 'Qc', 'Qz'}) > Error: DELAUNAY では、Qhull 固有のオプションはサポートされなくなったか、不要になりました。errormatlabsphereptsDelaunaymatlab MATLAB > fileIO > load() > save()の結果と同じ書式でなくても読込むことができるmatlabsphereptsfileIOmatlab MATLAB > ./ > rdivideを使うと良い > ソースリーディングの観点から | リンク: 省略表記 | ドキュメント和訳matlabsphereptsdifficultToSearchmatlab MATLAB > rdivide > right array divisionとは | Right-array division > rdivide(A,B) > Bによる除算 | 二つのAPIを用意する欠点と利点apimatlabsphereptskeywordmatlab MATLAB > rdivide > arrayと一変数での右除算 | arrayと列ベクトルでの右除算 | sum(A.^2, 2)apimatlabsphereptsmatlab MATLAB: repmat(inpx, [1 3]) > Numpy: np.tile(inpx, (1,3))apimatlabsphereptsmatlab MATLAB: acos(44) > 0.0000 + 4.4772i | Numpy: np.arccos(44) > nannumpymatlabsphereptsmatlab MATLAB: dot(x1, x2, 2) > dot product > Numpy: 実装 v0.1, v0.2 > 一次元と二次元は確認numpymatlabsphereptsmatlab MATLAB > unique(vs, 'rows')を使ってみた > [res, ~, idx] = unique(vs, 'rows')matlabsphereptsmatlab MATLAB: unique(vs, 'rows') > Numpy: 実装 > v0.1,..v0.3 | np.where(unqs == alist)は符号に関係なくインデックスを返すようだnumpymatlabsphereptstrapmatlab MATLAB:fix() ゼロ方向への丸め > Numpy:実装 > 正: floor(), 負: ceil()numpymatlabsphereptsmatlab MATLAB: fliplr() > Numpy:実装 > 1次元とN(>1)次元で別処理が必要 | np.flipud(), np.fliplr(), [::-1], ascontiguousarraynumpymatlabsphereptsmatlab Numpy > MATLABのsph2cart()の実装 > sph2cart_180204.py > 直交座標から球座標への変換 v0.1..v0.3sphereptsnumpygeometrymatlabpySphereptsmatlab