← Back to branch summary

~ubuntu-branches/ubuntu/natty/deal.ii/natty

« back to all changes in this revision

Viewing changes to contrib/umfpack/UMFPACK/Source/umfpack_load_symbolic.c

  • Committer: Bazaar Package Importer
  • Author(s): Adam C. Powell, IV
  • Date: 2009-05-08 23:13:50 UTC
  • Revision ID: james.westby@ubuntu.com-20090508231350-rrh1ltgi0tifabwc
Tags: upstream-6.2.0
Import upstream version 6.2.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
contrib/umfpack/UMFPACK/Source/umfpack_load_symbolic.c
 
1
/* ========================================================================== */
 
2
/* === UMFPACK_load_symbolic ================================================ */
 
3
/* ========================================================================== */
 
4
 
 
5
/* -------------------------------------------------------------------------- */
 
6
/* UMFPACK Copyright (c) Timothy A. Davis, CISE,                              */
 
7
/* Univ. of Florida.  All Rights Reserved.  See ../Doc/License for License.   */
 
8
/* web: http://www.cise.ufl.edu/research/sparse/umfpack                       */
 
9
/* -------------------------------------------------------------------------- */
 
10
 
 
11
/*
 
12
    User-callable.  Loads a Symbolic object from a file created by
 
13
    umfpack_*_save_symbolic.
 
14
*/
 
15
 
 
16
#include "umf_internal.h"
 
17
#include "umf_valid_symbolic.h"
 
18
#include "umf_malloc.h"
 
19
#include "umf_free.h"
 
20
 
 
21
#define READ(object,type,n) \
 
22
{ \
 
23
    object = (type *) UMF_malloc (n, sizeof (type)) ; \
 
24
    if (object == (type *) NULL) \
 
25
    { \
 
26
        UMFPACK_free_symbolic ((void **) &Symbolic) ; \
 
27
        fclose (f) ; \
 
28
        return (UMFPACK_ERROR_out_of_memory) ; \
 
29
    } \
 
30
    if (fread (object, sizeof (type), n, f) != n) \
 
31
    { \
 
32
        UMFPACK_free_symbolic ((void **) &Symbolic) ; \
 
33
        fclose (f) ; \
 
34
        return (UMFPACK_ERROR_file_IO) ; \
 
35
    } \
 
36
    if (ferror (f)) \
 
37
    { \
 
38
        UMFPACK_free_symbolic ((void **) &Symbolic) ; \
 
39
        fclose (f) ; \
 
40
        return (UMFPACK_ERROR_file_IO) ; \
 
41
    } \
 
42
}
 
43
 
 
44
/* ========================================================================== */
 
45
/* === UMFPACK_load_symbolic ================================================ */
 
46
/* ========================================================================== */
 
47
 
 
48
GLOBAL Int UMFPACK_load_symbolic
 
49
(
 
50
    void **SymbolicHandle,
 
51
    char *user_filename
 
52
)
 
