Skip to content

Instantly share code, notes, and snippets.

View DanHickstein's full-sized avatar

Danhickstein DanHickstein

  • Octave Photonics
  • Boulder, CO
View GitHub Profile
clear all;
N=1001;
Rm = fix(N/2)+1;
NBF=500;
I=(1:N)';
R2 = ((I-Rm).^2);
R = I-Rm;
M = zeros(N,NBF);
Mc =zeros(N,NBF);
Mc(:,1)=exp(-R2);
#include "stdosl.h"
shader volumeTexture(
point Vector = P,
string directory="//",
string prefix="",
int framePositions=3,
string extension=".tif",
float zMax = 2.0,
int fromFrame = 1,
{
"action": "published",
"organization": {
"issues_url": "https://api.github.com/orgs/PyAbel/issues",
"members_url": "https://api.github.com/orgs/PyAbel/members{/member}",
"description": "A Python package for performing inverse abel transforms",
"public_members_url": "https://api.github.com/orgs/PyAbel/public_members{/member}",
"url": "https://api.github.com/orgs/PyAbel",
"events_url": "https://api.github.com/orgs/PyAbel/events",
"avatar_url": "https://avatars3.githubusercontent.com/u/15894666?v=4",
@DanHickstein
DanHickstein / comb.py
Created April 11, 2019 14:16
comb comparison
import abel
from abel.tools.analytical import PiecewisePolynomial
from itertools import chain
import matplotlib.pyplot as plt
import numpy as np
hw = 1 # peak half-width
step = 10 # center-to-center distance between peaks
n = 10 # number of peaks
@DanHickstein
DanHickstein / comb.py
Created April 11, 2019 14:16
comb comparison
import abel
from abel.tools.analytical import PiecewisePolynomial
from itertools import chain
import matplotlib.pyplot as plt
import numpy as np
hw = 1 # peak half-width
step = 10 # center-to-center distance between peaks
n = 10 # number of peaks
import abel
from abel.tools.analytical import PiecewisePolynomial
from itertools import chain
import matplotlib.pyplot as plt
import numpy as np
hw = 1 # peak half-width - Or is this the full-width?
step = 10 # center-to-center distance between peaks
n = 10 # number of peaks
import abel
from abel.tools.analytical import PiecewisePolynomial
from itertools import chain
import matplotlib.pyplot as plt
import numpy as np
hw = 1 # peak half-width - Or is this the full-width?
step = 10 # center-to-center distance between peaks
n = 10 # number of peaks
import numpy as np
import matplotlib.pyplot as plt
import sys, os
import abel
import scipy.ndimage
import bz2
transforms = [
("basex" ,abel.basex.basex_transform),
("direct" ,abel.direct.direct_transform),
@DanHickstein
DanHickstein / basex(var).dat
Created April 26, 2019 00:38
i7-8700 benchmarks
# basex(var) iabel(ms)
5 0.52261
9 0.49675
17 0.50111
33 0.54974
65 0.69175
129 0.81117
257 1.53354
513 4.21259
1025 24.40259
import numpy as np
import matplotlib.pyplot as plt
from pynlo.interactions.FourWaveMixing import SSFM
from pynlo.media.fibers import fiber
from pynlo.light.DerivedPulses import SechPulse
########## Parameters:
Pulse1 = .200 # pulse duration (ps)
Pulse2 = .200
Pulse3 = .200