Return-Path: X-Original-To: apmail-systemml-dev-archive@minotaur.apache.org Delivered-To: apmail-systemml-dev-archive@minotaur.apache.org Received: from mail.apache.org (hermes.apache.org [140.211.11.3]) by minotaur.apache.org (Postfix) with SMTP id 92E9918F44 for ; Tue, 15 Dec 2015 08:03:14 +0000 (UTC) Received: (qmail 7303 invoked by uid 500); 15 Dec 2015 08:03:14 -0000 Delivered-To: apmail-systemml-dev-archive@systemml.apache.org Received: (qmail 7247 invoked by uid 500); 15 Dec 2015 08:03:14 -0000 Mailing-List: contact dev-help@systemml.incubator.apache.org; run by ezmlm Precedence: bulk List-Help: List-Unsubscribe: List-Post: List-Id: Reply-To: dev@systemml.incubator.apache.org Delivered-To: mailing list dev@systemml.incubator.apache.org Received: (qmail 7228 invoked by uid 99); 15 Dec 2015 08:03:14 -0000 Received: from Unknown (HELO spamd1-us-west.apache.org) (209.188.14.142) by apache.org (qpsmtpd/0.29) with ESMTP; Tue, 15 Dec 2015 08:03:14 +0000 Received: from localhost (localhost [127.0.0.1]) by spamd1-us-west.apache.org (ASF Mail Server at spamd1-us-west.apache.org) with ESMTP id C41E8C2FF9 for ; Tue, 15 Dec 2015 08:03:13 +0000 (UTC) X-Virus-Scanned: Debian amavisd-new at spamd1-us-west.apache.org X-Spam-Flag: NO X-Spam-Score: 5.667 X-Spam-Level: ***** X-Spam-Status: No, score=5.667 tagged_above=-999 required=6.31 tests=[HTML_IMAGE_ONLY_20=0.7, HTML_MESSAGE=3, MIME_HEADER_CTYPE_ONLY=1.996, MSGID_FROM_MTA_HEADER=0.001, RCVD_IN_MSPIKE_H3=-0.01, RCVD_IN_MSPIKE_WL=-0.01, SPF_PASS=-0.001, TVD_FW_GRAPHIC_NAME_MID=0.001, T_RP_MATCHES_RCVD=-0.01] autolearn=disabled Received: from mx1-us-east.apache.org ([10.40.0.8]) by localhost (spamd1-us-west.apache.org [10.40.0.7]) (amavisd-new, port 10024) with ESMTP id 0qsUIZ5CHAJ7 for ; Tue, 15 Dec 2015 08:03:02 +0000 (UTC) Received: from e34.co.us.ibm.com (e34.co.us.ibm.com [32.97.110.152]) by mx1-us-east.apache.org (ASF Mail Server at mx1-us-east.apache.org) with ESMTPS id 39E68429DD for ; Tue, 15 Dec 2015 08:03:00 +0000 (UTC) Received: from localhost by e34.co.us.ibm.com with IBM ESMTP SMTP Gateway: Authorized Use Only! Violators will be prosecuted for from ; Tue, 15 Dec 2015 01:03:00 -0700 Received: from d03dlp01.boulder.ibm.com (9.17.202.177) by e34.co.us.ibm.com (192.168.1.134) with IBM ESMTP SMTP Gateway: Authorized Use Only! Violators will be prosecuted; Tue, 15 Dec 2015 01:02:57 -0700 X-IBM-Helo: d03dlp01.boulder.ibm.com X-IBM-MailFrom: mboehm@us.ibm.com X-IBM-RcptTo: dev@systemml.incubator.apache.org Received: from b03cxnp08025.gho.boulder.ibm.com (b03cxnp08025.gho.boulder.ibm.com [9.17.130.17]) by d03dlp01.boulder.ibm.com (Postfix) with ESMTP id 0F0B41FF0043 for ; Tue, 15 Dec 2015 00:51:08 -0700 (MST) Received: from d03av03.boulder.ibm.com (d03av03.boulder.ibm.com [9.17.195.169]) by b03cxnp08025.gho.boulder.ibm.com (8.14.9/8.14.9/NCO v10.0) with ESMTP id tBF82v4l31129658 for ; Tue, 15 Dec 2015 01:02:57 -0700 Received: from d03av03.boulder.ibm.com (localhost [127.0.0.1]) by d03av03.boulder.ibm.com (8.14.4/8.14.4/NCO v10.0 AVout) with ESMTP id tBF82ung031911 for ; Tue, 15 Dec 2015 01:02:56 -0700 Received: from d50lp02.ny.us.ibm.com (d50lp02.pok.ibm.com [146.89.104.208]) by d03av03.boulder.ibm.com (8.14.4/8.14.4/NCO v10.0 AVin) with ESMTP id tBF82urF031883 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=OK) for ; Tue, 15 Dec 2015 01:02:56 -0700 Message-Id: <201512150802.tBF82urF031883@d03av03.boulder.ibm.com> Received: from /spool/local by d50lp02.ny.us.ibm.com with IBM ESMTP SMTP Gateway: Authorized Use Only! Violators will be prosecuted for from ; Tue, 15 Dec 2015 03:02:56 -0500 Received: from smtp.notes.na.collabserv.com (192.155.248.81) by d50lp02.ny.us.ibm.com (158.87.18.21) with IBM ESMTP SMTP Gateway: Authorized Use Only! Violators will be prosecuted; (version=TLSv1/SSLv3 cipher=AES256-SHA bits=256/256) Tue, 15 Dec 2015 03:02:55 -0500 Received: from /spool/local by smtp.notes.na.collabserv.com with smtp.notes.na.collabserv.com ESMTP for from ; Tue, 15 Dec 2015 08:02:54 -0000 Received: from us1a3-smtp04.a3.dal06.isc4sb.com (10.106.154.237) by smtp.notes.na.collabserv.com (10.106.227.88) with smtp.notes.na.collabserv.com ESMTP; Tue, 15 Dec 2015 08:02:51 -0000 Received: from us1a3-mail149.a3.dal06.isc4sb.com ([10.146.38.84]) by us1a3-smtp04.a3.dal06.isc4sb.com with ESMTP id 2015121508025009-55146 ; Tue, 15 Dec 2015 08:02:50 +0000 In-Reply-To: Subject: Re: Is there any equivalent of kron function in DML ? To: dev@systemml.incubator.apache.org From: "Matthias Boehm" Date: Tue, 15 Dec 2015 09:02:46 +0100 References: X-KeepSent: 9BB55274:D908F8D8-00257F1C:002A2E33; type=4; name=$KeepSent X-Mailer: IBM Notes Release 9.0.1FP2 SHF37 August 25, 2014 X-LLNOutbound: False X-Disclaimed: 56215 X-TNEFEvaluated: 1 Content-type: multipart/related; Boundary="0__=8FBBF58FDFB9A8A38f9e8a93df938690918c8FBBF58FDFB9A8A3" x-cbid: 15121508-0017-0000-0000-000010580549 X-IBM-ISS-SpamDetectors: Score=0.373977; BY=0.110255; FL=0; FP=0; FZ=0; HX=0; KW=0; PH=0; SC=0.373977; ST=0; TS=0; UL=0; ISC= X-IBM-ISS-DetailInfo: BY=3.00004670; HX=3.00000236; KW=3.00000007; PH=3.00000004; SC=3.00000125; SDB=6.00630859; UDB=6.00282367; UTC=2015-12-15 08:02:52 x-cbparentid: 15121508-4536-0000-0000-0000055A8594 X-TM-AS-MML: disable X-Content-Scanned: Fidelis XPS MAILER --0__=8FBBF58FDFB9A8A38f9e8a93df938690918c8FBBF58FDFB9A8A3 Content-type: multipart/alternative; Boundary="1__=8FBBF58FDFB9A8A38f9e8a93df938690918c8FBBF58FDFB9A8A3" --1__=8FBBF58FDFB9A8A38f9e8a93df938690918c8FBBF58FDFB9A8A3 Content-Transfer-Encoding: quoted-printable Content-type: text/plain; charset=US-ASCII Hi Sourav, well, we do not support Kronecker products K =3D kron(X,Y) yet. However, as= a workaround you could indeed create a custom permutation matrix P from X, where P[1:nrow(Y),] =3D diag(matrix(as.scalar(X[1,1]),nrow(Y),1)) P[nrow(Y)+1:2*nrow(Y),] =3D diag(matrix(as.scalar(X[1,2]),nrow(Y),1)) ... and subsequently create a linearized output with K =3D P %*% Y Finally, we provide a reshape matrix constructor matrix(K, rows=3D..., cols=3D...) that allows you to rearrange your entries as needed. If you run this on large data for X, however, we might want to put this on the roadmap for 2016. Could you share the characteristics of the workloads you have in mind? Regards, Matthias From: Sourav Mazumder To: dev@systemml.incubator.apache.org Date: 12/14/2015 11:56 PM Subject: Is there any equivalent of kron function in DML ? Hi, Is there any equivalent to kron function (of Matlab) in DML ? Can I use the matrix multiplication %*% for the same ? Regards, Sourav --1__=8FBBF58FDFB9A8A38f9e8a93df938690918c8FBBF58FDFB9A8A3 Content-Transfer-Encoding: quoted-printable Content-type: text/html; charset=US-ASCII Content-Disposition: inline

