← Back to branch summary

~ubuntu-branches/ubuntu/natty/suitesparse/natty

« back to all changes in this revision

Viewing changes to CXSparse_newfiles/Demo/cs_idemo.c

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme
  • Date: 2006-12-22 10:16:15 UTC
  • Revision ID: james.westby@ubuntu.com-20061222101615-2ohaj8902oix2rnk
Tags: upstream-2.3.1
Import upstream version 2.3.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
CXSparse_newfiles/Demo/cs_idemo.c
 
1
#include "cs.h"
 
2
/* test real/complex conversion routines (int version) */
 
3
int main (void)
 
4
{
 
5
    cs_ci *T, *A, *A1, *A2, *B ;
 
6
    cs_di *C1, *C2, *Treal, *Timag ;
 
7
 
 
8
    printf ("\n--- cs_idemo\n") ;
 
9
 
 
10
    T = cs_ci_load (stdin) ;            /* load a complex triplet matrix, T */
 
11
    printf ("\nT:\n") ;
 
12
    cs_ci_print (T, 0) ;
 
13
 
 
14
    Treal = cs_i_real (T, 1) ;          /* Treal = real part of T */
 
15
    printf ("\nTreal:\n") ;
 
16
    cs_di_print (Treal, 0) ;
 
17
 
 
18
    Timag = cs_i_real (T, 0) ;          /* Treal = imaginary part of T */
 
19
    printf ("\nTimag:\n") ;
 
20
    cs_di_print (Timag, 0) ;
 
21
 
 
22
    A = cs_ci_compress (T) ;            /* A = compressed-column form of T */
 
23
    printf ("\nA:\n") ;
 
24
    cs_ci_print (A, 0) ;
 
25
 
 
26
    C1 = cs_i_real (A, 1) ;             /* C1 = real (A) */
 
27
    printf ("\nC1 = real(A):\n") ;
 
28
    cs_di_print (C1, 0) ;
 
29
 
 
30
    C2 = cs_i_real (A, 0) ;             /* C2 = imag (A) */
 
31
    printf ("\nC2 = imag(A):\n") ;
 
32
    cs_di_print (C2, 0) ;
 
33
 
 
34
    A1 = cs_i_complex (C1, 1) ;         /* A1 = complex version of C1 */
 
35
    printf ("\nA1:\n") ;
 
36
    cs_ci_print (A1, 0) ;
 
37
 
 
38
    A2 = cs_i_complex (C2, 0) ;         /* A2 = complex version of C2 (imag.) */
 
39
    printf ("\nA2:\n") ;
 
40
    cs_ci_print (A2, 0) ;
 
41
 
 
42
    B = cs_ci_add (A1, A2, 1., -1.) ;   /* B = A1 - A2 */
 
43
    printf ("\nB = conj(A):\n") ;
 
44
    cs_ci_print (B, 0) ;
 
45
 
 
46
    cs_ci_spfree (T) ;
 
47
    cs_ci_spfree (A) ;
 
48
    cs_ci_spfree (A1) ;
 
49
    cs_ci_spfree (A2) ;
 
50
    cs_ci_spfree (B) ;
 
51
    cs_di_spfree (C1) ;
 
52
    cs_di_spfree (C2) ;
 
53
    cs_di_spfree (Treal) ;
 
54
    cs_di_spfree (Timag) ;
 
55
 
 
56
    return (0) ;
 
57
}