From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 28985 invoked by alias); 3 Dec 2001 04:07:34 -0000 Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm Precedence: bulk List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sources.redhat.com Received: (qmail 28972 invoked from network); 3 Dec 2001 04:07:32 -0000 Received: from unknown (HELO sys713.kaist.ac.kr) (143.248.51.23) by sources.redhat.com with SMTP; 3 Dec 2001 04:07:32 -0000 Received: (qmail 37163 invoked from network); 3 Dec 2001 04:03:45 -0000 Received: from casad13.kaist.ac.kr (HELO sys713.kaist.ac.kr) (143.248.59.74) by 0 with SMTP; 3 Dec 2001 04:03:45 -0000 Message-ID: <3C0AFA50.4060604@sys713.kaist.ac.kr> Date: Tue, 27 Nov 2001 12:54:00 -0000 From: "Dan, Ho-Jin" User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:0.9.6) Gecko/20011120 X-Accept-Language: ko, en-us MIME-Version: 1.0 To: gsl-discuss Subject: Re: matrix multiplication References: <1007346875.25380.0.camel@mrburns> Content-Type: multipart/mixed; boundary="------------030801030402030409030309" X-SW-Source: 2001-q4/txt/msg00081.txt.bz2 Message-ID: <20011127125400.63VtTLaV1iem6JSid5fHZICJJwZFDwd6840ka-wqUlU@z> This is a multi-part message in MIME format. --------------030801030402030409030309 Content-Type: text/plain; charset=EUC-KR Content-Transfer-Encoding: 7bit Content-length: 483 try to use level 3 of blas. gsl supports blas interfaces. A simple example is attatched. gsl_gemm.c LDFLAGS=-lgsl -lgslcblas -lm yours Dan, Ho-Jin Matthew J. Doller wrote: >someone correct me if i am wrong, but is there really no way of >multiplying two matrices, or a matrix and a vector with gsl >i know i can go element by element, but i'm looking for something of the >form: > >return_matrix = gsl_matrix_mul_matrix_by_matrix ( matrix_one , >matrix_two ); > >thanks! >matt > --------------030801030402030409030309 Content-Type: text/plain; name="gsl_gemm.c" Content-Transfer-Encoding: base64 Content-Disposition: inline; filename="gsl_gemm.c" Content-length: 887 I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxnc2wvZ3NsX21hdHJpeC5o PgojaW5jbHVkZSA8Z3NsL2dzbF9ibGFzLmg+CgppbnQgbWFpbigpCnsKICBk b3VibGUgQVtdID0geyAKICAgIDEuLCAyLiwgCiAgICAzLiwgNC4gCiAgfTsK ICBkb3VibGUgQltdID0gewogICAgMi4sIDMuLAogICAgNC4sIDUuCiAgfTsK ICBkb3VibGUgYWxwaGEgPSAxLiwgYmV0YSA9IDAuOwoKICBnc2xfbWF0cml4 X3ZpZXcgQV9tID0gZ3NsX21hdHJpeF92aWV3X2FycmF5KEEsIDIsIDIpOwog IGdzbF9tYXRyaXhfdmlldyBCX20gPSBnc2xfbWF0cml4X3ZpZXdfYXJyYXko QiwgMiwgMik7CiAgZ3NsX21hdHJpeCAqQyA9IGdzbF9tYXRyaXhfYWxsb2Mo MiwyKTsKCiAgZ3NsX2JsYXNfZGdlbW0oQ2JsYXNOb1RyYW5zLCBDYmxhc05v VHJhbnMsIGFscGhhLCAmQV9tLm1hdHJpeCwgCiAgICAgICZCX20ubWF0cml4 LCBiZXRhLCBDKTsKCiAgcHJpbnRmKCJbW1xuJWZcdCVmXG4lZlx0JWZdXVxu IiwgCiAgICAgIGdzbF9tYXRyaXhfZ2V0KEMsIDAsIDApLCBnc2xfbWF0cml4 X2dldChDLCAwLCAxKSwgCiAgICAgIGdzbF9tYXRyaXhfZ2V0KEMsIDEsIDAp LCBnc2xfbWF0cml4X2dldChDLCAxLCAxKSk7CgogIGdzbF9tYXRyaXhfZnJl ZShDKTsKICByZXR1cm4gMDsKfQoKCg== --------------030801030402030409030309--