lxsd / fvn

Blame view

fvn_sparse/UMFPACK/MATLAB/umfpack_simple.m.out 1.03 KB
edit raw normal view history
422234dc3   daniau   git-svn-id: https... 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
 
  umfpack_simple
   UMFPACK_SIMPLE a simple demo
   
    Example:
      umfpack_simple
   
    Copyright 1995-2007 by Timothy A. Davis.
   
    UMFPACK License:
   
        Your use or distribution of UMFPACK or any modified version of
        UMFPACK implies that you agree to this License.  UMFPACK is
        is free software; you can redistribute it and/or
        modify it under the terms of the GNU Lesser General Public
        License as published by the Free Software Foundation; either
        version 2.1 of the License, or (at your option) any later version.
  
  A = 
       2     3     0     0     0
       3     0     4     0     6
       0    -1    -3     2     0
       0     0     1     0     0
       0     4     2     0     1
  b = 
       8
      45
      -3
       3
      19
  Solution to Ax=b via UMFPACK:
  x1 = umfpack2 (A, '\', b)
  x1 = 
      1.0000
      2.0000
      3.0000
      4.0000
      5.0000
  Solution to Ax=b via MATLAB:
  x2 = A\b
  x2 = 
      1.0000
      2.0000
      3.0000
      4.0000
      5.0000
  norm (x1-x2) should be small: 1.28037e-15
  Type 'umfpack_demo' for a full demo of UMFPACK
  diary off