~ [ source navigation ] ~ [ diff markup ] ~ [ identifier search ] ~ [ freetext search ] ~ [ file search ] ~

Linux Cross Reference
Tina6/tina-libs/tina/file/fileMatrix_r.c

Version: ~
Architecture: ~ [ i386 ] ~ [ alpha ] ~ [ m68k ] ~ [ mips ] ~ [ ppc ] ~ [ sparc ] ~ [ sparc64 ] ~

  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 

~ [ source navigation ] ~ [ diff markup ] ~ [ identifier search ] ~ [ freetext search ] ~ [ file search ] ~

This page was automatically generated by the LXR engine.
Visit the LXR main site for more information.