Changeset 802 for codes/icosagcm/devel/Python/src/unstructured.pyx
- Timestamp:
- 01/21/19 11:34:04 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/src/unstructured.pyx
r792 r802 9 9 cimport numpy as np 10 10 import numpy as np 11 12 import getargs 13 log_master, log_world = getargs.getLogger(__name__) 14 INFO, DEBUG, ERROR = log_master.info, log_master.debug, log_world.error 15 INFO_ALL, DEBUG_ALL = log_world.info, log_world.debug 11 16 12 17 #-------------- choose precision of kernel computations ------------# … … 58 63 59 64 try: 60 61 65 kernels = wrap.SharedLib(vars(ker), libicosa, check_args=check_args) 66 setvar, setvars, getvar, getvars = kernels.setvar, kernels.setvars, kernels.getvar, kernels.getvars 62 67 except OSError: 63 print """ 64 Unable to load shared library 'libicosa.so' ! 65 """ 66 raise 68 ERROR("Unable to load shared library 'libicosa.so' !") 69 raise 67 70 68 71 # providing a full prototype enables type-checking when calling … … 90 93 91 94 kernels.addvars( 92 c_bool,'hydrostatic','debug_hevi_solver', 95 c_bool,'hydrostatic','debug_hevi_solver', 'is_mpi_master', 93 96 c_int,'llm','nqdyn','primal_num','max_primal_deg', 94 97 'dual_num','max_dual_deg','edge_num','max_trisk_deg', 95 'caldyn_thermo','caldyn_eta','nb_threads', 'dynamico_mpi_rank',98 'caldyn_thermo','caldyn_eta','nb_threads', 96 99 c_double,'elapsed', 97 100 c_num, 'g', 'ptop', 'cpp', 'cppv', … … 119 122 return data 120 123 121 cdef check_ptr(name, num_ptr p, ndarray data):122 if p != ptr(data) : print name, 'p <> ptr(data) !!'123 124 124 cdef class Caldyn_step: 125 125 # number of time steps to do at each invocation of advance() … … 300 300 Ai, Av, fv, le_de, Riv2, wee): 301 301 302 print 'Types of Ai, Av, fv, le_de, Riv2, wee : ', Ai.dtype,Av.dtype,fv.dtype,le_de.dtype,Riv2.dtype,wee.dtype 302 DEBUG('Types of Ai, Av, fv, le_de, Riv2, wee : %s' % 303 ((Ai.dtype,Av.dtype,fv.dtype,le_de.dtype,Riv2.dtype,wee.dtype),) ) 303 304 for var,varname in zip((Ai,Av,fv,le_de,Riv2,wee), ('Ai','Av','fv','le_de','Riv2','wee')): 304 305 assert var.dtype == np.float64, '%s must be double precision'%varname … … 308 309 (llm, nqdyn, edge_num, primal_num, dual_num, 309 310 max_trisk_deg, max_primal_deg, max_dual_deg) ) 310 print('init_mesh ...')311 INFO('init_mesh ...') 311 312 ker.dynamico_init_mesh(primal_nb,primal_edge,primal_ne, 312 313 dual_nb,dual_edge,dual_ne,dual_vertex, 313 314 left,right,down,up,trisk_deg,trisk) 314 print('...done')315 print('init_metric ...')315 INFO('...done') 316 INFO('init_metric ...') 316 317 ker.dynamico_init_metric(Ai,Av,fv,le_de,Riv2,wee) 317 print('...done')318 INFO('...done') 318 319 319 320 #------------------------ Mesh partitioning ------------------------
Note: See TracChangeset
for help on using the changeset viewer.