1 /**********
2 *
3 * Copyright (c) 2003, Division of Imaging Science and Biomedical Engineering,
4 * University of Manchester, UK. All rights reserved.
5 *
6 * Redistribution and use in source and binary forms, with or without modification,
7 * are permitted provided that the following conditions are met:
8 *
9 * . Redistributions of source code must retain the above copyright notice,
10 * this list of conditions and the following disclaimer.
11 *
12 * . Redistributions in binary form must reproduce the above copyright notice,
13 * this list of conditions and the following disclaimer in the documentation
14 * and/or other materials provided with the distribution.
15 *
16 * . Neither the name of the University of Manchester nor the names of its
17 * contributors may be used to endorse or promote products derived from this
18 * software without specific prior written permission.
19 *
20 *
21 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
22 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
23 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
24 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
25 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
26 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
27 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
28 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
29 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
30 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
31 * POSSIBILITY OF SUCH DAMAGE.
32 *
33 **********
34 *
35 * Program : TINA
36 * File : $Source: /home/tina/cvs/tina-libs/tina/file/fileMatrix_r.c,v $
37 * Date : $Date: 2003/09/22 16:09:01 $
38 * Version : $Revision: 1.3 $
39 * CVS Id : $Id: fileMatrix_r.c,v 1.3 2003/09/22 16:09:01 tony Exp $
40 *
41 * Author : Legacy TINA
42 *
43 * Notes :
44 *
45 *
46 *********
47 */
48
49 #include "fileMatrix_r.h"
50
51 #if HAVE_CONFIG_H
52 #include <config.h>
53 #endif
54
55 #include <stdio.h>
56 #include <string.h>
57 #include <tina/sys/sysPro.h>
58 #include <tina/sys/sysDef.h>
59 #include <tina/math/mathPro.h>
60 #include <tina/math/mathDef.h>
61
62 Matrix *dmatrix_read_fp(FILE * fp, int n, int m)
63 {
64 int i, j;
65 Matrix *a;
66 Matrix *matrix_alloc();
67
68 if ((a = matrix_alloc(m, n, matrix_full, double_v)) == NULL)
69 return (NULL);
70 for (i = 0; i < a->m; i++)
71 for (j = 0; j < a->n; j++)
72 (void) fscanf(fp, "%lf ", &(a->el.double_v[i][j]));
73 return (a);
74 }
75
76 /**
77 * @brief allocate and read a full matrix. File contains, eg, "double full 44 44 ......"
78 */
79 Matrix* dmatrix_read (FILE* f)
80 {
81 Matrix* result;
82 char dtype[64], mtype[64];
83
84 fscanf(f, "%s %s", dtype, mtype);
85 if (strcmp(dtype, "double" ) != 0 || strcmp(mtype, "full" ) != 0)
86 {
87 error("Wrong data type in file", non_fatal);
88 return NULL;
89 }
90 {
91 int m,n;
92 fscanf(f, "%d %d", &m,&n);
93 result = dmatrix_read_fp (f, m, n);
94 }
95 return result;
96 }
97
98 /**
99 * @brief allocate and read a double matrix. File contains, eg "double 42 ...."
100 */
101 Vector* dvector_read (FILE* f)
102 {
103 Vector* result;
104 int i,n;
105 char type[64];
106
107 fscanf(f, "%s", type);
108 if (strcmp(type, "double") != 0)
109 {
110 error("Wrong data type in file", non_fatal);
111 return NULL;
112 }
113 fscanf (f, "%d", &n);
114 result = vector_alloc(n, double_v);
115
116 {
117 double *el = (double *) result->data;
118 for (i = 0; i < result->n; i++)
119 fscanf(f, "%lf", &el[i]);
120 }
121
122
123 return result;
124 }
125
126
This page was automatically generated by the
LXR engine.
Visit the LXR main site for more
information.