scipy - Calling MKL from Python : DSTEVR -


dstevr computes eigensolutions of triangular symmetric matrix. cool. except not 1 of routines ported wrapper scipy. i've followed instructions on how call mkl directly python , attached seems give correct answer. gosh.... there someway clean up?!

  import numpy np    scipy import linalg   ctypes import *     c_double_p = pointer(c_double)   c_int_p = pointer(c_int)   c_char_p = pointer(c_char)    mkl = cdll('mkl_rt.dll')   dstevr = mkl.dstevr    #subroutine dstevr( jobz, range, n, d, e, vl, vu, il, iu, abstol, m,    # *      w, z, ldz, isuppz, work, lwork, iwork, liwork, info)   #  character * 1 jobz, range   #  integer n, il, iu, m, ldz, lwork, liwork, info   #  integer isuppz(*), iwork(*)   #  double precision vl, vu, abstol   #  double precision d(*), e(*), w(*), z(ldz,*), work(*)    dstevr.argtypes = [ c_char_p, c_char_p, c_int_p, c_double_p, c_double_p,  c_double_p, c_double_p, c_int_p, c_int_p, c_double_p,  \  c_int_p, c_double_p, c_double_p, c_int_p, c_int_p, c_double_p, c_int_p, c_int_p, c_int_p, c_int_p]     sv = "v"   eig_j = c_char_p(c_char(sv[0]))   sr = "a"   eig_r = c_char_p(c_char(sr[0]))    vl = c_double(0.0)   vu = c_double(0.0)   il = c_int(0)   iu = c_int(0)   abstol = c_double(0.0)   cm = c_int(0)   eig_info = c_int(0)    n = 6   cn = c_int(n)   ldz = cn     lwork = c_int(20*n)   liwork = c_int(10*n)    diag = np.ascontiguousarray(np.ones(n)*2)   diag_p = diag.ctypes.data_as(c_double_p)    offdiag = np.ascontiguousarray(np.ones(n-1)*(-1))   offdiag_p = offdiag.ctypes.data_as(c_double_p)    isuppz = np.ascontiguousarray(np.ones(n*2),dtype=int)   isuppz_p = isuppz.ctypes.data_as(c_int_p)    eigw = np.ascontiguousarray(np.zeros(n))   eigw_p = eigw.ctypes.data_as(c_double_p)    workz = np.ascontiguousarray(np.ones(20*n))   workz_p = workz.ctypes.data_as(c_double_p)    iworkz = np.ascontiguousarray(np.ones(10*n),dtype=int)   iworkz_p = iworkz.ctypes.data_as(c_int_p)    eigz = np.ascontiguousarray(np.ones(n*n))   eigz_p = eigz.ctypes.data_as(c_double_p)     dstevr( eig_j, eig_r, byref(cn), diag_p, offdiag_p, byref(vl),byref(vu),byref(il),byref(iu),byref(abstol),byref(cm),\   eigw_p, eigz_p, byref(ldz), isuppz_p, workz_p, byref(lwork), iworkz_p, byref(liwork), byref(eig_info))     print "eig_info", eig_info    print eigz    = np.eye(n,n,k=-1)*(-1) + np.eye(n,n)*2 + np.eye(n,n,k=1)*(-1)    w,v= linalg.eigh(a)    print v.t

thanks,

you can use instead cython wrapper in cython_lapack, http://docs.scipy.org/doc/scipy/reference/linalg.cython_lapack.html

you'll need cython though.


-

Comments

Post a Comment

Popular posts from this blog

amazon web services - S3 Pre-signed POST validate file type? -

i know it's possible put limit on file size content-length-range header. possible validate file type? https://docs.aws.amazon.com/amazons3/latest/dev/httppostforms.html#policyconditions i see there content-type header, if set say, audio/mp3 allow mp3 files , return error if file not mp3? i found previous question answers mention validating file size: s3 direct upload restricting file size , type you can specify content-type @ post request. you can specify @ signature, post must done content-type: { "expiration": "2007-12-01t12:00:00.000z", "conditions": [ {"acl": "public-read" }, {"bucket": "johnsmith" }, {"content-type: "audio/mp3"} ] } creating html form (using aws signature version 4) edit: @ previous question found, checking content-type.
Read more

c# - Check Keyboard Input Winforms -

i wondering if instead of doing this protected override void onkeydown(keyeventargs e) { if (e.keycode == keys.a) console.writeline("the key down."); } i set bool method , this: if(keydown(keys.a)) // whatever here i've been sitting here ages trying figure out how it. can't wrap head around it. in case wondering, plan call bool inside different method, check input. since want perform action after pressing key, using keydown event enough. but in cases suppose want check if specific key down in middle of process, can use getkeystate method way: [dllimport("user32.dll", charset = charset.auto, exactspelling = true)] public static extern short getkeystate(int keycode); public const int key_pressed = 0x8000; public static bool iskeydown(keys key) { return convert.toboolean(getkeystate((int)key) & key_pressed); } you should know, each time check key state using example iskeydown(keys.a) method returns true if ke...
Read more