Update to v0.4.2:v0.4.2

  get simpleq-Kv to work again

20 files changed, 36 insertions, 42 deletions

diff --git a/Changelog b/Changelog
index ce743bc..e8ae755 100644
--- a/Changelog
+++ b/Changelog
@@ -1,3 +1,6 @@
+0.4.2:
+  get simpleq-Kv to work again
+
 0.4.1:
   Updated the references in the documentation
 
diff --git a/doc/Makefile b/doc/Makefile
index d0c4c64..8221979 100644
--- a/doc/Makefile
+++ b/doc/Makefile
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/doc/libs/code.sty b/doc/libs/code.sty
index e067dfe..8bf3c41 100644
--- a/doc/libs/code.sty
+++ b/doc/libs/code.sty
@@ -1,4 +1,4 @@
-%% Copyright 2021-2023 Ian Jauslin
+%% Copyright 2021-2024 Ian Jauslin
 %% 
 %% Licensed under the Apache License, Version 2.0 (the "License");
 %% you may not use this file except in compliance with the License.
diff --git a/doc/libs/dlmf.sty b/doc/libs/dlmf.sty
index f66157b..6422fa4 100644
--- a/doc/libs/dlmf.sty
+++ b/doc/libs/dlmf.sty
@@ -1,4 +1,4 @@
-%% Copyright 2021-2023 Ian Jauslin
+%% Copyright 2021-2024 Ian Jauslin
 %% 
 %% Licensed under the Apache License, Version 2.0 (the "License");
 %% you may not use this file except in compliance with the License.
diff --git a/doc/libs/iantheo.sty b/doc/libs/iantheo.sty
index 31842a4..0c20ea8 100644
--- a/doc/libs/iantheo.sty
+++ b/doc/libs/iantheo.sty
@@ -1,4 +1,4 @@
-%% Copyright 2021-2023 Ian Jauslin
+%% Copyright 2021-2024 Ian Jauslin
 %% 
 %% Licensed under the Apache License, Version 2.0 (the "License");
 %% you may not use this file except in compliance with the License.
diff --git a/doc/libs/largearray.sty b/doc/libs/largearray.sty
index 3a8762e..a8a63c3 100644
--- a/doc/libs/largearray.sty
+++ b/doc/libs/largearray.sty
@@ -1,4 +1,4 @@
-%% Copyright 2021-2023 Ian Jauslin
+%% Copyright 2021-2024 Ian Jauslin
 %% 
 %% Licensed under the Apache License, Version 2.0 (the "License");
 %% you may not use this file except in compliance with the License.
diff --git a/doc/libs/point.sty b/doc/libs/point.sty
index 796e0d2..e4fb119 100644
--- a/doc/libs/point.sty
+++ b/doc/libs/point.sty
@@ -1,4 +1,4 @@
-%% Copyright 2021-2023 Ian Jauslin
+%% Copyright 2021-2024 Ian Jauslin
 %% 
 %% Licensed under the Apache License, Version 2.0 (the "License");
 %% you may not use this file except in compliance with the License.
diff --git a/doc/simplesolv-doc.tex b/doc/simplesolv-doc.tex
index 7e6b397..0b6e06d 100644
--- a/doc/simplesolv-doc.tex
+++ b/doc/simplesolv-doc.tex
@@ -1,4 +1,4 @@
-%% Copyright 2021-2023 Ian Jauslin
+%% Copyright 2021-2024 Ian Jauslin
 %% 
 %% Licensed under the Apache License, Version 2.0 (the "License");
 %% you may not use this file except in compliance with the License.
@@ -27,7 +27,7 @@
 \hfil{\bf\Large
 {\tt simplesolv}\par
 \medskip
-\hfil \large v0.4.1
+\hfil \large v0.4.2
 }
 
 \vfill
