0
1

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

More than 5 years have passed since last update.

pySpherepts > sphereptsのNumpy、Scipy実装 > 結果: sphereptsと同じものが得られた

Last updated at Posted at 2017-11-25
動作環境
GeForce GTX 1070 (8GB)
ASRock Z170M Pro4S [Intel Z170chipset]
Ubuntu 16.04 LTS desktop amd64
TensorFlow v1.2.1
cuDNN v5.1 for Linux
CUDA v8.0
Python 3.5.2
IPython 6.0.0 -- An enhanced Interactive Python.
gcc (Ubuntu 5.4.0-6ubuntu1~16.04.4) 5.4.0 20160609
GNU bash, version 4.3.48(1)-release (x86_64-pc-linux-gnu)
scipy v0.19.1
geopandas v0.3.0
MATLAB R2017b (Home Edition)

関連

pySpherepts

上記のMATLABコードの一部であるgetIcosNodes.mをNumpy+Scipy実装した。
名称は「pySpherepts」とする。

コードは中身を整理した上でGitHubに上げる予定。

pySpherepts結果

使ってみた。

test_getIcosNode_171126.py
import numpy as np
import getIcosNodes_171126 as gIN

xs, tris = gIN.getIcosNodes(4,0)
print(xs)
np.savetxt('res_IcodsNodes_4_0_171126.txt', xs, fmt="%.5f")

上記の結果をMATLAB上で読込み、オリジナルのsphereptsの表示プログラムplotSphNodesで表示してみた。

test_numpyIcosNodes_171126.m
function test_numpyIcosNodes_171126()
    xs = load('res_IcodsNodes_4_0_171126.txt')
    plotSphNodes(xs)

qiita.png

spherepts結果

オリジナルのsphereptsの結果

>> [x0,tri0] = getIcosNodes(4,0);
>> plotSphNodes(x0)

qiita.png

同じものが得られた。

MATLAB分からなかったので、勉強しながら20時間はかかっただろうか。

0
1
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
0
1

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?