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 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75
|
>> umfpack_simple
umfpack_simple: a simple demo of UMFPACK
UMFPACK Version 4.4, Copyright (c) 2005 by Timothy A. Davis.
All Rights Reserved. Type umfpack_details for License.
UMFPACK License:
Your use or distribution of UMFPACK or any modified version of
UMFPACK implies that you agree to this License.
THIS MATERIAL IS PROVIDED AS IS, WITH ABSOLUTELY NO WARRANTY
EXPRESSED OR IMPLIED. ANY USE IS AT YOUR OWN RISK.
Permission is hereby granted to use or copy this program, provided
that the Copyright, this License, and the Availability of the original
version is retained on all copies. User documentation of any code that
uses UMFPACK or any modified version of UMFPACK code must cite the
Copyright, this License, the Availability note, and "Used by permission."
Permission to modify the code and to distribute modified code is granted,
provided the Copyright, this License, and the Availability note are
retained, and a notice that the code was modified is included. This
software was developed with support from the National Science Foundation,
and is provided to you free of charge.
Availability: http://www.cise.ufl.edu/research/sparse/umfpack
See also: umfpack, umfpack_details
Hit enter to agree to the above License:
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 = umfpack (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
|