diff --git a/src/anyeq.jl b/src/anyeq.jl
index fee77d8..79d1ddd 100644
--- a/src/anyeq.jl
+++ b/src/anyeq.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/chebyshev.jl b/src/chebyshev.jl
index af4be40..f3103d0 100644
--- a/src/chebyshev.jl
+++ b/src/chebyshev.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/easyeq.jl b/src/easyeq.jl
index 2dbbd1a..ea44d9e 100644
--- a/src/easyeq.jl
+++ b/src/easyeq.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/integration.jl b/src/integration.jl
index 223d6cc..4d10138 100644
--- a/src/integration.jl
+++ b/src/integration.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/interpolation.jl b/src/interpolation.jl
index 066bc20..d3a8161 100644
--- a/src/interpolation.jl
+++ b/src/interpolation.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/main.jl b/src/main.jl
index 28fc2be..9de9a0a 100644
--- a/src/main.jl
+++ b/src/main.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
@@ -439,7 +439,7 @@ function main()
     if command=="2pt"
       simpleq_Kv_2pt(minlrho_init,nlrho_init,taus,P,N,J,rho,a0,v,maxiter,tolerance,xmin,xmax,nx)
     elseif command=="condensate_fraction_rho"
-      simpleq_Kv_condensate_fraction(rhos,taus,P,N,J,a0,v,maxiter,tolerance)
+      simpleq_Kv_condensate_fraction(minlrho_init,nlrho_init,rhos,taus,P,N,J,a0,v,maxiter,tolerance)
     elseif command=="Kv"
       simpleq_Kv_Kv(minlrho_init,nlrho_init,taus,P,N,J,rho,a0,v,maxiter,tolerance,xmin,xmax,nx)
     else
