From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 24631 invoked by alias); 27 Nov 2003 09:14:04 -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 24623 invoked from network); 27 Nov 2003 09:14:03 -0000 Received: from unknown (HELO smtp0.libero.it) (193.70.192.33) by sources.redhat.com with SMTP; 27 Nov 2003 09:14:03 -0000 Received: from libero.it (193.70.192.57) by smtp0.libero.it (7.0.020-DD01) id 3F6F1CE30081F705 for gsl-discuss@sources.redhat.com; Thu, 27 Nov 2003 10:14:03 +0100 Date: Thu, 27 Nov 2003 09:14:00 -0000 Message-Id: Subject: linear system MIME-Version: 1.0 X-Sensitivity: 3 Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: quoted-printable From: "lauradiara\@libero\.it" To: "gsl-discuss" X-XaM3-API-Version: 4.1 (B19) X-type: 0 X-SenderIP: 193.76.233.91 X-SW-Source: 2003-q4/txt/msg00097.txt.bz2 Hi all, I've implemented a function which solves linear system A x =3D b with the u= ser input. Reading the docs I have noticed that in gsl_linalg.h there are two very sim= ilar methods which differ for the type of input parameters: gsl_linalg_LU_solve = and gsl_linalg_complex_LU_solve.=0D What are the differences? (I know that the second one works on complex numb= ers.) Could I use only gsl_linalg_complex_LU_solve? If the user insert double num= bers, I insert the img part of complex number=3D0. Is it true? What are the real usages of gsl_linalg_complex_LU_solve? Other question: how can I solve linear system A x =3D b with indinite solut= ions? For example: we suppose to have: x+2y+3z=3D3 2x+3y+8z=3D4 3x+2y+17z=3D1 Thi system have infinite solutions: x=3D-1-7t y=3D2+2t=0D z=3Dt t in R Can i use gsl_linalg_LU_solve or gsl_linalg_complex_LU_solve to solve this = kind of linear system? Cheers and thanks for your support --Laura