Unify XSpectrum types
#101
icweaverAnnotations
10 errors
|
build:
src/Spectra.jl#L194
doctest failure in src/Spectra.jl:194-208
```jldoctest
julia> using Unitful, UnitfulAstro, Measurements
julia> wave = range(1, 4, length=1000)u"μm";
julia> sigma = randn(size(wave));
julia> flux = (100 .± sigma)u"erg/cm^2/s/angstrom";
julia> spec = spectrum(wave, flux)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 4.0 μm)
flux: (100.0 ± 1.2 erg Å^-1 cm^-2 s^-1, 100.0 ± 1.1 erg Å^-1 cm^-2 s^-1)
meta: Dict{Symbol, Any}()
```
Subexpression:
spec = spectrum(wave, flux)
Evaluated output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave ((1000,)): 1.0 μm .. 4.0 μm
flux ((1000,)): 100.0 ± 1.2 erg Å^-1 cm^-2 s^-1 .. 100.0 ± 1.1 erg Å^-1 cm^-2 s^-1
meta: Dict{Symbol, Any}()
Expected output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 4.0 μm)
flux: (100.0 ± 1.2 erg Å^-1 cm^-2 s^-1, 100.0 ± 1.1 erg Å^-1 cm^-2 s^-1)
meta: Dict{Symbol, Any}()
diff =
Warning: Diff output requires color.
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, wave ((1000,)): 1.0 μm .. 4.0 μm)
flux: (100.0 μm
flux ((1000,)): 100.0 ± 1.2 erg Å^-1 cm^-2 s^-1, s^-1 .. 100.0 ± 1.1 erg Å^-1 cm^-2 s^-1)
s^-1
meta: Dict{Symbol, Any}()
|
|
build:
src/Spectra.jl#L170
doctest failure in src/Spectra.jl:170-189
```jldoctest
julia> wave = range(1e4, 4e4, length=1000);
julia> flux = 100 .* ones(size(wave));
julia> spec = spectrum(wave, flux)
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
meta: Dict{Symbol, Any}()
julia> spec = spectrum(wave, flux, name="Just Noise")
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
meta: Dict{Symbol, Any}(:name => "Just Noise")
julia> spec.name
"Just Noise"
```
Subexpression:
spec = spectrum(wave, flux, name="Just Noise")
Evaluated output:
SingleSpectrum(Float64, Float64)
wave ((1000,)): 10000.0 .. 40000.0
flux ((1000,)): 100.0 .. 100.0
meta: Dict{Symbol, Any}(:name => "Just Noise")
Expected output:
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
meta: Dict{Symbol, Any}(:name => "Just Noise")
diff =
Warning: Diff output requires color.
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
wave ((1000,)): 10000.0 .. 40000.0
flux ((1000,)): 100.0 .. 100.0
meta: Dict{Symbol, Any}(:name => "Just Noise")
|
|
build:
src/Spectra.jl#L170
doctest failure in src/Spectra.jl:170-189
```jldoctest
julia> wave = range(1e4, 4e4, length=1000);
julia> flux = 100 .* ones(size(wave));
julia> spec = spectrum(wave, flux)
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
meta: Dict{Symbol, Any}()
julia> spec = spectrum(wave, flux, name="Just Noise")
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
meta: Dict{Symbol, Any}(:name => "Just Noise")
julia> spec.name
"Just Noise"
```
Subexpression:
spec = spectrum(wave, flux)
Evaluated output:
SingleSpectrum(Float64, Float64)
wave ((1000,)): 10000.0 .. 40000.0
flux ((1000,)): 100.0 .. 100.0
meta: Dict{Symbol, Any}()
Expected output:
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
meta: Dict{Symbol, Any}()
diff =
Warning: Diff output requires color.
SingleSpectrum(Float64, Float64)
wave: (10000.0, 40000.0)
flux: (100.0, 100.0)
wave ((1000,)): 10000.0 .. 40000.0
flux ((1000,)): 100.0 .. 100.0
meta: Dict{Symbol, Any}()
|
|
build:
src/utils.jl#L17
doctest failure in src/utils.jl:17-40
```jldoctest
julia> using Spectra, Unitful, UnitfulAstro
julia> wave = range(1, 3, length=100)u"μm"
(1.0:0.020202020202020204:3.0) μm
julia> bb = blackbody(wave, 2000u"K")
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (49010.54557924032 W μm^-1 m^-2, 131058.60552995963 W μm^-1 m^-2)
meta: Dict{Symbol, Any}(:T => 2000 K, :name => "Blackbody")
julia> blackbody(ustrip.(u"angstrom", wave), 6000)
SingleSpectrum(Float64, Float64)
wave: (10000.0, 30000.0)
flux: (40.04325690910415, 1190.9562575755397)
meta: Dict{Symbol, Any}(:T => 6000, :name => "Blackbody")
julia> bb.wave[argmax(bb)]
1.4444444444444444 μm
julia> 2898u"μm*K" / bb.T # See if it matches up with Wien's law
1.449 μm
```
Subexpression:
blackbody(ustrip.(u"angstrom", wave), 6000)
Evaluated output:
SingleSpectrum(Float64, Float64)
wave ((100,)): 10000.0 .. 30000.0
flux ((100,)): 1190.9562575755397 .. 40.04325690910415
meta: Dict{Symbol, Any}(:T => 6000, :name => "Blackbody")
Expected output:
SingleSpectrum(Float64, Float64)
wave: (10000.0, 30000.0)
flux: (40.04325690910415, 1190.9562575755397)
meta: Dict{Symbol, Any}(:T => 6000, :name => "Blackbody")
diff =
Warning: Diff output requires color.
SingleSpectrum(Float64, Float64)
wave: (10000.0, 30000.0)
flux: (40.04325690910415, 1190.9562575755397)
wave ((100,)): 10000.0 .. 30000.0
flux ((100,)): 1190.9562575755397 .. 40.04325690910415
meta: Dict{Symbol, Any}(:T => 6000, :name => "Blackbody")
|
|
build:
src/utils.jl#L17
doctest failure in src/utils.jl:17-40
```jldoctest
julia> using Spectra, Unitful, UnitfulAstro
julia> wave = range(1, 3, length=100)u"μm"
(1.0:0.020202020202020204:3.0) μm
julia> bb = blackbody(wave, 2000u"K")
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (49010.54557924032 W μm^-1 m^-2, 131058.60552995963 W μm^-1 m^-2)
meta: Dict{Symbol, Any}(:T => 2000 K, :name => "Blackbody")
julia> blackbody(ustrip.(u"angstrom", wave), 6000)
SingleSpectrum(Float64, Float64)
wave: (10000.0, 30000.0)
flux: (40.04325690910415, 1190.9562575755397)
meta: Dict{Symbol, Any}(:T => 6000, :name => "Blackbody")
julia> bb.wave[argmax(bb)]
1.4444444444444444 μm
julia> 2898u"μm*K" / bb.T # See if it matches up with Wien's law
1.449 μm
```
Subexpression:
bb = blackbody(wave, 2000u"K")
Evaluated output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave ((100,)): 1.0 μm .. 3.0 μm
flux ((100,)): 89534.30930426194 W μm^-1 m^-2 .. 49010.54557924032 W μm^-1 m^-2
meta: Dict{Symbol, Any}(:T => 2000 K, :name => "Blackbody")
Expected output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (49010.54557924032 W μm^-1 m^-2, 131058.60552995963 W μm^-1 m^-2)
meta: Dict{Symbol, Any}(:T => 2000 K, :name => "Blackbody")
diff =
Warning: Diff output requires color.
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, wave ((100,)): 1.0 μm .. 3.0 μm)
flux: (49010.54557924032 μm
flux ((100,)): 89534.30930426194 W μm^-1 m^-2, 131058.60552995963 m^-2 .. 49010.54557924032 W μm^-1 m^-2)
m^-2
meta: Dict{Symbol, Any}(:T => 2000 K, :name => "Blackbody")
|
|
build:
docs/src/index.md#L54
doctest failure in docs/src/index.md:54-62
```jldoctest guide
julia> cont_fit = continuum(spec)
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (3815.0483f0 Å, 9206.613f0 Å)
flux: (0.4582525451059156 erg Å^-1 cm^-2 s^-1, 4.402153019328562 erg Å^-1 cm^-2 s^-1)
meta: Dict{Symbol, Any}(:coeffs => Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}}[1.983152216046405e-15 erg Å^-1 cm^-2 s^-1, -1.8822245369267038e-16 erg Å^-1 cm^-2 s^-1, -1.0422750370065006e-16 erg Å^-1 cm^-2 s^-1, 4.8112282273206135e-17 erg Å^-1 cm^-2 s^-1], :normalized => true)
julia> plot(cont_fit, xlims=(6545, 6600));
```
Subexpression:
cont_fit = continuum(spec)
Evaluated output:
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave ((3827,)): 3815.0483f0 Å .. 9206.613f0 Å
flux ((3827,)): 1.0808438837160355 erg Å^-1 cm^-2 s^-1 .. 1.0098373106940344 erg Å^-1 cm^-2 s^-1
meta: Dict{Symbol, Any}(:coeffs => Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}}[1.983152216046405e-15 erg Å^-1 cm^-2 s^-1, -1.8822245369267038e-16 erg Å^-1 cm^-2 s^-1, -1.0422750370065006e-16 erg Å^-1 cm^-2 s^-1, 4.8112282273206135e-17 erg Å^-1 cm^-2 s^-1], :normalized => true)
Expected output:
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (3815.0483f0 Å, 9206.613f0 Å)
flux: (0.4582525451059156 erg Å^-1 cm^-2 s^-1, 4.402153019328562 erg Å^-1 cm^-2 s^-1)
meta: Dict{Symbol, Any}(:coeffs => Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}}[1.983152216046405e-15 erg Å^-1 cm^-2 s^-1, -1.8822245369267038e-16 erg Å^-1 cm^-2 s^-1, -1.0422750370065006e-16 erg Å^-1 cm^-2 s^-1, 4.8112282273206135e-17 erg Å^-1 cm^-2 s^-1], :normalized => true)
diff =
Warning: Diff output requires color.
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (3815.0483f0 Å, wave ((3827,)): 3815.0483f0 Å .. 9206.613f0 Å)
flux: (0.4582525451059156 Å
flux ((3827,)): 1.0808438837160355 erg Å^-1 cm^-2 s^-1, 4.402153019328562 s^-1 .. 1.0098373106940344 erg Å^-1 cm^-2 s^-1)
s^-1
meta: Dict{Symbol, Any}(:coeffs => Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}}[1.983152216046405e-15 erg Å^-1 cm^-2 s^-1, -1.8822245369267038e-16 erg Å^-1 cm^-2 s^-1, -1.0422750370065006e-16 erg Å^-1 cm^-2 s^-1, 4.8112282273206135e-17 erg Å^-1 cm^-2 s^-1], :normalized => true)
|
|
build:
docs/src/index.md#L23
doctest failure in docs/src/index.md:23-50
```jldoctest guide
julia> using Spectra, FITSIO, Unitful, UnitfulAstro, Plots
julia> fitsurl = "https://dr14.sdss.org/optical/spectrum/view/data/format=fits/spec=lite?plateid=1323&mjd=52797&fiberid=12";
julia> # download(fitsurl, "sdss.fits");
julia> f = FITS("sdss.fits")
File: sdss.fits
Mode: "r" (read-only)
HDUs: Num Name Type
1 Image
2 COADD Table
3 SPECOBJ Table
4 SPZLINE Table
julia> wave = (10 .^ read(f[2], "loglam"))u"angstrom";
julia> flux = (read(f[2], "flux") .* 1e-17)u"erg/s/cm^2/angstrom";
julia> spec = spectrum(wave, flux)
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (3815.0483f0 Å, 9206.613f0 Å)
flux: (9.416705322265625e-16 erg Å^-1 cm^-2 s^-1, 9.156107177734375e-15 erg Å^-1 cm^-2 s^-1)
meta: Dict{Symbol, Any}()
julia> plot(spec);
```
Subexpression:
spec = spectrum(wave, flux)
Evaluated output:
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave ((3827,)): 3815.0483f0 Å .. 9206.613f0 Å
flux ((3827,)): 2.182261505126953e-15 erg Å^-1 cm^-2 s^-1 .. 1.7559197998046877e-15 erg Å^-1 cm^-2 s^-1
meta: Dict{Symbol, Any}()
Expected output:
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (3815.0483f0 Å, 9206.613f0 Å)
flux: (9.416705322265625e-16 erg Å^-1 cm^-2 s^-1, 9.156107177734375e-15 erg Å^-1 cm^-2 s^-1)
meta: Dict{Symbol, Any}()
diff =
Warning: Diff output requires color.
SingleSpectrum(Quantity{Float32, 𝐋, Unitful.FreeUnits{(Å,), 𝐋, nothing}}, Quantity{Float64, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(Å^-1, erg, cm^-2, s^-1), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (3815.0483f0 Å, wave ((3827,)): 3815.0483f0 Å .. 9206.613f0 Å)
flux: (9.416705322265625e-16 Å
flux ((3827,)): 2.182261505126953e-15 erg Å^-1 cm^-2 s^-1, 9.156107177734375e-15 s^-1 .. 1.7559197998046877e-15 erg Å^-1 cm^-2 s^-1)
s^-1
meta: Dict{Symbol, Any}()
|
|
build:
docs/src/transforms.md#L14
doctest failure in docs/src/transforms.md:14-66
```jldoctest
julia> using Unitful, Measurements, Random
julia> rng = Random.seed!(0);
julia> wave = (1:0.5:3)u"μm"
(1.0:0.5:3.0) μm
julia> sigma = randn(rng, size(wave))
5-element Vector{Float64}:
0.942970533446119
0.13392275765318448
1.5250689085124804
0.12390123120559722
-1.205772284259936
julia> flux = (100 .± sigma)u"W/m^2/μm"
5-element Vector{Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}}}:
100.0 ± 0.94 W μm^-1 m^-2
100.0 ± 0.13 W μm^-1 m^-2
100.0 ± 1.5 W μm^-1 m^-2
100.0 ± 0.12 W μm^-1 m^-2
100.0 ± -1.2 W μm^-1 m^-2
julia> spec = spectrum(wave, flux)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± -1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red = redden(spec, 0.3)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (89.44 ± 0.84 W μm^-1 m^-2, 98.1 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red.flux
5-element Vector{Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}}}:
89.44 ± 0.84 W μm^-1 m^-2
94.35 ± 0.13 W μm^-1 m^-2
96.4 ± 1.5 W μm^-1 m^-2
97.48 ± 0.12 W μm^-1 m^-2
98.1 ± 1.2 W μm^-1 m^-2
julia> deredden!(red, 0.3)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red.flux ≈ spec.flux
true
```
Subexpression:
deredden!(red, 0.3)
Evaluated output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave ((5,)): 1.0 μm .. 3.0 μm
flux ((5,)): 100.0 ± 0.94 W μm^-1 m^-2 .. 100.0 ± 1.2 W μm^-1 m^-2
meta: Dict{Symbol, Any}()
Expected output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
diff =
Warning: Diff output requires color.
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, wave ((5,)): 1.0 μm .. 3.0 μm)
flux: (100.0 μm
flux ((5,)): 100.0 ± 0.94 W μm^-1 m^-2, m^-2 .. 100.0 ± 1.2 W μm^-1 m^-2)
m^-2
meta: Dict{Symbol, Any}()
|
|
build:
docs/src/transforms.md#L14
doctest failure in docs/src/transforms.md:14-66
```jldoctest
julia> using Unitful, Measurements, Random
julia> rng = Random.seed!(0);
julia> wave = (1:0.5:3)u"μm"
(1.0:0.5:3.0) μm
julia> sigma = randn(rng, size(wave))
5-element Vector{Float64}:
0.942970533446119
0.13392275765318448
1.5250689085124804
0.12390123120559722
-1.205772284259936
julia> flux = (100 .± sigma)u"W/m^2/μm"
5-element Vector{Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}}}:
100.0 ± 0.94 W μm^-1 m^-2
100.0 ± 0.13 W μm^-1 m^-2
100.0 ± 1.5 W μm^-1 m^-2
100.0 ± 0.12 W μm^-1 m^-2
100.0 ± -1.2 W μm^-1 m^-2
julia> spec = spectrum(wave, flux)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± -1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red = redden(spec, 0.3)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (89.44 ± 0.84 W μm^-1 m^-2, 98.1 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red.flux
5-element Vector{Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}}}:
89.44 ± 0.84 W μm^-1 m^-2
94.35 ± 0.13 W μm^-1 m^-2
96.4 ± 1.5 W μm^-1 m^-2
97.48 ± 0.12 W μm^-1 m^-2
98.1 ± 1.2 W μm^-1 m^-2
julia> deredden!(red, 0.3)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red.flux ≈ spec.flux
true
```
Subexpression:
red = redden(spec, 0.3)
Evaluated output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave ((5,)): 1.0 μm .. 3.0 μm
flux ((5,)): 89.44 ± 0.84 W μm^-1 m^-2 .. 98.1 ± 1.2 W μm^-1 m^-2
meta: Dict{Symbol, Any}()
Expected output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (89.44 ± 0.84 W μm^-1 m^-2, 98.1 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
diff =
Warning: Diff output requires color.
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, wave ((5,)): 1.0 μm .. 3.0 μm)
flux: (89.44 μm
flux ((5,)): 89.44 ± 0.84 W μm^-1 m^-2, m^-2 .. 98.1 ± 1.2 W μm^-1 m^-2)
m^-2
meta: Dict{Symbol, Any}()
|
|
build:
docs/src/transforms.md#L14
doctest failure in docs/src/transforms.md:14-66
```jldoctest
julia> using Unitful, Measurements, Random
julia> rng = Random.seed!(0);
julia> wave = (1:0.5:3)u"μm"
(1.0:0.5:3.0) μm
julia> sigma = randn(rng, size(wave))
5-element Vector{Float64}:
0.942970533446119
0.13392275765318448
1.5250689085124804
0.12390123120559722
-1.205772284259936
julia> flux = (100 .± sigma)u"W/m^2/μm"
5-element Vector{Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}}}:
100.0 ± 0.94 W μm^-1 m^-2
100.0 ± 0.13 W μm^-1 m^-2
100.0 ± 1.5 W μm^-1 m^-2
100.0 ± 0.12 W μm^-1 m^-2
100.0 ± -1.2 W μm^-1 m^-2
julia> spec = spectrum(wave, flux)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± -1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red = redden(spec, 0.3)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (89.44 ± 0.84 W μm^-1 m^-2, 98.1 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red.flux
5-element Vector{Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}}}:
89.44 ± 0.84 W μm^-1 m^-2
94.35 ± 0.13 W μm^-1 m^-2
96.4 ± 1.5 W μm^-1 m^-2
97.48 ± 0.12 W μm^-1 m^-2
98.1 ± 1.2 W μm^-1 m^-2
julia> deredden!(red, 0.3)
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± 1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
julia> red.flux ≈ spec.flux
true
```
Subexpression:
spec = spectrum(wave, flux)
Evaluated output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave ((5,)): 1.0 μm .. 3.0 μm
flux ((5,)): 100.0 ± 0.94 W μm^-1 m^-2 .. 100.0 ± -1.2 W μm^-1 m^-2
meta: Dict{Symbol, Any}()
Expected output:
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, 3.0 μm)
flux: (100.0 ± 0.94 W μm^-1 m^-2, 100.0 ± -1.2 W μm^-1 m^-2)
meta: Dict{Symbol, Any}()
diff =
Warning: Diff output requires color.
SingleSpectrum(Quantity{Float64, 𝐋, Unitful.FreeUnits{(μm,), 𝐋, nothing}}, Quantity{Measurement{Float64}, 𝐌 𝐋^-1 𝐓^-3, Unitful.FreeUnits{(μm^-1, m^-2, W), 𝐌 𝐋^-1 𝐓^-3, nothing}})
wave: (1.0 μm, wave ((5,)): 1.0 μm .. 3.0 μm)
flux: (100.0 μm
flux ((5,)): 100.0 ± 0.94 W μm^-1 m^-2, m^-2 .. 100.0 ± -1.2 W μm^-1 m^-2)
m^-2
meta: Dict{Symbol, Any}()
|