CG_Iterations Namespace Reference

Functions
def	Parse_Nektar_Output

Variables
string	Mesh = 'cyl-small.xml'
string	Max_N_Z = 'output_72.txt'
string	Conditions_File = 'conditions_80.xml'
string	Scheme = 'IterativeStaticCond'
list	Consider_Modes = [80]
int	Num_Constants = 2
	Compare_Serial = False
	Parallel = True
string	Parallelisation = 'Hybrid_Socket'
	Compare_Parallel = True
int	Num_Node = 1
string	output_path = 'Output/Figures/Max_Min_CG/'
string	cmd_string_clear = 'rm -r Output/Figures/Max_Min_CG/ \n'
tuple	process = Popen([cmd_string_clear],shell=True, stdout=PIPE, stdin=PIPE)
list	Pressure = []
list	Velocity_1 = []
list	Velocity_2 = []
list	Velocity_3 = []
list	Pressure_Key = []
list	Velocity_1_Key = []
list	Velocity_2_Key = []
list	Velocity_3_Key = []
list	Max_Pressure_CG = []
list	Max_Velocity_1_CG = []
list	Max_Velocity_2_CG = []
list	Max_Velocity_3_CG = []
list	Min_Pressure_CG = []
list	Min_Velocity_1_CG = []
list	Min_Velocity_2_CG = []
list	Min_Velocity_3_CG = []
int	Skip_Steps = 20
tuple	key_P = map(int, Pressure[i].keys())
tuple	key_V_1 = map(int, Velocity_1[i].keys())
tuple	key_V_2 = map(int, Velocity_2[i].keys())
tuple	key_V_3 = map(int, Velocity_3[i].keys())
list	CG_Min_Pressure = []
list	CG_Max_Pressure = []
list	CG_Min_Velocity_1 = []
list	CG_Max_Velocity_1 = []
list	CG_Min_Velocity_2 = []
list	CG_Max_Velocity_2 = []
list	CG_Min_Velocity_3 = []
list	CG_Max_Velocity_3 = []
string	savefile = 'Output/Figures/Max_Min_CG/Min_CG'
string	title = 'Minimum CG Iterations: '
tuple	max_modes = len(Pressure[len(Timing_Files) - 1])
list	modes = []
list	data_mean = [[], [], [], []]
list	data_std_dev = [[], [], [], []]
list	data_var = [[], [], [], []]
list	Data_P = []
list	Data_V_1 = []
list	Data_V_2 = []
list	Data_V_3 = []
list	Settling_Deviation_P = []
list	Settling_Deviation_V_1 = []
list	Settling_Deviation_V_2 = []
list	Settling_Deviation_V_3 = []
list	Settling_Time = []
	Error = True
list	Data_2_P = []
list	Data_2_V_1 = []
list	Data_2_V_2 = []
list	Data_2_V_3 = []

Function Documentation

def CG_Iterations.Parse_Nektar_Output	(		Input_Filename,
			Skip_Steps
	)

Definition at line 31 of file CG_Iterations.py.

Referenced by parallel.Run_Parallel_Comparison().

def Parse_Nektar_Output(Input_Filename, Skip_Steps):

    # Open file to be pharsed
    f = open(Input_Filename, "r")
    
    # List for outputs
    Times = []

    # We use dictionaries to generate look up tables of the CG information
    Pressure = {}
    Velocity_1 = {}
    Velocity_2 = {}
    Velocity_3 = {}

    # Count timestep we're on
    count = 0

    # Pharse desired data
    for line in f:
        a = line.split()
        for i in range(0, len(a)):
            
            # Iterate the timestep counter every time the file hits 'Pressure'
            if (a[i] == 'Pressure'):
                Velocity = False

                # Var is iterated depending on pressure or velocity being read in
                var = 1
                count += 1
                continue
            
            # Skip these time steps as they contain all the start up computations
            if (count <= Skip_Steps):
                continue

            if (a[i] == 'Velocity'):
                Velocity = True
                var = 1
                continue
            
            # Record the plane number
            if (a[i] == 'Plane'):
                plane = a[i + 1]

                if Velocity is True and plane == '0':
                    var += 1

                continue

            # Append each value of CG to the dictionaries
            if (a[i] == 'CG'):
                if (var == 1):
                    if plane in Pressure.keys():
                        Pressure[plane].append(int(a[i + 4]))
                        continue
                    else:
                        Pressure[plane] = [int(a[i + 4])]

                if (var == 2):
                    if plane in Velocity_1.keys():
                        Velocity_1[plane].append(int(a[i + 4]))
                        continue
                    else:
                        Velocity_1[plane] = [int(a[i + 4])]

                if (var == 3):
                    if plane in Velocity_2.keys():
                        Velocity_2[plane].append(int(a[i + 4]))
                        continue
                    else:
                        Velocity_2[plane] = [int(a[i + 4])]

                if (var == 4):
                    if plane in Velocity_3.keys():
                        Velocity_3[plane].append(int(a[i + 4]))
                        continue
                    else:
                        Velocity_3[plane] = [int(a[i + 4])]

    
    
    # Return the dictionaries of CG iteration
    return(Pressure, Velocity_1, Velocity_2, Velocity_3)

#------------------------------------
# Functions end here
#------------------------------------

#------------------------------------
# Program begins here
#------------------------------------