This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
double free or corruption (out) | |
signal (6): Aborted | |
in expression starting at /home/antoine/.julia/dev/DFTK/examples/silicon.jl:50 | |
in expression starting at /home/antoine/.julia/dev/DFTK/examples/silicon.jl:50 | |
in expres11): Segmentation fault | |
in expression starting at /home/antoine/.julia/dev/DFTK/examples/silicon.jl:50 | |
fftw_twiddle_awake at /home/antoine/.julia/artifacts/e40697527cebb56d421346210295905df6e421dc/lib/libfftw3.so (unknown line) | |
fftw_plan_awake at /home/antoine/.julia/artifacts/e40697527cebb56d421346210295905df6e421dc/lib/libfftw3.so (unknown line) | |
fftw_plan_awake at /home/antoine/.julia/artifacts/e40697527cebb56d421346210295905df6e421dc/lib/libfftw3.so (unknown line) |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
double free or corruption (fasttop) | |
signal (6): Aborted | |
in expression starting at /home/antoine/.julia/dev/DFTK/examples/silicon.jl:50 | |
__libc_signal_restore_set at /build/glibc-OTsEL5/glibc-2.27/signal/../sysdeps/unix/sysv/linux/nptl-signals.h:80 [inlined] | |
raise at /build/glibc-OTsEL5/glibc-2.27/signal/../sysdeps/unix/sysv/linux/raise.c:48 | |
abort at /build/glibc-OTsEL5/glibc-2.27/stdlib/abort.c:79 | |
__libc_message at /build/glibc-OTsEL5/glibc-2.27/libio/../sysdeps/posix/libc_fatal.c:181 | |
malloc_printerr at /build/glibc-OTsEL5/glibc-2.27/malloc/malloc.c:5350 | |
_int_free at /build/glibc-OTsEL5/glibc-2.27/malloc/malloc.c:4230 [inlined] |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using StaticArrays | |
using PyPlot | |
using QuadGK | |
const latt = @SMatrix[1 0; | |
0 1] | |
const α = .1 | |
f(r) = exp(-α*r^2) | |
fp(r) = -α * 2 * r * exp(-α*r^2) | |
Rmax = 100 |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using PseudoArcLengthContinuation, LinearAlgebra, Plots, PyPlot, Optim | |
const PALC = PseudoArcLengthContinuation | |
# fonction à annuler | |
function f(x, alpha) | |
resu = zeros(1) | |
resu[1]= alpha*x[1] + x[1]^3 | |
return resu |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using PseudoArcLengthContinuation, LinearAlgebra, Plots, PyPlot, Optim | |
const PALC = PseudoArcLengthContinuation | |
using ForwardDiff | |
function f(x, alpha) | |
[alpha*x[1] + x[1]^3 ] | |
end | |
function Jf(x, alpha) | |
Jf = zeros(1, 1) | |
Jf[1] = alpha+3x[1]^2 |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
ERROR: LoadError: MethodError: no method matching getTangent!(::BorderedArray{Array{Float64,1},Float64}, ::BorderedArray{Array{Float64,1},Float64}, ::BorderedArray{Array{Float64,1},Int64}, ::PALCIterable{typeof(F),PseudoArcLengthContinuation.var"#175#177"{typeof(F)},Array{Float64,1},NamedTuple{(:C,),Tuple{Int64}},Setfield.PropertyLens{:C},Float64,DefaultLS,DefaultEig{typeof(real)},SecantPred,BorderingBLS{DefaultLS,Float64},PseudoArcLengthContinuation.var"#155#169"{PseudoArcLengthContinuation.var"#155#156#170"},PseudoArcLengthContinuation.var"#153#167",typeof(norm),PseudoArcLengthContinuation.DotTheta{PseudoArcLengthContinuation.var"#154#168"},PseudoArcLengthContinuation.var"#157#171",PseudoArcLengthContinuation.var"#158#172"{PseudoArcLengthContinuation.var"#158#159#173"},String}, ::Float64, ::Float64, ::SecantPred, ::Int64) | |
Closest candidates are: | |
getTangent!(::M, ::M, ::M, ::PALCIterable, ::Any, ::Any, ::Talgo, ::Any) where {T, vectype, M<:BorderedArray{vectype,T}, Talgo<:PseudoArcLengthContinuation.Abstra |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
C_beg = 0. | |
C_end = -6. | |
opt_newton = PALC.NewtonPar(tol = 1e-8, verbose = true, maxIter = 20) | |
opts_br = ContinuationPar(dsmin = 0.001, dsmax = 0.05, ds = -0.01, pMax = C_beg, pMin = C_end, | |
detectBifurcation = 2, nev = 30, plotEveryNsteps = 10, newtonOptions = opt_newton, | |
maxSteps = 100, precisionStability = 1e-6, nInversion = 4, dsminBisection = 1e-7, maxBisectionSteps = 25) | |
ρ0 = zeros(prod(basis.fft_size)) |
OlderNewer