RTModel 2.0__py3-none-any.whl

RTModel/RTModel.py

@@ -0,0 +1,293 @@
+
+import subprocess
+import os
+import sys
+import glob
+import time
+import inspect
+from tqdm import tqdm
+import shutil
+
+
+class RTModel:
+    def __init__(self, event = None):
+        # Directory preliminaries
+        print('*********************')
+        print('****   RTModel   ****')
+        print('*********************')
+        self.pathtoRTM = inspect.getfile(RTModel)
+        self.bindir = os.path.dirname(self.pathtoRTM) + '/bin/'
+        if os.name =='nt':
+            self.readerexe = 'Reader.exe'
+            self.initcondexe = 'InitCond.exe'
+            self.levmarexe = 'LevMar.exe'
+            self.modelselectorexe = 'ModelSelector.exe'
+            self.finalizerexe = 'Finalizer.exe'
+        else:
+            self.readerexe = 'Reader'
+            self.initcondexe = 'InitCond'
+            self.levmarexe = 'LevMar'
+            self.modelselectorexe = 'ModelSelector'
+            self.finalizerexe = 'Finalizer'
+        if(event == None):
+            self.eventname = os.getcwd()
+        else:
+            self.eventname = os.path.realpath(event)
+            print("Event name: " + self.eventname)
+        self.inidir = "ini"
+        self.modelcodes = ['PS', 'PX', 'BS', 'BO', 'LS', 'LX', 'LO']
+        self.endphase =  len(self.modelcodes)*2+3
+        self.eventinifile = 'event.ini'        
+        self.nprocessors = os.cpu_count()
+        print('Number of processors: {}'.format(self.nprocessors))
+        self.config_Reader()
+        self.config_InitCond()
+        self.config_LevMar()
+        self.config_ModelSelector()
+        self.satellitedir = '.'
+
+    def set_processors(self, nprocessors):
+        self.nprocessors = nprocessors
+
+    def set_event(self, event):
+        self.eventname = os.path.realpath(event)
+
+    def set_satellite_dir(self, satellitedir):
+        self.satellitedir = satellitedir
+
+    def config_Reader(self, tau = 0.1, binning = 4000, otherseasons = 1, renormalize = 1, thresholdoutliers = 10):
+        self.Reader_tau= tau # conventional correlation time for consecutive points
+        self.Reader_binning = binning # maximum number of points left after re-binning
+        self.Reader_otherseasons = otherseasons # How to use other seasons (0 = Yes, 1 = decrease significance, 2 = remove)
+        self.Reader_renormalize = renormalize # Re-normalize error bars if non-zero
+        self.Reader_thresholdoutliers = thresholdoutliers # Threshold in sigmas for removing outliers
+        
+    def Reader(self):
+        if(not os.path.exists(self.eventname + '/' + self.inidir)):
+            os.makedirs(self.eventname + '/' + self.inidir)
+        with open(self.eventname + '/' + self.inidir + '/Reader.ini','w') as f:
+            f.write('tau = ' + str(self.Reader_tau) + '\n')
+            f.write('binning = ' + str(self.Reader_binning) + '\n')
+            f.write('otherseasons = ' + str(self.Reader_otherseasons) + '\n')
+            f.write('renormalize = ' + str(self.Reader_renormalize) + '\n')
+            f.write('thresholdoutliers = ' + str(self.Reader_thresholdoutliers) + '\n')
+        print('- Launching: Reader')
+        print('  Pre-processing data...')
+        completedprocess = subprocess.run([self.bindir+self.readerexe,self.eventname], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL)
+        if(completedprocess.returncode != 0):
+            print('! Error in pre-processing. Please check your data!')            
+            self.done = True
+        else:
+            print('  OK')
+
+    def config_InitCond(self, npeaks = 2, peakthreshold = 10.0, oldmodels = 4, override = None, nostatic = False, onlyorbital = False, usesatellite = 0):
+        self.InitCond_npeaks = npeaks # Number of peaks in the observed light curve to be considered for setting initial conditions.
+        self.InitCond_peakthreshold = peakthreshold # Number of sigmas necessary for a deviation to be identified as a maximum or a minimum.
+        self.InitCond_oldmodels = oldmodels # Maximum number of old models to include in new run as initial conditions
+        self.InitCond_override = override # Override peak identification and manually set peak times
+        self.InitCond_nostatic = nostatic or onlyorbital # No static models will be calculated.
+        self.InitCond_noparallax = onlyorbital; # Only orbital motion models will be calculated.
+        self.InitCond_usesatellite = usesatellite; # Satellite to be used for initial conditions. Ground telescopes by default.
+        
+    def InitCond(self):
+        if(not os.path.exists(self.eventname + '/' + self.inidir)):
+            os.makedirs(self.eventname + '/' + self.inidir)
+        with open(self.eventname + '/' + self.inidir + '/InitCond.ini','w') as f:
+            f.write('npeaks = ' + str(self.InitCond_npeaks) + '\n')
+            f.write('peakthreshold = ' + str(self.InitCond_peakthreshold) + '\n')
+            f.write('oldmodels = ' + str(self.InitCond_oldmodels) + '\n')
+            f.write('usesatellite = ' + str(self.InitCond_usesatellite) + '\n')
+            if(self.InitCond_nostatic):            
+                f.write('nostatic = 1\n')
+            if(self.InitCond_noparallax):            
+                f.write('noparallax = 1\n')
+            if(self.InitCond_override != None):
+                f.write('override = ' + str(self.InitCond_override[0])+ ' ' + str(self.InitCond_override[1]) + '\n')            
+        print('- Launching: InitCond')
+        print('  Setting initial conditions...')
+        completedprocess = subprocess.run([self.bindir+self.initcondexe,self.eventname], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL)
+        if(completedprocess.returncode != 0):
+            print('! Error in setting initial conditions!')
+            self.done = True
+        else:
+            print('  OK')
+
+    def config_LevMar(self, nfits = 5, timelimit = 600.0, maxsteps = 50, bumperpower = 2.0):
+        self.LevMar_nfits = nfits # Number of models to be calculated from the same initial condition using the bumper method
+        self.LevMar_maxsteps = maxsteps # Maximum number of steps in each fit
+        self.LevMar_timelimit = timelimit # Maximum time in seconds for total execution
+        self.LevMar_bumperpower = bumperpower # Repulsion factor of bumpers
+    
+    def LevMar(self,strmodel):
+        if(not os.path.exists(self.eventname + '/' + self.inidir)):
+            os.makedirs(self.eventname + '/' + self.inidir)
+        with open(self.eventname + '/' + self.inidir + '/LevMar.ini','w') as f:
+            f.write('nfits = ' + str(self.LevMar_nfits) + '\n')
+            f.write('maxsteps = ' + str(self.LevMar_maxsteps) + '\n')
+            f.write('timelimit = ' + str(self.LevMar_timelimit) + '\n')
+            f.write('bumperpower = ' + str(self.LevMar_bumperpower) + '\n')
+        print('- Launching: LevMar')
+        print('  Fitting ' + strmodel + ' ...')
+        completedprocess = subprocess.run([self.bindir+self.levmarexe,self.eventname, strmodel,self.satellitedir], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL)
+        if(completedprocess.returncode != 0):
+            print('! Error in fit!')
+            self.done = True
+        else:
+            print('  OK')
+  
+    def launch_fits(self,modelcode):
+        if(not os.path.exists(self.eventname + '/' + self.inidir)):
+            os.makedirs(self.eventname + '/' + self.inidir)
+        with open(self.eventname + '/' + self.inidir + '/LevMar.ini','w') as f:
+            f.write('nfits = ' + str(self.LevMar_nfits) + '\n')
+            f.write('maxsteps = ' + str(self.LevMar_maxsteps) + '\n')
+            f.write('timelimit = ' + str(self.LevMar_timelimit) + '\n')
+            f.write('bumperpower = ' + str(self.LevMar_bumperpower) + '\n')
+        stringfits = {'PS' : '- Single-lens-Single-source fits',
+                      'PX' : '- Single-lens-Single-source fits with parallax',
+                      'BS' : '- Single-lens-Binary-source fits',
+                      'BO' : '- Single-lens-Binary-source fits with xallarap',
+                      'LS' : '- Binary-lens-Single-source fits',
+                      'LX' : '- Binary-lens-Single-source fits with parallax',
+                      'LO' : '- Binary-lens-Single-source fits with orbital motion'}       
+        print(stringfits[modelcode])
+        initcondfile = self.eventname + '/InitCond/' + 'InitCond'+ modelcode + '.txt'
+        if(os.path.exists(initcondfile)):
+            with open(self.eventname + '/InitCond/' + 'InitCond'+ modelcode + '.txt') as f:
+                line = f.readline().split()
+                npeaks = int(line[0])
+                ninitconds = int(line[1])   
+            processes = []
+            procnumbers = []
+            iinitcond = 0
+            finitcond = 0
+            finitcondold = -1       
+            pbar = tqdm(total = ninitconds,desc = 'Fits completed',file=sys.stdout, colour='GREEN', smoothing = 0)
+            while(finitcond < ninitconds):
+                i=0
+                while i < len(processes):
+                    if(processes[i].poll() != None):
+                        processes.pop(i)
+                        procnumbers.pop(i)
+                        finitcond += 1
+                    else:
+                        i += 1
+                while(iinitcond < ninitconds and len(processes) < self.nprocessors):
+                    strmodel =  modelcode + '{:0>4}'.format(str(iinitcond))
+                    if(glob.glob(self.eventname +'/PreModels/' + strmodel + '/t' + strmodel + '.dat')==[]):
+                        processes.append(subprocess.Popen([self.bindir+self.levmarexe,self.eventname, strmodel,self.satellitedir], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL))
+                        procnumbers.append(iinitcond)
+                    else:
+                        finitcond += 1
+                    iinitcond += 1
+                if(finitcond != finitcondold):
+                    #print('  Fits launched: {}; completed: {}/{}'.format(iinitcond, finitcond, ninitconds))
+                    pbar.update(finitcond - max(finitcondold,0))
+                    finitcondold =finitcond
+                time.sleep(0.1)
+            pbar.close()
+        else:
+            print('- No initial conditions for this class')
+ 
+    def config_ModelSelector(self, sigmasoverlap = 3.0, sigmachisquare = 1.0, maxmodels = 10):
+        self.ModelSelector_sigmasoverlap = sigmasoverlap # factor multiplying the inverse covariance in search for superpositions (models are incompatible if farther than sigmasoverlap*sigma)
+        self.ModelSelector_sigmachisquare = sigmachisquare # number of sigmas in the chi square distribution for accepting alternative models after the best one
+        self.ModelSelector_maxmodels = maxmodels # maximum number of models returned
+    
+    def ModelSelector(self, modelcode):
+        if(not os.path.exists(self.eventname + '/' + self.inidir)):
+            os.makedirs(self.eventname + '/' + self.inidir)
+        with open(self.eventname + '/' + self.inidir + '/ModelSelector.ini','w') as f:
+            f.write('sigmasoverlap = ' + str(self.ModelSelector_sigmasoverlap) + '\n')
+            f.write('sigmachisquare = ' + str(self.ModelSelector_sigmachisquare) + '\n')
+            f.write('maxmodels = ' + str(self.ModelSelector_maxmodels) + '\n')
+        stringmodels = {'PS' : '- Selecting models for Single-lens-Single-source fits',
+                        'PX' : '- Selecting models for Single-lens-Single-source fits with parallax',
+                        'BS' : '- Selecting models for Single-lens-Binary-source fits',
+                        'BO' : '- Selecting models for Single-lens-Binary-source fits with xallarap',
+                        'LS' : '- Selecting models for Binary-lens-Single-source fits',
+                        'LX' : '- Selecting models for Binary-lens-Single-source fits with parallax',
+                        'LO' : '- Selecting models for Binary-lens-Single-source fits with orbital motion'}
+        print(stringmodels[modelcode])
+        completedprocess = subprocess.run([self.bindir+self.modelselectorexe,self.eventname, modelcode], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL)
+        if(completedprocess.returncode != 0):
+            print('! Error in model selection!')
+            self.done = True
+        else:
+            print('  OK')
+            
+    def Finalizer(self):
+        print('- Launching: Finalizer')
+        print('  Making final assessment for this event')
+        completedprocess = subprocess.run([self.bindir+self.finalizerexe,self.eventname], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL)
+        if(completedprocess.returncode != 0):
+            print('! Error in finalization. Maybe there are problems with models')
+            self.done = True
+        else:
+            with open(self.eventname + '/Nature.txt') as f:
+                for line in f.readlines():
+                    print("  " + line,end='')
+                print("  OK")  
+
+    def run(self, event = None):
+        phase =0
+        if(event!= None):
+            self.eventname = os.path.realpath(event)
+        self.done = False        
+        while not(self.done):
+            print("o " + time.asctime())
+            # Check that event directory exists
+            if phase == 0:
+                if(os.path.exists(self.eventname + '/Data')):
+                    print('- Analyzing event: ',self.eventname)
+                    phase = 1
+                else:
+                    print('! Event data for ' + self.eventname + ' not found !')
+                    self.done = True
+            # Launch Reader
+            elif phase == 1:
+                self.Reader()
+                phase = 2
+            # Launch InitCond
+            elif phase == 2:
+                self.InitCond()
+                phase = 3
+            # Launch Finalizer
+            elif phase == self.endphase:
+                self.Finalizer()
+                phase += 1
+            # Conclude analysis
+            elif phase > self.endphase:
+                print("- Analysis of " + self.eventname + " successfully completed!")
+                self.done = True
+            # Launch LevMar for next class
+            elif phase%2 == 1:
+                self.launch_fits(self.modelcodes[phase//2-1]) 
+                phase += 1            
+            # Launch ModelSelector for this class
+            else:
+                self.ModelSelector(self.modelcodes[phase//2-2])
+                phase += 1     
+                
+    def archive_run(self, destination = None):
+        olddir = os.getcwd()
+        os.chdir(self.eventname)
+        if(os.path.exists('LCToFit.txt')):
+            previousrunslist = glob.glob('run-*')
+            previousrunslist.sort()
+            if(destination == None):
+                if(len(previousrunslist)>0):
+                    lastrun = int(previousrunslist[-1].split('-')[-1])
+                else:
+                    lastrun = 0
+                rundir = 'run-' + str(lastrun+1). zfill(4)
+            else:
+                rundir = destination
+            alllist = glob.glob('*')
+            alllist.remove('Data')
+            filelist = list(set(alllist) - set(previousrunslist))
+            os.mkdir(rundir)
+            shutil.copytree('Data',rundir+'/Data')
+            for nam in filelist:
+                shutil.move(nam,rundir)
+        os.chdir(olddir)

RTModel/__init__.py

@@ -0,0 +1,5 @@
+__version__ = "2.0.0"
+__author__ = 'Valerio Bozza'
+__credits__ = 'University of Salerno, Italy'
+
+from .RTModel import RTModel

RTModel/data/TemplateLibrary.txt

@@ -0,0 +1,114 @@
+113						
+0.7	0.5	0.15	3.5	0.01	-0.183	0.016
+0.7	0.5	0.15	3.5	0.01	-0.183	0.086
+0.7	0.5	0.15	3.5	0.01	0.016	0.086
+0.7	0.1	0.	5.38	0.01	-0.066	0.025
+0.7	0.1	0.	5.38	0.01	-0.066	0.305
+0.7	0.1	0.	5.38	0.01	0.025	0.305
+0.7	0.5	0.	2.	0.01	-1.202	-0.113
+0.7	0.5	0.	2.	0.01	-1.202	0.138
+0.7	0.5	0.	2.	0.01	-0.113	0.138
+0.7	0.5	-0.1	1.57	0.01	-0.508	-0.013
+0.7	0.5	-0.1	1.57	0.01	-0.508	0.013
+0.7	0.5	-0.1	1.57	0.01	-0.013	0.013
+0.7	0.5	0.1	1.35	0.01	-0.082	0.075
+0.7	0.5	0.1	1.35	0.01	-0.082	1.121
+0.7	0.5	0.1	1.35	0.01	0.075	1.121
+0.65	0.1	0.52	5.	0.01	-0.291	0.871
+0.65	0.1	0.52	5.	0.01	-0.291	0.879
+0.65	0.1	0.52	5.	0.01	0.871	0.879
+0.7	0.5	-0.35	1.5	0.01	-1.017	-0.975
+0.7	0.5	-0.35	1.5	0.01	-1.017	0.015
+0.7	0.5	-0.35	1.5	0.01	-0.975	0.015
+0.7	0.5	0.8	4.	0.01	0.045	0.552
+0.7	0.5	0.8	4.	0.01	0.045	0.605
+0.7	0.5	0.8	4.	0.01	0.552	0.605
+0.7	0.1	0.1	4.6	0.01	-0.208	-0.007
+0.7	0.1	0.1	4.6	0.01	-0.208	0.149
+0.7	0.1	0.1	4.6	0.01	-0.007	0.149
+0.7	0.5	-0.2	1.4	0.01	-0.516	0.031
+0.7	0.5	-0.2	1.4	0.01	-0.516	1.129
+0.7	0.5	-0.2	1.4	0.01	0.031	1.129
+0.7	0.5	0.3	1.93	0.01	-1.119	0.106
+0.7	0.1	0.25	5.3	0.01	-0.102	0.914
+0.7	0.5	1.12	3.2	0.01	-0.326	-0.095
+0.7	0.1	0.75	4.6	0.01	-0.743	0.032
+0.7	0.1	0.75	4.6	0.01	-0.743	0.54
+0.7	0.1	0.75	4.6	0.01	0.032	0.54
+0.9	0.03	-0.15	2.	0.01	-0.277	0.066
+0.9	0.03	-0.15	2.	0.01	-0.277	0.09
+0.9	0.03	-0.15	2.	0.01	0.066	0.09
+0.9	0.03	-0.1	1.6	0.01	-0.151	-0.122
+0.9	0.03	-0.1	1.6	0.01	-0.151	0.112
+0.9	0.03	-0.1	1.6	0.01	-0.122	0.112
+0.9	0.03	-0.2	2.3	0.01	-0.472	0.01
+0.9	0.03	-0.2	2.3	0.01	-0.472	0.04
+0.9	0.03	-0.2	2.3	0.01	0.01	0.04
+1.5	0.5	0.1	3.2	0.01	-0.459	-0.325
+1.5	0.5	0.1	3.2	0.01	-0.459	0.632
+1.5	0.5	0.1	3.2	0.01	-0.325	0.632
+1.5	0.5	0.	3.	0.01	-0.446	-0.348
+1.5	0.5	0.	3.	0.01	-0.446	0.656
+1.5	0.5	0.	3.	0.01	-0.348	0.656
+0.9	0.03	-0.26	1.65	0.01	-0.222	0.167
+0.9	0.03	-0.26	1.65	0.01	-0.222	0.194
+0.9	0.03	-0.26	1.65	0.01	0.167	0.194
+1.5	0.5	0.3	2.8	0.01	0.314	0.782
+1.5	0.5	0.3	2.8	0.01	0.588	0.782
+0.9	0.03	-0.3	1.8	0.01	-0.285	0.129
+0.9	0.03	-0.3	1.8	0.01	-0.285	0.142
+0.9	0.03	-0.3	1.8	0.01	0.129	0.142
+1.5	0.5	0.3	2.9	0.01	-0.491	0.392
+1.5	0.5	0.3	2.9	0.01	-0.491	0.562
+1.5	0.5	0.3	2.9	0.01	0.392	0.562
+1.5	0.5	0.25	2.75	0.01	0.208	0.738
+1.5	0.5	0.25	2.75	0.01	0.668	0.738
+1.5	0.5	0.35	2.2	0.01	0.06	0.646
+1.5	0.5	0.35	2.2	0.01	0.538	0.646
+1.5	0.5	0.25	2.5	0.01	-0.458	0.096
+1.5	0.5	0.25	2.5	0.01	-0.458	0.587
+1.5	0.5	0.25	2.5	0.01	0.096	0.587
+1.5	0.5	0.25	2.	0.01	-0.041	0.296
+1.5	0.5	0.25	2.	0.01	-0.041	0.472
+1.5	0.5	0.25	2.	0.01	0.296	0.472
+1.5	0.5	0.7	1.5	0.01	-0.08	-0.016
+1.5	0.5	0.7	1.5	0.01	-0.08	0.346
+1.5	0.5	0.7	1.5	0.01	-0.016	0.346
+1.5	0.5	0.45	4.6	0.01	-0.409	-0.049
+1.5	0.5	0.45	4.6	0.01	-0.409	0.277
+1.5	0.5	0.45	4.6	0.01	-0.049	0.277
+1.5	0.5	0.33	3.2	0.01	-0.367	0.652
+0.9	0.1	-0.45	1.65	0.01	-0.466	-0.001
+0.9	0.1	-0.45	1.65	0.01	-0.466	0.372
+0.9	0.1	-0.45	1.65	0.01	-0.001	0.372
+2.2	0.5	0.1	3.	0.01	0.746	1.069
+2.2	0.5	0.1	3.	0.01	0.746	1.189
+2.2	0.5	0.1	3.	0.01	1.069	1.189
+2.2	0.5	0.35	2.9	0.01	1.073	1.142
+2.2	0.5	0.35	2.9	0.01	1.073	1.287
+2.2	0.5	0.35	2.9	0.01	1.142	1.287
+2.2	0.5	0.4	2.8	0.01	0.985	1.169
+2.2	0.5	0.4	2.8	0.01	0.985	1.215
+2.2	0.5	0.4	2.8	0.01	1.169	1.215
+2.2	0.5	1.05	2.	0.01	0.394	0.563
+2.2	0.5	1.05	2.	0.01	0.394	0.657
+2.2	0.5	1.05	2.	0.01	0.563	0.657
+2.2	0.5	0.45	2.7	0.01	-0.675	0.896
+2.2	0.5	0.45	2.7	0.01	-0.675	1.107
+2.2	0.5	0.45	2.7	0.01	0.896	1.107
+2.2	0.5	1.	1.9	0.01	0.3	0.398
+2.2	0.5	1.	1.9	0.01	0.3	0.526
+2.2	0.5	1.	1.9	0.01	0.398	0.526
+2.2	0.5	1.23	1.6	0.01	-0.17	0.018
+2.2	0.5	1.23	1.6	0.01	-0.17	0.053
+2.2	0.5	1.23	1.6	0.01	0.018	0.053
+2.2	0.5	0.35	2.7	0.01	-0.667	0.745
+2.2	0.5	0.35	2.7	0.01	-0.667	1.146
+2.2	0.5	0.35	2.7	0.01	0.745	1.146
+2.2	0.5	0.05	2.95	0.01	-0.645	0.231
+2.2	0.5	0.05	2.95	0.01	-0.645	1.201
+2.2	0.5	0.05	2.95	0.01	0.231	1.201
+2.2	0.5	0.35	3.	0.01	-0.69	1.154
+1.8	0.03	1.	2.3	0.01	-0.084	0.72
+1.8	0.03	1.	2.3	0.01	-0.084	0.889
+1.8	0.03	1.	2.3	0.01	0.72	0.889

RTModel/include/LevMarFit.h

@@ -0,0 +1,81 @@
+// LevMarFit.h
+// Definition of the LevMar and bumper classes for Levenberg-Marquardt fitting
+
+#include <stdio.h>
+#include "bumper.h"
+#include <regex>
+#include <filesystem>
+
+using namespace std;
+using namespace std::filesystem;
+
+#ifndef _LevMarFit
+#define _LevMarFit
+#define __unmanaged
+
+#include <VBBinaryLensingLibrary.h>
+
+class LevMar{
+
+		VBBinaryLensing *VBBL;
+
+		char eventname[512],filename[30],filnum[10],outdir[30],satellitedir[256];
+		char modelcode[16];
+		int error;
+		int flagblending;
+		path exedir;
+
+		double tim0,tm;
+
+		double (VBBinaryLensing::*model)(double *,double);
+		int nps;
+		double *sigmapr,*leftlim,*rightlim;
+
+		int *filter,*satel,nfil,np,OGLE;
+		double *t,*y,*w,*delta,*maxdelta,*Curv,*A,*B,*B0,*Cov,*fb,**Gr,*dFdp,*errs;
+		double *pr,*prn,*sumy, *sumy2, *sumsigma,*sumfy,*sumf,*sumf2,*limbdarks;
+
+		double Tol;
+
+		void (LevMar::*PrintOut)(double *);
+		void (LevMar::*PrintFile)(FILE *,double, bool);
+		
+		bumper *stepchain,*bumperlist,*laststep;
+
+
+	public:
+		LevMar(int,char **);
+		~LevMar();
+
+		void ReadFiles(int,char **);
+		int InitCond(double *presigmapr, double *preleftlim,double *prerightlim);
+		void ReadCurve();
+		void Run();
+		double ChiSquared(double *);
+		void Grad();
+		void Covariance();
+
+		void PrintOutPS(double *);
+		void PrintOutPX(double *);
+		void PrintOutBS(double *);
+		void PrintOutBO(double *);
+		void PrintOutLS(double *);
+		void PrintOutLX(double *);
+		void PrintOutLO(double *);
+		void PrintOutLK(double*);
+
+		void PrintFilePS(FILE *,double, bool);
+		void PrintFilePX(FILE *,double, bool);
+		void PrintFileBS(FILE *,double, bool);
+		void PrintFileBO(FILE *,double, bool);
+		void PrintFileLS(FILE *,double, bool);
+		void PrintFileLX(FILE *,double, bool);
+		void PrintFileLO(FILE *,double, bool);
+		void PrintFileLK(FILE*, double, bool);
+
+};
+
+double Determinant(double *,int);
+void Inverse(double*, double*, int);
+
+#endif