Stefanos Carlström

Complex derivative of power fails

7

I'm trying to calculate complex derivatives of analytic functions, using [a trick](https://discourse.julialang.org/t/automatic-differentiation-of-complex-valued-functions/30263/3) by @tpapp: ```julia using ForwardDiff ∂(f, x::Real) = ForwardDiff.derivative(f, x) function ∂(f, z::Complex) # https://discourse.julialang.org/t/automatic-differentiation-of-complex-valued-functions/30263/3 ff = ((x,y),)...

Banded matrix multiplication does not LoopVectorize

5

I have the slightly unusual use-case dense matrix = dense matrix * banded matrix (actually more complicated than this, but this shows the issue). With the help of @mcabbott I...

Banded concretization of derivatives with boundary conditions gives dense matrices

1

Consider the following example: ```julia julia> L = 10.0 10.0 julia> N = 9 9 julia> δx = L/(N+1) 1.0 julia> Δ = CenteredDifference(2, 2, δx, N) DerivativeOperator{Float64,1,false,Float64,StaticArrays.SArray{Tuple{3},Float64,1,3},StaticArrays.SArray{Tuple{0},StaticArrays.SArray{Tuple{4},Float64,1,4},1,0},Nothing,Nothing}(2, 2, 1.0,...

[WIP] Correct concretization BCs and ghost derivative operators

2

Fixes #182. So far, the different concretizations of boundary conditions have been implemented, ghost derivative operators are up next. Please check `u` for `PeriodicBC`s, I was unsure if they should...

Added @materialize convenience DSL

13

As sketched in https://github.com/JuliaApproximation/ContinuumArrays.jl/pull/37

Small numbers formatted incorrectly or crash

2

```julia julia> using Formatting julia> fmt = FormatExpr("{1:+11.3e}") FormatExpr("", "", Formatting.FormatEntry[Formatting.FormatEntry(Formatting.ArgSpec(1, false, identity), FormatSpec cls = f typ = e fill = align = > sign = + width =...

fullcite command to print citation in-place

1

This package seems awesome! Would it be nice to have a [`fullcite`](https://tex.stackexchange.com/a/2922/9742) command that prints the whole reference in-place?

findall for Date(Time) ranges

1

Many thanks for getting this done. I just found another edge case, which I'm not really sure this ought to try to support, but anyway here it is: ```julia julia>...

Other algorithms

4

Is this package the place where we'd like to put other algorithms as well? Such as - Chebyshev (requires prior knowledge of spectral range of the operator) - [Padé](https://en.wikipedia.org/wiki/Padé_table#An_example_–_the_exponential_function) (explicit...

Intermittent crashes

1

Every so often, a simple script like ```julia using Blender using GeometryTypes using FileIO using AffineTransforms Blender.clear_objects() b = load("bunny.obj") trf = tformscale(15.f0,3)*tformrotate([1.f0,0.f0,0.f0], Float32(π)/2) ob = Blender.add(b, "bunny", trf) println(ob)...