Hi Sourav,

well, we do not support Kronecker products= K =3D kron(X,Y) yet. However, as a workaround you could indeed create a cu= stom permutation matrix P from X, where

= P[1:nrow(Y),] =3D diag(matrix(as.scalar(X[1,1]),nrow(Y),1))
P[nrow(Y)+1:2*nrow(Y),] =3D diag(matrix(as.scalar(X[1= ,2]),nrow(Y),1))
...

and subsequently create a linearized = output with

K =3D P %*% Y

= Finally, we provide a reshape matrix constructor matrix(K, rows=3D..., cols= =3D...) that allows you to rearrange your entries as needed. If you run thi= s on large data for X, however, we might want to put this on the roadmap fo= r 2016. Could you share the characteristics of the workloads you have in mi= nd?


Regards,
Matthias

3D"InactiveSourav Mazumder ---12/14/2015 11:56:52 PM---Hi, Is there any e= quivalent to kron function (of Matlab) in DML ?

From: Sourav Mazumder = <sourav.mazumder00@gmail.com>
To: dev@systemml.incubator.apache.org<= /font>
Date: 12/14/2015 11:56 PM
Sub= ject: Is there any equivalent of kron functi= on in DML ?





Hi,

Is there any equivalen= t to kron function (of Matlab) in DML ?

Can I use the matrix multipl= ication %*% for the same ?

Regards,
Sourav


--1__=8FBBF58FDFB9A8A38f9e8a93df938690918c8FBBF58FDFB9A8A3-- --0__=8FBBF58FDFB9A8A38f9e8a93df938690918c8FBBF58FDFB9A8A3--