nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/matter/setup_nep.py

@@ -0,0 +1,478 @@
+import os
+import sys
+import numpy as np  # 1.15.0
+from scipy.interpolate import CubicSpline
+
+#nucleardatapy_tk = os.getenv('NUCLEARDATAPY_TK')
+#sys.path.insert(0, nucleardatapy_tk)
+
+import nucleardatapy as nuda
+
+def nep_models():
+    """
+    Return a list of models available in this toolkit and print them all on the prompt.
+
+    :return: The list of models with can be 'Skyrme', 'ESkyrme', 'NLRH', 'DDRH', 'DDRHF'.
+    :rtype: list[str].
+    """
+    #
+    if nuda.env.verb: print("\nEnter pheno_models()")
+    #
+    models = [ 'Skyrme', 'Skyrme2', 'ESkyrme', 'Gogny', 'Fayans', 'NLRH', 'DDRH', 'DDRHF' ]
+    #models = [ 'Skyrme', 'Skyrme2', 'ESkyrme', 'Gogny', 'Fayans', 'NLRH', 'DDRH', 'DDRHF' ]
+    #print('Phenomenological models available in the toolkit:',models)
+    models_lower = [ item.lower() for item in models ]
+    #
+    if nuda.env.verb: print("Exit pheno_models()")
+    #
+    return models, models_lower
+
+def nep_params( model ):
+    """
+    Return a list with the parameterizations available in 
+    this toolkit for a given model and print them all on the prompt.
+
+    :param model: The type of model for which there are parametrizations. They should be chosen among the following options: 'Skyrme', 'NLRH', 'DDRH', 'DDRHF'.
+    :type model: str.
+    :return: The list of parametrizations. \
+    If `models` == 'Skyrme': 'BSK14', \
+    'BSK16', 'BSK17', 'BSK27', 'F-', 'F+', 'F0', 'FPL', 'LNS', 'LNS1', 'LNS5', \
+    'NRAPR', 'RATP', 'SAMI', 'SGII', 'SIII', 'SKGSIGMA', 'SKI2', 'SKI4', 'SKMP', \
+    'SKMS', 'SKO', 'SKOP', 'SKP', 'SKRSIGMA', 'SKX', 'Skz2', 'SLY4', 'SLY5', \
+    'SLY230A', 'SLY230B', 'SV', 'T6', 'T44', 'UNEDF0', 'UNEDF1'. \
+    If `models` == 'ESkyrme': 'BSk22', 'BSk24', 'BSk25', 'BSk26', 'BSk31', 'BSk32', \
+    'BSkG1', 'BSkG2', 'BSkG3'. \
+    If `models` == 'Fayans': 'SLy4', 'SkM*', 'Fy(IVP)', 'Fy(Dr,HDB)', 'Fy(std)', \
+    'SV-min', 'SV-bas', 'SV-K218', 'SV-K226', 'SV-K241', 'SV-mas07', 'SV-mas08', 'SV-mas10',
+    'SV-sym28', 'SV-sym32', 'SV-sym34', 'SV-kap00', 'SV-kap20', 'SV-kap60'.
+    If `models` == 'NLRH': 'NL-SH', 'NL3', 'NL3II', 'PK1', 'PK1R', 'TM1'. \
+    If `models` == 'DDRH': 'DDME1', 'DDME2', 'DDMEd', 'PKDD', 'TW99'. \
+    If `models` == 'DDRHF': 'PKA1', 'PKO1', 'PKO2', 'PKO3'.
+    :rtype: list[str].
+    """
+    #
+    if nuda.env.verb: print("\nEnter nep_params()")
+    #
+    #print('For model:',model)
+    if model.lower() == 'skyrme':
+        params = [ 'BSK14', 'BSK16', 'BSK17', 'BSK27', 'F-', \
+            'F+', 'F0', 'FPL', 'LNS', 'LNS1', 'LNS5', 'NRAPR', 'RATP', \
+            'SAMI', 'SGII', 'SIII', 'SKGSIGMA', 'SKI2', 'SKI4', 'SKMP', \
+            'SKMS', 'SKO', 'SKOP', 'SKP', 'SKRSIGMA', 'SKX', 'Skz2', \
+            'SLY4', 'SLY5', 'SLY230A', 'SLY230B', 'SV', 'T6', 'T44', \
+            'UNEDF0', 'UNEDF1' ]
+    elif model.lower() == 'skyrme2':
+        params = [ 'SLy4', 'SkM*', 'SV-min', 'SV-bas', 'SV-K218', \
+            'SV-K226', 'SV-K241', 'SV-mas07', 'SV-mas08', 'SV-mas10',\
+            'SV-sym28', 'SV-sym32', 'SV-sym34', 'SV-kap00', 'SV-kap20', 'SV-kap60']
+    elif model.lower() == 'eskyrme':
+        params = [ 'BSk22', 'BSk24', 'BSk25', 'BSk26', 'BSk31', 'BSk32', 'BSkG1', 'BSkG2', 'BSkG3' ]
+    elif model.lower() == 'gogny':
+        params = [ 'D1S', 'D1', 'D250', 'D260', 'D280', 'D300' ]
+    elif model.lower() == 'fayans':
+        params = [ 'Fy(IVP)', 'Fy(Dr,HFB)', 'Fy(std)']
+    elif model.lower() == 'nlrh':
+        params = [ 'NL-SH', 'NL3', 'NL3II', 'PK1', 'PK1R', 'TM1' ]
+    elif model.lower() == 'ddrh':
+        params = [ 'DDME1', 'DDME2', 'DDMEd', 'PKDD', 'TW99' ]
+    elif model.lower() == 'ddrhf':
+        params = [ 'PKA1', 'PKO1', 'PKO2', 'PKO3' ]
+    #print('For model:',model)
+    #print('Parameters available in the toolkit:',params)
+    params_lower = [ item.lower() for item in params ]
+    #
+    if nuda.env.verb: print("Exit nep_params()")
+    #
+    return params, params_lower
+
+class setupNEP():
+    """
+    Instantiate the object with results based on phenomenological\
+    interactions and choosen by the toolkit practitioner. \
+    This choice is defined in the variables `model` and `param`.
+
+    If `models` == 'Skyrme', `param` can be: 'BSK14', \
+    'BSK16', 'BSK17', 'BSK27', 'F-', 'F+', 'F0', 'FPL', 'LNS', 'LNS1', 'LNS5', \
+    'NRAPR', 'RATP', 'SAMI', 'SGII', 'SIII', 'SKGSIGMA', 'SKI2', 'SKI4', 'SKMP', \
+    'SKMS', 'SKO', 'SKOP', 'SKP', 'SKRSIGMA', 'SKX', 'Skz2', 'SLY4', 'SLY5', \
+    'SLY230A', 'SLY230B', 'SV', 'T6', 'T44', 'UNEDF0', 'UNEDF1'. 
+
+    If `models` == 'Skyrme2', `param` can be: 'SLy4', 'SkM*', \
+    'SV-min', 'SV-bas', 'SV-K218', 'SV-K226', 'SV-K241', 'SV-mas07', \
+    'SV-mas08', 'SV-mas10','SV-sym28', 'SV-sym32', 'SV-sym34', 'SV-kap00', \
+    'SV-kap20', 'SV-kap60'.
+
+    If `models` == 'ESkyrme', `param` can be: 'BSk22', 'BSk24', 'BSk25', \
+    'BSk26', 'BSk31', 'BSk32', 'BSkG1', 'BSkG2', 'BSkG3'.
+
+    If `models` == 'Fayans', `param` can be: 'Fy(IVP)', 'Fy(Dr,HDB)', 'Fy(std)'.
+
+    If `models` == 'NLRH', `param` can be: 'NL-SH', 'NL3', 'NL3II', 'PK1', 'PK1R', 'TM1'. 
+
+    If `models` == 'DDRH', `param` can be: 'DDME1', 'DDME2', 'DDMEd', 'PKDD', 'TW99'. 
+
+    If `models` == 'DDRHF', `param` can be: 'PKA1', 'PKO1', 'PKO2', 'PKO3'. 
+    
+    :param model: Fix the name of model: 'Skyrme', 'NLRH', \
+    'DDRH', 'DDRHF'. Default value: 'Skyrme'.
+    :type model: str, optional. 
+    :param param: Fix the parameterization associated to model. \
+    Default value: 'SLY5'.
+    :type param: str, optional. 
+
+    **Attributes:**
+    """
+    #
+    def __init__( self, model = 'Skyrme', param = 'SLY5' ):
+        #
+        if nuda.env.verb: print("\nEnter setupNEP()")
+        #
+        #: Attribute model.
+        self.model = model
+        if nuda.env.verb: print("model:",model)
+        #print("-> model:",model)
+        #: Attribute param.
+        self.param = param
+        if nuda.env.verb: print("param:",param)
+        #print("-> param:",param)
+        #
+        self = setupNEP.init_self( self )
+        #
+        models, models_lower = nep_models( )
+        #
+        if model.lower() not in models_lower:
+            print('setup_nep: The model name ',model,' is not in the list of models.')
+            print('setup_nep: list of models:',models)
+            print('setup_nep: -- Exit the code --')
+            exit()
+        #
+        params, params_lower = nep_params( model = model )
+        #
+        if param.lower() not in params_lower:
+            print('setup_nep: The param set ',param,' is not in the list of param.')
+            print('setup_nep: list of param:',params)
+            print('setup_nep: -- Exit the code --')
+            exit()
+        #
+        if model.lower() == 'skyrme':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPSkyrme.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            #print('Reads file:',file_in)
+            #: Attribute providing the full reference to the paper to be citted.
+            #self.ref = ''
+            #: Attribute providing the label the data is references for figures.
+            self.label = 'SKY-'+param
+            #: Attribute providing additional notes about the data.
+            self.note = "write here notes about this EOS."
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype = str )
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            #print('name:',name)
+            #
+            #nsat, Esat, Ksat, Qsat, Esym, Lsym, Ksym, msat \
+            #   = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8), comments='#', unpack = True, dtype = float )
+            #kappas = 1./msat - 1.0
+            #kappav = np.zeros( kappas.size )
+            #Zsat = kappav.copy()
+            #Qsym = kappav.copy()
+            #Zsym = kappav.copy()
+            nsat, Esat, Ksat, Qsat, Zsat, Esym, Lsym, Ksym, Qsym, Zsym, \
+                msat, kappas, kappav = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8,9,10,11,12,13), comments='#', unpack = True )
+            kappasym = kappas - kappav
+            Dmsat = -2*kappasym/( (1+kappas)**2-kappasym**2)
+            #
+            #print('param:',param)
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.nsat = nsat[ind]; self.Esat = Esat[ind]; self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind]; self.Zsat = Zsat[ind]; 
+                self.Esym = Esym[ind]; self.Lsym = Lsym[ind]; self.Ksym = Ksym[ind]; self.Qsym = Qsym[ind]; self.Zsym = Zsym[ind];
+                self.msat = msat[ind]; self.kappas = kappas[ind]; self.kappav = kappav[ind];
+                self.kappasym = kappasym[ind]; self.Dmsat = Dmsat[ind]
+            else:
+                self.nep = False
+        #
+        elif model.lower() == 'skyrme2':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPSkyrme2.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            self.label = 'Skyrme2-'+param
+            self.note = "write here notes about this EOS."
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype=str )
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            nsat, Esat, Ksat, Qsat, msat, Esym, Lsym, kappav \
+                = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8), comments='#', unpack = True )
+            kappas = 1.0/msat - 1.0
+            kappasym = kappas - kappav
+            Dmsat = -2*kappasym/( (1+kappas)**2-kappasym**2)
+            #
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.nsat = nsat[ind]; self.Esat = Esat[ind]; self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind];
+                self.Esym = Esym[ind]; self.Lsym = Lsym[ind];
+                self.msat = msat[ind]; self.kappas = kappas[ind]; self.kappav = kappav[ind];
+                self.kappasym = kappasym[ind]; self.Dmsat = Dmsat[ind]
+            else:
+                self.nep = False
+            #
+        elif model.lower() == 'eskyrme':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPESkyrme.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            #: Attribute providing the full reference to the paper to be citted.
+            #self.ref = ''
+            #: Attribute providing the label the data is references for figures.
+            self.label = 'ESKY-'+param
+            #: Attribute providing additional notes about the data.
+            self.note = "write here notes about this EOS."
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype=str )
+            nsat, Esat, Ksat, Qsat, Esym, Lsym, Ksym, msat = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8), comments='#', unpack = True )
+            kappas = 1.0/msat - 1.0
+            kappav = np.zeros( kappas.size )
+            Zsat = kappav.copy()
+            Qsym = kappav.copy()
+            Zsym = kappav.copy()
+            kappasym = kappas - kappav
+            Dmsat = -2*kappasym/( (1+kappas)**2-kappasym**2)
+            #
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.nsat = nsat[ind]; self.Esat = Esat[ind]; self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind]; self.Zsat = None
+                self.Esym = Esym[ind]; self.Lsym = Lsym[ind]; self.Ksym = Ksym[ind]; self.Qsym = None; self.Zsym = None
+                self.msat = msat[ind]; self.kappas = kappas[ind]; self.kappav = None;
+                self.kappasym = None; self.Dmsat = None
+            else:
+                self.nep = False
+            #
+        elif model.lower() == 'gogny':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPGogny.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            self.label = 'Gogny-'+param
+            self.note = "write here notes about this EOS."
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype=str )
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            Ksat, Qsat = np.loadtxt( file_in, usecols=(1,2), comments='#', unpack = True )
+            #
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind];
+                self.nsat = None; self.Esat = None; self.Zsat = None
+                self.Esym = None; self.Lsym = None; self.Ksym = None; self.Qsym = None; self.Zsym = None
+                self.msat = None; self.kappas = None; self.kappav = None;
+                self.kappasym = None; self.Dmsat = None
+            else:
+                self.nep = False
+#            pass
+            #
+        elif model.lower() == 'fayans':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPFayans.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            self.label = 'Fayans-'+param
+            self.note = "write here notes about this EOS."
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype=str )
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            nsat, Esat, Ksat, Qsat, msat, Esym, Lsym, kappav \
+                = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8), comments='#', unpack = True )
+            kappas = 1.0/msat - 1.0
+            kappasym = kappas - kappav
+            Dmsat = -2*kappasym/( (1+kappas)**2-kappasym**2)
+            #
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.nsat = nsat[ind]; self.Esat = Esat[ind]; self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind];
+                self.Esym = Esym[ind]; self.Lsym = Lsym[ind];
+                self.msat = msat[ind]; self.kappas = kappas[ind]; self.kappav = kappav[ind];
+                self.kappasym = kappasym[ind]; self.Dmsat = Dmsat[ind]
+            else:
+                self.nep = False
+            #
+        elif model.lower() == 'nlrh':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPnlrh.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            #self.ref = ''
+            self.label = 'NLRH-'+param
+            self.note = "write here notes about this EOS."
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype=str )
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            nsat, Esat, Ksat, Qsat, Zsat, Esym, Lsym, Ksym, Qsym, Zsym, \
+                msat, kappas, kappav = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8,9,10,11,12,13), comments='#', unpack = True )
+            kappasym = kappas - kappav
+            Dmsat = -2*kappasym/( (1+kappas)**2-kappasym**2)
+            #
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.nsat = nsat[ind]; self.Esat = Esat[ind]; self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind]; self.Zsat = Zsat[ind]; 
+                self.Esym = Esym[ind]; self.Lsym = Lsym[ind]; self.Ksym = Ksym[ind]; self.Qsym = Qsym[ind]; self.Zsym = Zsym[ind];
+                self.msat = msat[ind]; self.kappas = kappas[ind]; self.kappav = kappav[ind];
+                self.kappasym = kappasym[ind]; self.Dmsat = Dmsat[ind]
+            else:
+                self.nep = False
+            #
+        elif model.lower() == 'ddrh':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPddrh.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            #self.ref = ''
+            self.label = 'DDRH-'+param
+            self.note = "write here notes about this EOS."
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype=str )
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            nsat, Esat, Ksat, Qsat, Zsat, Esym, Lsym, Ksym, Qsym, Zsym, \
+                msat, kappas, kappav = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8,9,10,11,12,13), comments='#', unpack = True )
+            kappasym = kappas - kappav
+            Dmsat = -2*kappasym/( (1+kappas)**2-kappasym**2)
+            #
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.nsat = nsat[ind]; self.Esat = Esat[ind]; self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind]; self.Zsat = Zsat[ind];
+                self.Esym = Esym[ind]; self.Lsym = Lsym[ind]; self.Ksym = Ksym[ind]; self.Qsym = Qsym[ind]; self.Zsym = Zsym[ind];
+                self.msat = msat[ind]; self.kappas = kappas[ind]; self.kappav = kappav[ind];
+                self.kappasym = kappasym[ind]; self.Dmsat = Dmsat[ind]
+            else:
+                self.nep = False
+            #
+        elif model.lower() == 'ddrhf':
+            #
+            file_in = os.path.join(nuda.param.path_data,'matter/nep/NEPddrhf.dat')
+            if nuda.env.verb: print('Reads file:',file_in)
+            #self.ref = ''
+            self.label = 'DDRHF-'+param
+            self.note = "write here notes about this EOS."
+            #name = np.loadtxt( file_in, usecols=(0), comments='#', unpack = True, dtype=str )
+            name = []
+            with open(file_in,"r") as file:
+                for line in file:
+                    if '#' in line:
+                        continue
+                    name.append( line.split()[0] )
+            name = np.array( name, dtype = str )
+            nsat, Esat, Ksat, Qsat, Zsat, Esym, Lsym, Ksym, Qsym, Zsym, \
+                msat, kappas, kappav = np.loadtxt( file_in, usecols=(1,2,3,4,5,6,7,8,9,10,11,12,13), comments='#', unpack = True )
+            kappasym = kappas - kappav
+            Dmsat = -2*kappasym/( (1+kappas)**2-kappasym**2)
+            #
+            if param in name:
+                self.nep = True
+                ind = np.where( name == param )[0][0]
+                self.nsat = nsat[ind]; self.Esat = Esat[ind]; self.Ksat = Ksat[ind]; self.Qsat = Qsat[ind]; self.Zsat = Zsat[ind];
+                self.Esym = Esym[ind]; self.Lsym = Lsym[ind]; self.Ksym = Ksym[ind]; self.Qsym = Qsym[ind]; self.Zsym = Zsym[ind];
+                self.msat = msat[ind]; self.kappas = kappas[ind]; self.kappav = kappav[ind];
+                self.kappasym = kappasym[ind]; self.Dmsat = Dmsat[ind]
+            else:
+                self.nep = False
+            #
+        self.den_unit = 'fm$^{-3}$'
+        self.kfn_unit = 'fm$^{-1}$'
+        self.e2a_unit = 'MeV'
+        self.pre_unit = 'MeV fm$^{-3}$'
+        self.gap_unit = 'MeV'
+        #
+        if nuda.env.verb: print("Exit SetupNEP()")
+    #
+    def print_outputs( self ):
+        """
+        Method which print outputs on terminal's screen.
+        """
+        print("")
+        #
+        if nuda.env.verb: print("Enter print_outputs()")
+        #
+        print("- Print output:")
+        print("   model:",self.model)
+        print("   param:",self.param)
+        #print("   ref:",self.ref)
+        print("   label:",self.label)
+        #print("   note:",self.note)
+        if any(self.sm_den): print(f"   sm_den: {np.round(self.sm_den,2)} in {self.den_unit}")
+        if any(self.sm_kfn): print(f"   sm_kfn: {np.round(self.sm_kfn,2)} in {self.kfn_unit}")
+        if any(self.sm_e2a): print(f"   sm_e2a: {np.round(self.sm_e2a,2)} in {self.e2a_unit}")
+        if any(self.nm_den): print(f"   nm_den: {np.round(self.nm_den,2)} in {self.den_unit}")
+        if any(self.nm_kfn): print(f"   nm_kfn: {np.round(self.nm_kfn,2)} in {self.kfn_unit}")
+        if any(self.nm_e2a): print(f"   nm_e2a: {np.round(self.nm_e2a,2)} in {self.e2a_unit}")
+        if any(self.nm_gap): print(f"   nm_gap: {np.round(self.nm_gap,2)} in {self.gap_unit}")
+        print(' NEP:')
+        if self.nep:
+            print(' sat:',self.Esat,self.nsat,self.Ksat,self.Qsat,self.Zsat)
+            print(' sym:',self.Esym,self.Lsym,self.Ksym,self.Qsym,self.Zsym)
+            print(' ms:',self.msat,self.kappas,self.kappav)
+        #
+        if nuda.env.verb: print("Exit print_outputs()")
+        #
+    def init_self( self ):
+        """
+        Initialize variables in self.
+        """
+        #
+        if nuda.env.verb: print("Enter init_self()")
+        #
+        #: Attribute the plot linestyle.
+        self.linestyle = 'solid'
+        #: Attribute providing the full reference to the paper to be citted.
+        self.ref = ''
+        #: Attribute providing additional notes about the data.
+        self.note = ''
+        #: Attribute the plot label data.
+        self.label = ''
+        #: Attribute the plot marker.
+        self.marker = None
+        #: Attribute the NEP.
+        self.Esat = None; self.nsat = None; self.Ksat = None; self.Qsat = None; self.Zsat = None;
+        self.Esym = None; self.Lsym = None; self.Ksym = None; self.Qsym = None; self.Zsym = None;
+        self.msat = None; self.kappas = None; self.kappav = None; 
+        self.kappasym = None; self.Dmsat = None;
+        #
+        if nuda.env.verb: print("Exit init_self()")
+        #
+        return self

