MAGICL:EXPM doesn't work on SINGLE-FLOAT or (COMPLEX SINGLE-FLOAT) matrices

[danikino]

After installing magicl from git : in emacs with sbcl and sly : (ql:quickload :magicl) (magicl.backends:active-backends) --> (:EXPOKIT :LAPACK :BLAS :LISP) (defparameter a (magicl:from-list '(1.0 0.0 0.0 0.0 2.0 0.0 0.0 0.0 3.0) '(3 3))) a -->

<MAGICL:MATRIX/DOUBLE-FLOAT (3x3):

1.000 0.000 0.000 0.000 2.000 0.000 0.000 0.000 3.000>

(magicl:logm a) -->

<MAGICL:MATRIX/DOUBLE-FLOAT (3x3):

0.000 0.000 0.000 0.000 0.693 0.000 0.000 0.000 1.099> now ,if i do : (magicl:expm a) --> MAGICL:EXPM is not implemented in any of the current active backends: EXPOKIT, LAPACK, BLAS, LISP. [Condition of type MAGICL.BACKENDS:NO-APPLICABLE-IMPLEMENTATION] the question is how to calculate the exponential of a matrix ? What am I doing wrong ? I am trying to use magicl to learn robotics :) :)

[stylewarning]

HI @danikino, sorry for the late response.

The issue is that expm only works on complex double-float matrices.

CL-USER> (defparameter *a (magicl:from-list '(1.0 0.0 0.0 0.0 2.0 0.0 0.0 0.0 3.0) '(3 3) :type '(complex double-float)))
*A
CL-USER> (magicl:expm *a)
#<MAGICL:MATRIX/COMPLEX-DOUBLE-FLOAT (3x3):
   2.718 + 0.000j     0.000 + 0.000j     0.000 + 0.000j
   0.000 + 0.000j     7.389 + 0.000j     0.000 + 0.000j
   0.000 + 0.000j     0.000 + 0.000j    20.086 + 0.000j>

We haven't extended it to get it to work on just double or single float matrices, but shouldn't be too hard!

[stylewarning]

I merged support for double-float, now you can do:

CL-USER> (setf *read-default-float-format* 'double-float)
DOUBLE-FLOAT
CL-USER> (defparameter *a (magicl:from-list '(1.0 0.0 0.0 0.0 2.0 0.0 0.0 0.0 3.0) '(3 3)))
*A
CL-USER> *a
#<MAGICL:MATRIX/DOUBLE-FLOAT (3x3):
   1.000     0.000     0.000
   0.000     2.000     0.000
   0.000     0.000     3.000>
CL-USER> (magicl:expm *a)
#<MAGICL:MATRIX/DOUBLE-FLOAT (3x3):
   2.718     0.000     0.000
   0.000     7.389     0.000
   0.000     0.000    20.086>

Still no implementation for single-floats though.