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)) |
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
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
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 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
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
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
## 2D Gross-Pitaevskii equation, with magnetic field, and with | |
## multiple electrons (of course, it doesn't make physical sense, but | |
## why not) | |
## This is pretty WIP, and only serves as a very rough demo. Nothing | |
## has been checked properly, so do not use for any serious purposes. | |
using PyCall | |
using DFTK | |
using Printf |
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
## 2D Gross-Pitaevskii equation, with magnetic field, and with | |
## multiple electrons (of course, it doesn't make physical sense, but | |
## why not) | |
## This is pretty WIP, and only serves as a very rough demo. Nothing | |
## has been checked properly, so do not use for any serious purposes. | |
using PyCall | |
using DFTK | |
using Printf |
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 FFTW | |
n = 128 | |
m = 5 | |
X = randn(ComplexF64, n, n, n, m) | |
p = plan_fft!(X, (1, 2, 3)) | |
Y = copy(X) | |
GC.gc() | |
GC.enable(false) |
NewerOlder