jwscook

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())

jwscook

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

jwscook

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);

jwscook

Two stream instability

Electrostatic linear plasma dispersion relation

$$ 1 + \Sigma_s \frac{\Pi_s^2}{\omega}\int_{-\infty}^{\infty} \frac{\partial f_s(v_\parallel)}{\partial v_\parallel}\frac{v_\parallel}{\omega-k_\parallel v_\parallel} dv_\parallel=0 $$

where, in SI units:

jwscook

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

jwscook

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

jwscook

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

jwscook

using ForwardDiff: Dual, Tag, jacobian, value, tagtype, partials, gradient
using LinearAlgebra
using NLsolve
using Random
using Test

Random.seed!(0)

const jac = jacobian

jwscook

module BiComplexes

using AutoHashEquals

export BiComplex, derivative

@auto_hash_equals struct BiComplex{T} <: Number
  re::Complex{T}
  im::Complex{T}
end

jwscook

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