From e7faf7f1c7f04149627a2d3dd77edb23dc17a561 Mon Sep 17 00:00:00 2001
From: Aaron LI <aaronly.me@outlook.com>
Date: Tue, 5 Apr 2016 15:59:56 +0800
Subject: sbp_fit.py: allow fix parameter

---
 python/sbp_fit.py | 17 +++++++++++++----
 1 file changed, 13 insertions(+), 4 deletions(-)

(limited to 'python/sbp_fit.py')

diff --git a/python/sbp_fit.py b/python/sbp_fit.py
index 0ce8d4d..c9a67a5 100644
--- a/python/sbp_fit.py
+++ b/python/sbp_fit.py
@@ -3,9 +3,11 @@
 #
 # Aaron LI
 # Created: 2016-03-13
-# Updated: 2016-03-31
+# Updated: 2016-04-05
 #
 # Changelogs:
+# 2016-04-05:
+#   * Allow fix parameters
 # 2016-03-31:
 #   * Remove `ci_report()'
 #   * Add `make_results()' to orgnize all results as s Python dictionary
@@ -21,6 +23,9 @@
 #   * Add `ignore_data()' and `notice_data()' support
 #   * Add classes `FitModelSBetaNorm' and `FitModelDBetaNorm'
 #
+# TODO:
+#   * to allow fit the outer beta component, then fix it, and fit the inner one
+#
 
 """
 Fit the surface brightness profile (SBP) with the single-beta model:
@@ -66,8 +71,9 @@ bkg   = 1.0e-9, 0.0,  1.0e-7
 -------------------------------------------------
 """
 
-__version__ = "0.5.0"
-__date__    = "2016-03-31"
+__version__ = "0.5.1"
+__date__    = "2016-04-05"
+
 
 import numpy as np
 import lmfit
@@ -581,8 +587,11 @@ def make_model(config):
     # set initial values and bounds for the model parameters
     params = config.get(modelname)
     for p, value in params.items():
+        variable = True
+        if len(value) == 4 and value[3].upper() in ["FIXED", "FALSE"]:
+            variable = False
         model.set_param(name=p, value=float(value[0]),
-                min=float(value[1]), max=float(value[2]))
+                min=float(value[1]), max=float(value[2]), vary=variable)
     return model
 
 
-- 
cgit v1.2.2