This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
[etalon:.+Software/Python]$ git clone https://github.com/PyAbel/PyAbel.git | |
Cloning into 'PyAbel'... | |
remote: Enumerating objects: 7858, done. | |
remote: Counting objects: 100% (675/675), done. | |
remote: Compressing objects: 100% (372/372), done. | |
remote: Total 7858 (delta 401), reused 508 (delta 285), pack-reused 7183 | |
Receiving objects: 100% (7858/7858), 22.54 MiB | 9.50 MiB/s, done. | |
Resolving deltas: 100% (5611/5611), done. | |
[etalon:.+Software/Python]$ mv PyAbel PyAbel-master | |
[etalon:.+Software/Python]$ python -V |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
import cse | |
import matplotlib.pyplot as plt | |
from scipy.optimize import curve_fit | |
import scipy.constants as const | |
from scipy.special import spherical_jn, spherical_yn | |
au = 0.52917720859 | |
eh = 27.21138386 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Name: Intel Core i7-6700 | |
uArch: Skylake | |
Technology: 14nm | |
Max Freq: 4.000 GHz | |
Cores: 4 cores (8 threads) | |
AVX: AVX,AVX2 | |
FMA: FMA3 | |
L1i Size: 32KB (128KB Total) | |
L1d Size: 32KB (128KB Total) | |
L2 Size: 256KB (1MB Total) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
import abel | |
import matplotlib.pyplot as plt | |
from scipy.ndimage.interpolation import shift | |
fig, ax = plt.subplots(2, 2, figsize=(8, 8)) | |
n = 201 | |
n2 = n//2 | |
geo_origin = (n2, n2) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
01-28 19:28:19.603 18040 18040 W tech.ula: JIT profile information will not be recorded: profile file does not exits. | |
01-28 19:28:19.606 18040 18040 I Perf : Connecting to perf service. | |
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to ensure /data/user/0/tech.ula/cache: mkdir failed: EACCES (Permission denied) | |
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to update user.inode_cache: stat failed: EACCES (Permission denied) | |
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to ensure /data/user_de/0/tech.ula/code_cache: mkdir failed: ENOENT (No such file or directory) | |
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to update user.inode_code_cache: stat failed: ENOENT (No such file or directory) | |
01-28 19:28:19.619 602 602 I hwservicemanager: getTransport: Cannot find entry vendor.qti.hardware.iop@2.0::IIop/default in either framework or device manifest. | |
01-28 19:28:19.619 971 971 E ANDR-IOP: IIop:: Iop HAL Service is not available. | |
01-28 19:28:19.620 18040 18059 E Perf : Fail t |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
import matplotlib.pyplot as plt | |
from pbasex import pbasex, loadG | |
from quadrant import foldQuadrant, resizeFolded | |
import abel | |
# Load images | |
im = np.loadtxt('O2-ANU1024.txt.bz2') | |
raw = abel.tools.center.center_image(im, square=True) | |
n = raw.shape[-1] |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
import abel | |
import bz2 | |
import anisotro | |
from scipy.signal import find_peaks, peak_widths | |
import matplotlib.pylab as plt | |
# O2- PES VM-image | |
imagefile = bz2.BZ2File('O2-ANU1024.txt.bz2') |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
# basex(var) iabel(ms) | |
5 0.53462 | |
9 0.53325 | |
17 0.54541 | |
33 0.56431 | |
65 0.67819 | |
129 0.80625 | |
257 1.45382 | |
513 3.34333 | |
1025 14.87905 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import abel | |
import matplotlib.pyplot as plt | |
im = abel.tools.analytical.SampleImage(1001, 'Ominus').image | |
proj = abel.Transform(im, method='basex', verbose=False, | |
direction='forward').transform | |
bas = abel.Transform(proj, method='basex', verbose=False, | |
angular_integration=True) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
import abel | |
import matplotlib.pyplot as plt | |
n = 101 | |
fig, axs = plt.subplots(3, 2, sharex=True, sharey=True) | |
case = '1st-order-diff' | |
# case = 'gradient' |
NewerOlder