File: property.c

package info (click to toggle)
felt 3.02-4
  • links: PTS
  • area: main
  • in suites: slink
  • size: 16,460 kB
  • ctags: 6,885
  • sloc: ansic: 72,103; fortran: 3,614; yacc: 2,825; lex: 1,172; sh: 311; makefile: 279
file content (73 lines) | stat: -rw-r--r-- 1,916 bytes parent folder | download | duplicates (3)
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
/*
    This file is part of the FElt finite element analysis package.
    Copyright (C) 1993-1997 Jason I. Gobat and Darren C. Atkinson

    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program; if not, write to the Free Software
    Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
*/

/************************************************************************
 * File:	property.c
 *	
 * Description:	
 *		
 ************************************************************************/

# include <stdio.h>
# include "matrix.h"

int IsSymmetricMatrix (a)	/* Aij == Aji ? 1 : 0		*/
   Matrix	a;		/* matrix to check for symmetry */
{
   unsigned	i, j;
   int		status;
   
   if (!IsSquare(a))
      return 0;

   status = 1;
   for (i = 1 ; i <= Mrows(a) ; i++) {
      for (j = i ; j <= Mcols(a) ; j++) {
         if (mdata(a, i, j) != mdata(a, j, i)) {
            status = 0;
            break;
         }
      }
      if (!status)
         break;
   }

   return status;
}

int IsZeroMatrix (a)		/* Aij == 0 ? 1 : 0		*/
   Matrix	a;		/* matrix to check for symmetry */
{
   unsigned	i, j;
   int		status;

   status = 1;
   for (i = 1 ; i <= Mrows(a) ; i++) {
      for (j = 1 ; j <= Mcols(a) ; j++) {
         if (mdata(a,i,j) != 0) {
            status = 0;
            break;
         }
      }
      if (!status)
         break;
   }

   return status;
}