fortran - How to correctly use mkl_domatcopy from MKL? -


i need find faster way transpose matrix using mkl. , tried mkl_domatcopy mkl never right.

here test code(fortran):

        program main                  integer, parameter:: nrow = 3 !rows                 integer, parameter:: ncol = 3 !cols                  real*8, allocatable:: m(:,:)                 real*8, allocatable:: mt(:,:)                  integer:: i,j                  allocate(m(nrow,ncol))                 allocate(mt(nrow,ncol))                   = 1, nrow                 j = 1, ncol                         m(i,j)=i                 end                 end  call mkl_domatcopy("c","t",3,3,9,m,3,mt,3)                  print *,m                 print *,"************************"                 print *,mt          end

and output is:

   1.00000000000000        2.00000000000000        3.00000000000000    1.00000000000000        2.00000000000000        3.00000000000000    1.00000000000000        2.00000000000000        3.00000000000000  ************************            0

why mt 0?is beacause use wrong or something?

documents function: https://software.intel.com/en-us/node/520863

ps:i still don't 'alpha' means.

as user roygvib suggests in comments, including mkl.fi file give additional details.

this code

    program main              include 'mkl.fi'              integer, parameter:: nrow = 3 !rows             integer, parameter:: ncol = 3 !cols              real*8, allocatable:: m(:,:)             real*8, allocatable:: mt(:,:)             integer:: i,j              allocate(m(nrow,ncol))             allocate(mt(nrow,ncol))              = 1, nrow             j = 1, ncol                     m(i,j)=i             end             end              call mkl_domatcopy("c","t",3,3,9,m,3,mt,3)              print *,m             print *,"************************"             print *,mt      end

raises following error

test.f90(23): error #6633: type of actual argument differs type of dummy argument. [9] call mkl_domatcopy("c","t",3,3,9,m,3,mt,3) -----------------------------------------------^ compilation aborted test.f90 (code 1)

interestingly, if turn 9 double precision value (or variable) -- note here appended d0 suffix floating point value.

    program main              include 'mkl.fi'              integer, parameter:: nrow = 3 !rows             integer, parameter:: ncol = 3 !cols              real*8, allocatable:: m(:,:)             real*8, allocatable:: mt(:,:)             integer:: i,j              allocate(m(nrow,ncol))             allocate(mt(nrow,ncol))              = 1, nrow             j = 1, ncol                     m(i,j)=i             end             end              call mkl_domatcopy("c","t",3,3,9d0,m,3,mt,3)              print *,m             print *,"************************"             print *,mt      end

then application returns

$ ./test 1.00000000000000        2.00000000000000        3.00000000000000 1.00000000000000        2.00000000000000        3.00000000000000 1.00000000000000        2.00000000000000        3.00000000000000 ************************ 9.00000000000000        9.00000000000000        9.00000000000000 18.0000000000000        18.0000000000000        18.0000000000000 27.0000000000000        27.0000000000000        27.0000000000000

finally, respect alpha means, manual says

alpha parameter scales input matrix alpha.

and notice output tranposed , each element multiplied 9.


-

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