Skip to content
Draft
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 4 additions & 2 deletions transcendental/transcendental.lisp
Original file line number Diff line number Diff line change
Expand Up @@ -9,11 +9,13 @@
(let ((ideg 6)
(rows (magicl:nrows m))
(tcoef (coerce 1.0 'double-float))
(h (copy-seq (magicl::storage m)))
(h-tensor (magicl::deep-copy-tensor m))
(iexph 0)
(ns 0)
(iflag 0))
(let ((lwsp (+ (* 4 rows rows) ideg 1))
(when (eql :row-major (magicl::layout h-tensor)) (magicl:transpose! h-tensor))
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

transponse on new line

(let ((h (magicl::storage h-tensor))
(lwsp (+ (* 4 rows rows) ideg 1))
(ipiv (magicl::make-array (list rows) :element-type '(signed-byte 32))))
(let ((wsp (magicl::make-array (list lwsp) :element-type '(complex double-float))))
;; Requires direct foreign function call due to need to access a pointer
Expand Down