Electrostatic linear plasma dispersion relation
where, in SI units:
using Random, Distributed, Test | |
Random.seed!(0) | |
addprocs(4, exeflags="-t 2") | |
@everywhere begin | |
using LinearAlgebra, Random | |
using Distributed, DistributedArrays, SharedArrays | |
LinearAlgebra.BLAS.set_num_threads(Base.Threads.nthreads()) |
using LinearAlgebra, Random, Base.Threads | |
Random.seed!(0) | |
alphafactor(x::Real) = -sign(x) | |
alphafactor(x::Complex) = -exp(im * angle(x)) | |
function householder!(A::Matrix{T}) where T | |
m, n = size(A) | |
H = A | |
α = zeros(T, min(m, n)) # Diagonal of R |
using LinearAlgebra, Test, Random | |
Random.seed!(0) | |
@testset "understand raw lapack calls" begin | |
@show Threads.nthreads() | |
@show BLAS.get_num_threads() | |
m = 5 | |
n = 4 | |
A = rand(m, n); |
using SpecialFunctions, Bessels, LinearAlgebra, Plots | |
const spj = SpecialFunctions.besselj | |
const bej = Bessels.besselj | |
function foo!(A, B, f::T, ns, ms, N) where T | |
for (j, n) in enumerate(ns) | |
for (i, m) in enumerate(ms) | |
A[i, j] = 0 | |
t = @elapsed for _ in 1:N |
This file contains any messages produced by compilers while | |
running configure, to aid debugging if configure makes a mistake. | |
It was created by GNU Octave configure 6.4.0, which was | |
generated by GNU Autoconf 2.69. Invocation command line was | |
$ ./configure CFLAGS=-O3 -m64 -mavx2 -mfma -march=native -mtune=native -mfpmath=sse -malign-data=cacheline -ftree-vectorize -fno-semantic-interposition -fomit-frame-pointer -I/home/cookj/BUILT/include -L/home/cookj/BUILT/lib -Wl,-rpath=/home/cookj/BUILT/lib -pipe -fno-builtin-malloc -fno-builtin-calloc -fno-builtin-realloc -fno-builtin-free -lmimalloc CXXFLAGS=-O3 -m64 -mavx2 -mfma -march=native -mtune=native -mfpmath=sse -malign-data=cacheline -ftree-vectorize -fno-semantic-interposition -fomit-frame-pointer -I/home/cookj/BUILT/include -L/home/cookj/BUILT/lib -Wl,-rpath=/home/cookj/BUILT/lib -pipe -fno-builtin-malloc -fno-builtin-calloc -fno-builtin-realloc -fno-builtin-free -lmimalloc FFLAGS=-O3 -m64 -mavx2 -mfma -march=native -mtune=native -mfpmath=sse -malign-dat |
using OrdinaryDiffEq | |
using LinearAlgebra | |
using StaticArrays | |
using Trixi | |
import Trixi: AbstractEquations, varnames, flux | |
abstract type AbstractMaxwellEquations{NDIMS, NVARS} <: AbstractEquations{NDIMS, NVARS} end | |
struct MaxwellEquations2D{RealT<:Real} <: AbstractMaxwellEquations{2, 6} | |
μ₀::Float64 | |
ϵ₀::Float64 |
using ForwardDiff: Dual, Tag, jacobian, value, tagtype, partials, gradient | |
using LinearAlgebra | |
using NLsolve | |
using Random | |
using Test | |
Random.seed!(0) | |
const jac = jacobian |
module BiComplexes | |
using AutoHashEquals | |
export BiComplex, derivative | |
@auto_hash_equals struct BiComplex{T} <: Number | |
re::Complex{T} | |
im::Complex{T} | |
end |
using Random | |
const seed = hash(time()) | |
Random.seed!(seed) | |
searchstring = "const seed = " | |
open("$(@__FILE__)" * ".reproducible", "w") do io | |
for line in readlines(open(@__FILE__)) | |
if length(line) > length(searchstring) && line[1:length(searchstring)] == searchstring |