53
{
 
54
    SymbolicType *Symbolic ;
 
55
    char *filename ;
 
56
    FILE *f ;
 
57
 
 
58
    *SymbolicHandle = (void *) NULL ;
 
59
 
 
60
    /* ---------------------------------------------------------------------- */
 
61
    /* get the filename, or use the default name if filename is NULL */
 
62
    /* ---------------------------------------------------------------------- */
 
63
 
 
64
    if (user_filename == (char *) NULL)
 
65
    {
 
66
        filename = "symbolic.umf" ;
 
67
    }
 
68
    else
 
69
    {
 
70
        filename = user_filename ;
 
71
    }
 
72
    f = fopen (filename, "rb") ;
 
73
    if (!f)
 
74
    {
 
75
        return (UMFPACK_ERROR_file_IO) ;
 
76
    }
 
77
 
 
78
    /* ---------------------------------------------------------------------- */
 
79
    /* read the Symbolic header from the file, in binary */
 
80
    /* ---------------------------------------------------------------------- */
 
81
 
 
82
    Symbolic = (SymbolicType *) UMF_malloc (1, sizeof (SymbolicType)) ;
 
83
    if (Symbolic == (SymbolicType *) NULL)
 
84
    {
 
85
        fclose (f) ;
 
86
        return (UMFPACK_ERROR_out_of_memory) ;
 
87
    }
 
88
    if (fread (Symbolic, sizeof (SymbolicType), 1, f) != 1)
 
89
    {
 
90
        (void) UMF_free ((void *) Symbolic) ;
 
91
        fclose (f) ;
 
92
        return (UMFPACK_ERROR_file_IO) ;
 
93
    }
 
94
    if (ferror (f))
 
95
    {
 
96
        (void) UMF_free ((void *) Symbolic) ;
 
97
        fclose (f) ;
 
98
        return (UMFPACK_ERROR_file_IO) ;
 
99
    }
 
100
 
 
101
    if (Symbolic->valid != SYMBOLIC_VALID || Symbolic->n_row <= 0 ||
 
102
        Symbolic->n_col <= 0 || Symbolic->nfr < 0 || Symbolic->nchains < 0 ||
 
103
        Symbolic->esize < 0)
 
104
    {
 
105
        /* Symbolic does not point to a Symbolic object */
 
106
        (void) UMF_free ((void *) Symbolic) ;
 
107
        fclose (f) ;
 
108
        return (UMFPACK_ERROR_invalid_Symbolic_object) ;
 
109
    }
 
110
 
 
111
    Symbolic->Cperm_init         = (Int *) NULL ;
 
112
    Symbolic->Rperm_init         = (Int *) NULL ;
 
113
    Symbolic->Front_npivcol      = (Int *) NULL ;
 
114
    Symbolic->Front_parent       = (Int *) NULL ;
 
115
    Symbolic->Front_1strow       = (Int *) NULL ;
 
116
    Symbolic->Front_leftmostdesc = (Int *) NULL ;
 
117
    Symbolic->Chain_start        = (Int *) NULL ;
 
118
    Symbolic->Chain_maxrows      = (Int *) NULL ;
 
119
    Symbolic->Chain_maxcols      = (Int *) NULL ;
 
120
    Symbolic->Cdeg               = (Int *) NULL ;
 
121
    Symbolic->Rdeg               = (Int *) NULL ;
 
122
    Symbolic->Esize              = (Int *) NULL ;
 
123
    Symbolic->Diagonal_map       = (Int *) NULL ;
 
124
 
 
125
    /* umfpack_free_symbolic can now be safely called if an error occurs */
 
126
 
 
127
    /* ---------------------------------------------------------------------- */
 
128
    /* read the rest of the Symbolic object */
 
129
    /* ---------------------------------------------------------------------- */
 
130
 
 
131
    READ (Symbolic->Cperm_init,         Int, Symbolic->n_col+1) ;
 
132
    READ (Symbolic->Rperm_init,         Int, Symbolic->n_row+1) ;
 
133
    READ (Symbolic->Front_npivcol,      Int, Symbolic->nfr+1) ;
 
134
    READ (Symbolic->Front_parent,       Int, Symbolic->nfr+1) ;
 
135
    READ (Symbolic->Front_1strow,       Int, Symbolic->nfr+1) ;
 
136
    READ (Symbolic->Front_leftmostdesc, Int, Symbolic->nfr+1) ;
 
137
    READ (Symbolic->Chain_start,        Int, Symbolic->nchains+1) ;
 
138
    READ (Symbolic->Chain_maxrows,      Int, Symbolic->nchains+1) ;
 
139
    READ (Symbolic->Chain_maxcols,      Int, Symbolic->nchains+1) ;
 
140
    READ (Symbolic->Cdeg,               Int, Symbolic->n_col+1) ;
 
141
    READ (Symbolic->Rdeg,               Int, Symbolic->n_row+1) ;
 
142
    if (Symbolic->esize > 0)
 
143
    {
 
144
        /* only when dense rows are present */
 
145
        READ (Symbolic->Esize, Int, Symbolic->esize) ;
 
146
    }
 
147
    if (Symbolic->prefer_diagonal)
 
148
    {
 
149
        /* only when diagonal pivoting is prefered */
 
150
        READ (Symbolic->Diagonal_map, Int, Symbolic->n_col+1) ;
 
151
    }
 
152
 
 
153
    /* close the file */
 
154
    fclose (f) ;
 
155
 
 
156
    /* make sure the Symbolic object is valid */
 
157
    if (!UMF_valid_symbolic (Symbolic))
 
158
    {
 
159
        UMFPACK_free_symbolic ((void **) &Symbolic) ;
 
160
        return (UMFPACK_ERROR_invalid_Symbolic_object) ;
 
161
    }
 
162
 
 
163
    *SymbolicHandle = (void *) Symbolic ;
 
164
    return (UMFPACK_OK) ;
 
165
}