nucleardatapy/matter/setup_nep_dist.py

@@ -0,0 +1,187 @@
+import os
+import sys
+import numpy as np  # 1.15.0
+from scipy.interpolate import CubicSpline
+
+import nucleardatapy as nuda
+
+class setupNEPDist():
+    """
+    Instantiate the object with results based on phenomenological\
+    interactions and choosen by the toolkit practitioner. \
+    This choice is defined in the variables `model` and `param`.
+
+    :param model: Fix the name of model: 'Skyrme', 'NLRH', \
+    'DDRH', 'DDRHF'. Default value: 'Skyrme'.
+    :type model: str, optional. 
+    :param param: Fix the parameterization associated to model. \
+    Default value: 'SLY5'.
+    :type param: str, optional. 
+
+    **Attributes:**
+    """
+    #
+    def __init__( self, models = [ 'Skyrme' ] ):
+        #
+        if nuda.env.verb: print("\nEnter setupNEPDist()")
+        #
+        #: Attribute model.
+        self.models = models
+        if nuda.env.verb: print("models:",models)
+        #
+        self = setupNEPDist.init_self( self )
+        #
+        models, models_lower = nuda.matter.nep_models( )
+        #
+        for model in models:
+            #
+            if model.lower() not in models_lower:
+                print('The model name ',model,' is not in the list of models.')
+                print('list of models:',models)
+                print('-- Exit the code --')
+                exit()
+            #
+            dist = nuda.matter.setupNEPModelDist( model )
+            #        
+            if dist.nep:
+                if dist.params: self.params.extend( dist.params ); 
+                if dist.Esat: self.Esat.extend( dist.Esat ); 
+                if dist.nsat: self.nsat.extend( dist.nsat ); 
+                if dist.Ksat: self.Ksat.extend( dist.Ksat ); 
+                if dist.Qsat: self.Qsat.extend( dist.Qsat ); 
+                if dist.Zsat: self.Zsat.extend( dist.Zsat ); 
+                if dist.Esym: self.Esym.extend( dist.Esym ); 
+                if dist.Lsym: self.Lsym.extend( dist.Lsym ); 
+                if dist.Ksym: self.Ksym.extend( dist.Ksym ); 
+                if dist.Qsym: self.Qsym.extend( dist.Qsym ); 
+                if dist.Zsym: self.Zsym.extend( dist.Zsym ); 
+                if dist.kappas: self.kappas.extend( dist.kappas ); 
+                if dist.kappav: self.kappav.extend( dist.kappav ); 
+                if dist.kappasym: self.kappasym.extend( dist.kappasym ); 
+                if dist.msat: self.msat.extend( dist.msat ); 
+                if dist.Dmsat: self.Dmsat.extend( dist.Dmsat ); 
+        self.Esat = np.array( self.Esat, dtype = float )
+        self.nsat = np.array( self.nsat, dtype = float )
+        self.Ksat = np.array( self.Ksat, dtype = float )
+        self.Qsat = np.array( self.Qsat, dtype = float )
+        self.Zsat = np.array( self.Zsat, dtype = float )
+        self.Esym = np.array( self.Esym, dtype = float )
+        self.Lsym = np.array( self.Lsym, dtype = float )
+        self.Ksym = np.array( self.Ksym, dtype = float )
+        self.Qsym = np.array( self.Qsym, dtype = float )
+        self.Zsym = np.array( self.Zsym, dtype = float )
+        self.kappas = np.array( self.kappas, dtype = float )
+        self.kappav = np.array( self.kappav, dtype = float )
+        self.kappasym = np.array( self.kappasym, dtype = float )
+        self.msat = np.array( self.msat, dtype = float )
+        self.Dmsat = np.array( self.Dmsat, dtype = float )
+        #
+        #  Compute statistical properties: centroids and standard distributions
+        #
+        self.Esat_mean = np.mean(self.Esat); self.Esat_std = np.std(self.Esat)
+        self.nsat_mean = np.mean(self.nsat); self.nsat_std = np.std(self.nsat)
+        self.Ksat_mean = np.mean(self.Ksat); self.Ksat_std = np.std(self.Ksat)
+        self.Qsat_mean = np.mean(self.Qsat); self.Qsat_std = np.std(self.Qsat)
+        self.Zsat_mean = np.mean(self.Zsat); self.Zsat_std = np.std(self.Zsat)
+        self.Esym_mean = np.mean(self.Esym); self.Esym_std = np.std(self.Esym)
+        self.Lsym_mean = np.mean(self.Lsym); self.Lsym_std = np.std(self.Lsym)
+        self.Ksym_mean = np.mean(self.Ksym); self.Ksym_std = np.std(self.Ksym)
+        self.Qsym_mean = np.mean(self.Qsym); self.Qsym_std = np.std(self.Qsym)
+        self.Zsym_mean = np.mean(self.Zsym); self.Zsym_std = np.std(self.Zsym)
+        self.kappas_mean = np.mean(self.kappas); self.kappas_std = np.std(self.kappas)
+        self.kappav_mean = np.mean(self.kappav); self.kappav_std = np.std(self.kappav)
+        self.kappasym_mean = np.mean(self.kappasym); self.kappasym_std = np.std(self.kappasym)
+        self.msat_mean = np.mean(self.msat);   self.msat_std = np.std(self.msat)
+        self.Dmsat_mean = np.mean(self.Dmsat); self.Dmsat_std = np.std(self.Dmsat)
+        #
+        if nuda.env.verb: print("Exit SetupNEPDist()")
+    #
+    def print_outputs( self ):
+        """
+        Method which print outputs on terminal's screen.
+        """
+        print("")
+        #
+        if nuda.env.verb: print("Enter print_outputs()")
+        #
+        print("- Print output:")
+        print("   models:",self.models)
+        print(' NEP:')
+        if any(self.Esat): print(f"   Esat: {np.round(self.Esat,2)}")
+        if any(self.nsat): print(f"   nsat: {np.round(self.nsat,3)}")
+        if any(self.Ksat): print(f"   Ksat: {np.round(self.Ksat,2)}")
+        if any(self.Qsat): print(f"   Qsat: {np.round(self.Qsat,2)}")
+        if any(self.Zsat): print(f"   Zsat: {np.round(self.Zsat,2)}")
+        if any(self.Esym): print(f"   Esym: {np.round(self.Esym,2)}")
+        if any(self.Lsym): print(f"   Lsym: {np.round(self.Lsym,2)}")
+        if any(self.Ksym): print(f"   Ksym: {np.round(self.Ksym,2)}")
+        if any(self.Qsym): print(f"   Qsym: {np.round(self.Qsym,2)}")
+        if any(self.Zsym): print(f"   Zsym: {np.round(self.Zsym,2)}")
+        if any(self.msat): print(f"   msat: {np.round(self.msat,2)}")
+        if any(self.Dmsat): print(f"  Dmsat: {np.round(self.Dmsat,2)}")
+        #
+        print(f" Esat_mean: {self.Esat_mean:.2f} +- {self.Esat_std:.2f}")
+        print(f" nsat_mean: {self.nsat_mean:.2f} +- {self.nsat_std:.2f}")
+        print(f" Ksat_mean: {self.Ksat_mean:.1f} +- {self.Ksat_std:.1f}")
+        print(f" Qsat_mean: {self.Qsat_mean:.0f} +- {self.Qsat_std:.0f}")
+        print(f" Zsat_mean: {self.Zsat_mean:.0f} +- {self.Zsat_std:.0f}")
+        print(f" Esym_mean: {self.Esym_mean:.2f} +- {self.Esym_std:.2f}")
+        print(f" Lsym_mean: {self.Lsym_mean:.1f} +- {self.Lsym_std:.1f}")
+        print(f" Ksym_mean: {self.Ksym_mean:.0f} +- {self.Ksym_std:.0f}")
+        print(f" Qsym_mean: {self.Qsym_mean:.0f} +- {self.Qsym_std:.0f}")
+        print(f" Zsym_mean: {self.Zsym_mean:.0f} +- {self.Zsym_std:.0f}")
+        print(f" kappas_mean: {self.kappas_mean:.2f} +- {self.kappas_std:.2f}")
+        print(f" kappav_mean: {self.kappav_mean:.2f} +- {self.kappav_std:.2f}")
+        print(f" kappasym_mean: {self.kappasym_mean:.2f} +- {self.kappasym_std:.2f}")
+        print(f" msat_mean: {self.msat_mean:.2f} +- {self.msat_std:.2f}")
+        print(f" Dmsat_mean: {self.Dmsat_mean:.2f} +- {self.Dmsat_std:.2f}")
+        #
+        if nuda.env.verb: print("Exit print_outputs()")
+        #
+    def print_latex( self ):
+        """
+        Method which print outputs on terminal's screen.
+        """
+        print("")
+        #
+        if nuda.env.verb: print("Enter print_latex()")
+        #
+        print("- Print latex:")
+        print("   models:", len(self.models), self.models, len(self.params))
+        print(' table:')
+        print(r' NEP & $E_{\sat}$ & $n_{\sat}$ & $K_{\sat}$ & $Q_{\sat}$ & $Z_{\sat}$ &',
+            r' $E_{sym}$ & $L_{sym}$ & $K_{sym}$ & $Q_{sym}$ & $Z_{sym}$ &',
+            r' $m^*_{sat}/m$ & $\Delta m^*_{sat}/m$ \\\\')
+        print(r' & MeV & fm$^{-3}n$ & MeV & MeV & MeV &',
+            r' MeV & MeV & MeV & MeV & MeV &',
+            r'  &  \\\\')
+        print(rf' centroid & {self.Esat_mean:.2f} & {self.nsat_mean:.3f} &',
+            rf' {self.Ksat_mean:.1f} & {self.Qsat_mean:.0f} & {self.Zsat_mean:.0f} &',
+            rf' {self.Esym_mean:.2f} & {self.Lsym_mean:.1f} & {self.Ksym_mean:.0f} &',
+            rf' {self.Qsym_mean:.0f} & {self.Zsym_mean:.0f} & {self.msat_mean:.2f} &',
+            rf' {self.Dmsat_mean:.3f} \\\\')
+        print(rf' std.dev. & {self.Esat_std:.2f} & {self.nsat_std:.3f} &',
+            rf' {self.Ksat_std:.1f} & {self.Qsat_std:.0f} & {self.Zsat_std:.0f} &',
+            rf' {self.Esym_std:.2f} & {self.Lsym_std:.1f} & {self.Ksym_std:.0f} &',
+            rf' {self.Qsym_std:.0f} & {self.Zsym_std:.0f} & {self.msat_std:.2f} &',
+            rf' {self.Dmsat_std:.3f} \\\\')
+        #
+        if nuda.env.verb: print("Exit print_latex()")
+        #
+    def init_self( self ):
+        """
+        Initialize variables in self.
+        """
+        #
+        if nuda.env.verb: print("Enter init_self()")
+        #
+        #: Attribute the NEP.
+        self.params = []
+        self.nsat = []; self.Esat = []; self.Ksat = []; self.Qsat = []; self.Zsat = []
+        self.Esym = []; self.Lsym = []; self.Ksym = []; self.Qsym = []; self.Zsym = []
+        self.kappas = []; self.kappav = []; self.kappasym = []
+        self.msat = []; self.Dmsat = []; 
+        #
+        if nuda.env.verb: print("Exit init_self()")
+        #
+        return self