RE: python libraries
I should have said: MicroPython is not CPython. Cython is a tool for CPython, i.e. translating existing Python code to C for speed reasons and make the bindings. The reason behind is that Python interpreters are slow, for example compared to Javascript interpreters, and are really slow compared to native code. Cython is a way for Python programmers who do not want to learn a compiled language to improve the speed of their code. I'm not convinced it's a good idea, some developers like Python a lot because they have a lot of tools easily available, but the drawback is that complexity is hidden, and auto-translating code that has bad complexity to C will not make the complexity any better. And even if your code has the right complexity, it has probably a style that once translated will not be full optimized.
But anyway that's not what would have to be done here: numpy is a collection of CPython bindings to existing C code (blas+lapack). Porting it to MicroPython means a complete rewrite of the bindings.
Below is a partial listing of the code I'm using to emulate some Numpy commands. As you can see it is really MicroPython code, I don't think you can automatically get this from CPython code.
Code:
static mp_obj_t linalg_zerosones(size_t n_args, const mp_obj_t *args,int v0,int v1) {
size_t n1=0,m1=0;
if (MP_OBJ_IS_SMALL_INT(args[0]))
n1= MP_OBJ_SMALL_INT_VALUE(args[0]);
if (n_args==2 && MP_OBJ_IS_SMALL_INT(args[1]))
m1= MP_OBJ_SMALL_INT_VALUE(args[1]);
if (v0==0 && v1==1)
m1=n1;
mp_obj_t r = mp_obj_new_list(0, NULL);
bool alea=v0==-1;
for (int i=0;i<n1;++i){
if (m1){
mp_obj_t l = mp_obj_new_list(0, NULL);
for (int j=0;j<m1;++j){
mp_obj_list_append(l, mp_obj_new_float(alea?rand()/(RAND_MAX+1.0):(i==j?v1:v0)));
}
mp_obj_list_append(r, l);
}
else
mp_obj_list_append(r, mp_obj_new_float(alea?rand()/(RAND_MAX+1.0):v0));
}
return r;
}
static mp_obj_t linalg_zeros(size_t n_args, const mp_obj_t *args) {
return linalg_zerosones(n_args,args,0,0);
}
MP_DEFINE_CONST_FUN_OBJ_VAR_BETWEEN(linalg_zeros_obj, 1, 2, linalg_zeros);
static mp_obj_t linalg_ones(size_t n_args, const mp_obj_t *args) {
return linalg_zerosones(n_args,args,1,1);
}
MP_DEFINE_CONST_FUN_OBJ_VAR_BETWEEN(linalg_ones_obj, 1, 2, linalg_ones);
static mp_obj_t linalg_eye(size_t n_args, const mp_obj_t *args) {
return linalg_zerosones(n_args,args,0,1);
}
...
static const mp_rom_map_elem_t linalg_locals_dict_table[] = {
{ MP_ROM_QSTR(MP_QSTR_pi), mp_const_float_pi },
{ MP_ROM_QSTR(MP_QSTR_matrix), (mp_obj_t) &linalg_matrix_obj },
{ MP_ROM_QSTR(MP_QSTR_arange), (mp_obj_t) &linalg_arange_obj },
{ MP_ROM_QSTR(MP_QSTR_linspace), (mp_obj_t) &linalg_linspace_obj },
{ MP_ROM_QSTR(MP_QSTR_dot), (mp_obj_t) &linalg_dot_obj },
{ MP_ROM_QSTR(MP_QSTR_cross), (mp_obj_t) &linalg_cross_obj },
{ MP_ROM_QSTR(MP_QSTR_solve), (mp_obj_t) &linalg_solve_obj },
{ MP_ROM_QSTR(MP_QSTR_peval), (mp_obj_t) &linalg_peval_obj },
{ MP_ROM_QSTR(MP_QSTR_horner), (mp_obj_t) &linalg_peval_obj },
{ MP_ROM_QSTR(MP_QSTR_mul), (mp_obj_t) &linalg_dot_obj },
{ MP_ROM_QSTR(MP_QSTR_transpose), (mp_obj_t) &linalg_transpose_obj },
{ MP_ROM_QSTR(MP_QSTR_rref), (mp_obj_t) &linalg_rref_obj },
{ MP_ROM_QSTR(MP_QSTR_egv), (mp_obj_t) &linalg_egv_obj },
{ MP_ROM_QSTR(MP_QSTR_eigenvects), (mp_obj_t) &linalg_egv_obj },
{ MP_ROM_QSTR(MP_QSTR_eig), (mp_obj_t) &linalg_eig_obj },
{ MP_ROM_QSTR(MP_QSTR_det), (mp_obj_t) &linalg_det_obj },
{ MP_ROM_QSTR(MP_QSTR_zeros), (mp_obj_t) &linalg_zeros_obj },
{ MP_ROM_QSTR(MP_QSTR_ones), (mp_obj_t) &linalg_ones_obj },
{ MP_ROM_QSTR(MP_QSTR_eye), (mp_obj_t) &linalg_eye_obj },
{ MP_ROM_QSTR(MP_QSTR_identity), (mp_obj_t) &linalg_eye_obj },
{ MP_ROM_QSTR(MP_QSTR_idn), (mp_obj_t) &linalg_eye_obj },
{ MP_ROM_QSTR(MP_QSTR_rand), (mp_obj_t) &linalg_rand_obj },
{ MP_ROM_QSTR(MP_QSTR_ranv), (mp_obj_t) &linalg_rand_obj },
{ MP_ROM_QSTR(MP_QSTR_ranm), (mp_obj_t) &linalg_rand_obj },
{ MP_ROM_QSTR(MP_QSTR_shape), (mp_obj_t) &linalg_shape_obj },
{ MP_ROM_QSTR(MP_QSTR_size), (mp_obj_t) &linalg_size_obj },
{ MP_ROM_QSTR(MP_QSTR_inv), (mp_obj_t) &linalg_inv_obj },
{ MP_ROM_QSTR(MP_QSTR_re), (mp_obj_t) &linalg_re_obj },
{ MP_ROM_QSTR(MP_QSTR_proot), (mp_obj_t) &linalg_proot_obj },
{ MP_ROM_QSTR(MP_QSTR_pcoeff), (mp_obj_t) &linalg_pcoeff_obj },
{ MP_ROM_QSTR(MP_QSTR_fft), (mp_obj_t) &linalg_fft_obj },
{ MP_ROM_QSTR(MP_QSTR_ifft), (mp_obj_t) &linalg_ifft_obj },
{ MP_ROM_QSTR(MP_QSTR_apply), (mp_obj_t) &linalg_apply_obj },
{ MP_ROM_QSTR(MP_QSTR_add), (mp_obj_t) &linalg_add_obj },
{ MP_ROM_QSTR(MP_QSTR_sub), (mp_obj_t) &linalg_sub_obj },
{ MP_ROM_QSTR(MP_QSTR_re), (mp_obj_t) &linalg_re_obj },
{ MP_ROM_QSTR(MP_QSTR_real), (mp_obj_t) &linalg_re_obj },
{ MP_ROM_QSTR(MP_QSTR_im), (mp_obj_t) &linalg_im_obj },
{ MP_ROM_QSTR(MP_QSTR_imag), (mp_obj_t) &linalg_im_obj },
{ MP_ROM_QSTR(MP_QSTR_conj), (mp_obj_t) &linalg_conj_obj },
{ MP_ROM_QSTR(MP_QSTR_abs), (mp_obj_t) &linalg_abs_obj },
};
static MP_DEFINE_CONST_DICT(linalg_locals_dict, linalg_locals_dict_table);
const mp_obj_type_t linalg_type = {
{ &mp_type_type },
.name = MP_QSTR_linalg,
.locals_dict = (mp_obj_t)&linalg_locals_dict
};
STATIC const mp_rom_map_elem_t mp_module_linalg_globals_table[] = {
{ MP_ROM_QSTR(MP_QSTR___name__), MP_ROM_QSTR(MP_QSTR__linalg) },
{ MP_ROM_QSTR(MP_QSTR_pi), mp_const_float_pi },
{ MP_ROM_QSTR(MP_QSTR_matrix), (mp_obj_t) &linalg_matrix_obj },
{ MP_ROM_QSTR(MP_QSTR_arange), (mp_obj_t) &linalg_arange_obj },
{ MP_ROM_QSTR(MP_QSTR_linspace), (mp_obj_t) &linalg_linspace_obj },
{ MP_ROM_QSTR(MP_QSTR_dot), (mp_obj_t) &linalg_dot_obj },
{ MP_ROM_QSTR(MP_QSTR_cross), (mp_obj_t) &linalg_cross_obj },
{ MP_ROM_QSTR(MP_QSTR_solve), (mp_obj_t) &linalg_solve_obj },
{ MP_ROM_QSTR(MP_QSTR_peval), (mp_obj_t) &linalg_peval_obj },
{ MP_ROM_QSTR(MP_QSTR_horner), (mp_obj_t) &linalg_peval_obj },
{ MP_ROM_QSTR(MP_QSTR_mul), (mp_obj_t) &linalg_dot_obj },
{ MP_ROM_QSTR(MP_QSTR_transpose), (mp_obj_t) &linalg_transpose_obj },
{ MP_ROM_QSTR(MP_QSTR_rref), (mp_obj_t) &linalg_rref_obj },
{ MP_ROM_QSTR(MP_QSTR_egv), (mp_obj_t) &linalg_egv_obj },
{ MP_ROM_QSTR(MP_QSTR_eigenvects), (mp_obj_t) &linalg_egv_obj },
{ MP_ROM_QSTR(MP_QSTR_eig), (mp_obj_t) &linalg_eig_obj },
{ MP_ROM_QSTR(MP_QSTR_det), (mp_obj_t) &linalg_det_obj },
{ MP_ROM_QSTR(MP_QSTR_zeros), (mp_obj_t) &linalg_zeros_obj },
{ MP_ROM_QSTR(MP_QSTR_ones), (mp_obj_t) &linalg_ones_obj },
{ MP_ROM_QSTR(MP_QSTR_eye), (mp_obj_t) &linalg_eye_obj },
{ MP_ROM_QSTR(MP_QSTR_identity), (mp_obj_t) &linalg_eye_obj },
{ MP_ROM_QSTR(MP_QSTR_idn), (mp_obj_t) &linalg_eye_obj },
{ MP_ROM_QSTR(MP_QSTR_rand), (mp_obj_t) &linalg_rand_obj },
{ MP_ROM_QSTR(MP_QSTR_ranv), (mp_obj_t) &linalg_rand_obj },
{ MP_ROM_QSTR(MP_QSTR_ranm), (mp_obj_t) &linalg_rand_obj },
{ MP_ROM_QSTR(MP_QSTR_shape), (mp_obj_t) &linalg_shape_obj },
{ MP_ROM_QSTR(MP_QSTR_size), (mp_obj_t) &linalg_size_obj },
{ MP_ROM_QSTR(MP_QSTR_inv), (mp_obj_t) &linalg_inv_obj },
{ MP_ROM_QSTR(MP_QSTR_re), (mp_obj_t) &linalg_re_obj },
{ MP_ROM_QSTR(MP_QSTR_proot), (mp_obj_t) &linalg_proot_obj },
{ MP_ROM_QSTR(MP_QSTR_pcoeff), (mp_obj_t) &linalg_pcoeff_obj },
{ MP_ROM_QSTR(MP_QSTR_fft), (mp_obj_t) &linalg_fft_obj },
{ MP_ROM_QSTR(MP_QSTR_ifft), (mp_obj_t) &linalg_ifft_obj },
{ MP_ROM_QSTR(MP_QSTR_apply), (mp_obj_t) &linalg_apply_obj },
{ MP_ROM_QSTR(MP_QSTR_add), (mp_obj_t) &linalg_add_obj },
{ MP_ROM_QSTR(MP_QSTR_sub), (mp_obj_t) &linalg_sub_obj },
{ MP_ROM_QSTR(MP_QSTR_re), (mp_obj_t) &linalg_re_obj },
{ MP_ROM_QSTR(MP_QSTR_real), (mp_obj_t) &linalg_re_obj },
{ MP_ROM_QSTR(MP_QSTR_im), (mp_obj_t) &linalg_im_obj },
{ MP_ROM_QSTR(MP_QSTR_imag), (mp_obj_t) &linalg_im_obj },
{ MP_ROM_QSTR(MP_QSTR_conj), (mp_obj_t) &linalg_conj_obj },
{ MP_ROM_QSTR(MP_QSTR_abs), (mp_obj_t) &linalg_abs_obj },
};
STATIC const mp_obj_dict_t mp_module_linalg_globals = {
.base = {&mp_type_dict},
.map = {
.all_keys_are_qstrs = 1,
.is_fixed = 1,
.used = MP_ARRAY_SIZE(mp_module_linalg_globals_table),
.alloc = MP_ARRAY_SIZE(mp_module_linalg_globals_table),
.table = (mp_map_elem_t*)mp_module_linalg_globals_table,
},
};
const mp_obj_module_t mp_module_linalg = {
.base = { &mp_type_module },
.globals = (mp_obj_dict_t*)&mp_module_linalg_globals,
};
|
Search
Member List
Calendar
Help
