np.seterr(divide='ignore') if __name__ == '__main__': logger = logging.getLogger('LOG') logger.setLevel(logging.DEBUG) fh = logging.FileHandler('LOGS_SolAnVsNum.log') fh.setLevel(logging.DEBUG) logger.addHandler(fh) formatter = logging.Formatter( '%(asctime)s - %(name)s - %(levelname)s - %(message)s') fh.setFormatter(formatter) logger.addHandler(fh) pr = gp.get_parameters() ################################################################## # Load txt file of temperature ################################################################## logger.debug("Loading txt file of temperature") try: matriz = np.loadtxt('Results/T_.txt') except Exception as e: logger.critical("No se pudo cargar el archico: ") sys.exit() size = pr["lattice_x"] # size ydata = np.zeros([size]) for i in range(size): ydata[i] = matriz[int(pr["lattice_y"] // 2)][i] xdata = np.linspace(0, size - 1, size)
#! /usr/bin/env python # -*- coding: utf-8 -*- # Authors: Benjamín S. Noyola García, Suemi Rodriguez Romo import local as loc from numba import cuda, float64, float32 import numpy as np import get_parser_parameters as gp pr = gp.get_parameters() # pr is a dictionary of parameters T_h = pr["T_h"] T_c = pr["T_c"] # links to local.py program ###################### velocity field ####################################### getf = cuda.jit('void (f8[:,:,:], f8[::1], i8, i8)', device=True)(loc.getf) getux = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getux) getuy = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getuy) getFx = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getFx) getFy = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getFy) getden = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getden) getf_l = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getf_l) getf_2l = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getf_2l) getT = cuda.jit('void (f8[:,:], f8[::1], i8, i8)', device=True)(loc.getT) f2m = cuda.jit('void (f8[::1], f8[::1])', device=True)(loc.f2m) calc_fl_PCMloc = cuda.jit('void (f8[::1], f8[::1])', device=True)(loc.calc_fl_PCMloc) calc_noruloc = cuda.jit('void (f8[::1], f8[::1], f8[::1])', device=True)(loc.calc_noruloc) calc_m_eqloc = cuda.jit( 'void (f8[::1], f8[::1], f8[::1], f8[::1], f8[::1], f8[::1], f8[::1])', device=True)(loc.calc_m_eqloc)