diff --git a/src/potentials.jl b/src/potentials.jl
index b480db0..01ce845 100644
--- a/src/potentials.jl
+++ b/src/potentials.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/print.jl b/src/print.jl
index 3587728..b7f2215 100644
--- a/src/print.jl
+++ b/src/print.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/simpleq-Kv.jl b/src/simpleq-Kv.jl
index 5a6579c..e61e464 100644
--- a/src/simpleq-Kv.jl
+++ b/src/simpleq-Kv.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
@@ -14,8 +14,8 @@
 
 # Compute Kv=(-\Delta+v+4e(1-\rho u*))^{-1}v
 function anyeq_Kv(
-  minlrho::Float64,
-  nlrho::Int64,
+  minlrho_init::Float64,
+  nlrho_init::Int64,
   taus::Array{Float64,1},
   P::Int64,
   N::Int64,
@@ -30,15 +30,10 @@ function anyeq_Kv(
   nx::Int64
 )
   # init vectors
-  (weights,T,k,V,V0,A,Upsilon,Upsilon0)=anyeq_init(taus,P,N,J,v)
+  (weights,T,k,V,V0,A,Abar,Upsilon,Upsilon0)=anyeq_init(taus,P,N,J,v,Anyeq_approx(0.,0.,1.,0.,0.,0.,0.,0.,0.,0.,0.),"")
 
   # compute initial guess from medeq
-  rhos=Array{Float64,1}(undef,nlrho)
-  for j in 0:nlrho-1
-    rhos[j+1]=(nlrho==1 ? rho : 10^(minlrho+(log10(rho)-minlrho)/(nlrho-1)*j))
-  end
-  u0s=anyeq_init_medeq(rhos,N,J,k,a0,v,maxiter,tolerance)
-  u0=u0s[nlrho]
+  u0=anyeq_init_medeq([rho],minlrho_init,nlrho_init,N,J,k,a0,v,maxiter,tolerance)
 
   (u,E,error)=anyeq_hatu(u0,P,N,J,rho,a0,weights,k,taus,V,V0,A,Abar,Upsilon,Upsilon0,v,maxiter,tolerance,Anyeq_approx(0.,0.,1.,0.,0.,0.,0.,0.,0.,0.,0.))
 
@@ -60,6 +55,8 @@ end
 
 # Compute the condensate fraction for simpleq using Kv
 function simpleq_Kv_condensate_fraction(
+  minlrho_init::Float64,
+  nlrho_init::Float64,
   rhos::Array{Float64,1},
   taus::Array{Float64,1},
   P::Int64,
@@ -71,10 +68,10 @@ function simpleq_Kv_condensate_fraction(
   tolerance::Float64
 )
   # init vectors
-  (weights,T,k,V,V0,A,Upsilon,Upsilon0)=anyeq_init(taus,P,N,J,v)
+  (weights,T,k,V,V0,A,Abar,Upsilon,Upsilon0)=anyeq_init(taus,P,N,J,v,Anyeq_approx(0.,0.,1.,0.,0.,0.,0.,0.,0.,0.,0.),"")
 
   # compute initial guess from medeq
-  u0s=anyeq_init_medeq(rhos,N,J,k,a0,v,0.,maxiter,tolerance)
+  u0s=anyeq_init_medeq(rhos,minlrho_init,nlrho_init,N,J,k,a0,v,0.,maxiter,tolerance)
 
   for j in 1:length(rhos)
     (u,E,error)=anyeq_hatu(u0s[j],0.,P,N,J,rhos[j],a0,weights,k,taus,V,V0,A,Abar,Upsilon,Upsilon0,v,maxiter,tolerance,Anyeq_approx(0.,0.,1.,0.,0.,0.,0.,0.,0.,0.,0.))
@@ -93,8 +90,8 @@ end
 
 # Compute the two-point correlation function for simpleq using Kv
 function simpleq_Kv_2pt(
-  minlrho::Float64,
-  nlrho::Int64,
+  minlrho_init::Float64,
+  nlrho_init::Int64,
   taus::Array{Float64,1},
   P::Int64,
   N::Int64,
@@ -109,17 +106,11 @@ function simpleq_Kv_2pt(
   nx::Int64
 )
   # init vectors
-  (weights,T,k,V,V0,A,Upsilon,Upsilon0)=anyeq_init(taus,P,N,J,v)
+  (weights,T,k,V,V0,A,Abar,Upsilon,Upsilon0)=anyeq_init(taus,P,N,J,v,Anyeq_approx(0.,0.,1.,0.,0.,0.,0.,0.,0.,0.,0.),"")
 
   # compute initial guess from medeq
-  rhos=Array{Float64,1}(undef,nlrho)
-  for j in 0:nlrho-1
-    rhos[j+1]=(nlrho==1 ? rho : 10^(minlrho+(log10(rho)-minlrho)/(nlrho-1)*j))
-  end
-  u0s=anyeq_init_medeq(rhos,N,J,k,a0,v,maxiter,tolerance)
-  u0=u0s[nlrho]
+  u0=anyeq_init_medeq([rho],minlrho_init,nlrho_init,N,J,k,a0,v,maxiter,tolerance)
 
-  Abar=Array{Float64,5}(undef,0,0,0,0,0)
   (u,E,error)=anyeq_hatu(u0,P,N,J,rho,a0,weights,k,taus,V,V0,A,Abar,Upsilon,Upsilon0,v,maxiter,tolerance,Anyeq_approx(0.,0.,1.,0.,0.,0.,0.,0.,0.,0.,0.))
 
   # Kv in Fourier space
diff --git a/src/simpleq-hardcore.jl b/src/simpleq-hardcore.jl
index 398ac06..58fccb4 100644
--- a/src/simpleq-hardcore.jl
+++ b/src/simpleq-hardcore.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/simpleq-iteration.jl b/src/simpleq-iteration.jl
index 4c4de07..5cef5c8 100644
--- a/src/simpleq-iteration.jl
+++ b/src/simpleq-iteration.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.
diff --git a/src/tools.jl b/src/tools.jl
index 5d1678d..f2336f7 100644
--- a/src/tools.jl
+++ b/src/tools.jl
@@ -1,4 +1,4 @@
-## Copyright 2021-2023 Ian Jauslin
+## Copyright 2021-2024 Ian Jauslin
 ## 
 ## Licensed under the Apache License, Version 2.0 (the "License");
 ## you may not use this file except in compliance with the License.