From a550be45e065c5386a056173c6040ee05b1f7233 Mon Sep 17 00:00:00 2001
From: Yuichi Motoyama <yomichi@tsg.jp>
Date: Tue, 21 Sep 2021 13:17:45 +0900
Subject: [PATCH] update sample for spm-pade

---
 docs/manual/source/docs/tutorials.rst         |  4 ++++
 .../{param.in => param_spm.in}                |  2 +-
 samples/fermion_twopeak/param_spmpade.in      | 20 +++++++++++++++++++
 3 files changed, 25 insertions(+), 1 deletion(-)
 rename samples/fermion_twopeak/{param.in => param_spm.in} (95%)
 create mode 100644 samples/fermion_twopeak/param_spmpade.in

diff --git a/docs/manual/source/docs/tutorials.rst b/docs/manual/source/docs/tutorials.rst
index 90bc55a..ad50101 100644
--- a/docs/manual/source/docs/tutorials.rst
+++ b/docs/manual/source/docs/tutorials.rst
@@ -206,6 +206,10 @@ The strength of the Gaussian noise is specified by ``g_sigma`` or the pair of ``
 The coefficient of the Pade weight, :math:`eta`, is specified by the ``PadeEta`` parameter.
 If ``PadeEta = 0``, the original SpM method will be used.
 
+``param_spm.in`` and ``param_spmpade.in`` are input files for the SpM (:math:`\eta=0`) and SpM-Pade (:math:`\eta=1`) method, respectively.
+Note that optimal :math:`\lambda` for the SpM and SpM-Pade differ, and hence users may need to change the range of :math:`lambda` (``lambdalogbegin`` and ``lambdalogend``).
+The following figure shows the result:
+
 .. image:: figs/spectrum_spmpade.png
 
 The blue curve and red curve show the spectrum reconstructed by the SpM method (:math:`\eta=0`) and the SpM-Pade method (:math:`\eta=1`), respectively.
diff --git a/samples/fermion_twopeak/param.in b/samples/fermion_twopeak/param_spm.in
similarity index 95%
rename from samples/fermion_twopeak/param.in
rename to samples/fermion_twopeak/param_spm.in
index 736168d..fa92d5a 100644
--- a/samples/fermion_twopeak/param.in
+++ b/samples/fermion_twopeak/param_spm.in
@@ -16,5 +16,5 @@ lambdalogend=-10.0
 tolerance=1e-10
 maxiteration=10000
 # Pade
-PadeEta = 1.0
+PadeEta = 0.0
 g_sigma=1e-6
diff --git a/samples/fermion_twopeak/param_spmpade.in b/samples/fermion_twopeak/param_spmpade.in
new file mode 100644
index 0000000..cb0b710
--- /dev/null
+++ b/samples/fermion_twopeak/param_spmpade.in
@@ -0,0 +1,20 @@
+# INPUT/OUTPUT
+statistics="fermion"
+beta=100
+filein_G="Gtau.in"
+column=1  
+fileout_spec="spectrum.dat"
+fileout_pade="pade.dat"
+# OMEGA
+Nomega=1001
+omegamin=-4
+omegamax=4
+SVmin=1e-12
+# ADMM
+lambdalogbegin=-6.0
+lambdalogend=-12.0
+tolerance=1e-10
+maxiteration=10000
+# Pade
+PadeEta = 1.0
+g_sigma=1e-6