update sample for spm-pade

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.

samples/fermion_twopeak/param.in → 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

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