intmatmul.c File Reference

#include "stack-c.h"

Include dependency graph for intmatmul.c:

Go to the source code of this file.

Defines
#define	A(i, k)   a[i + k*n]
#define	B(k, j)   b[k + j*m]
#define	C(i, j)   c[i + j*n]
Functions
void matmul	__PARAMS ((double *a, int n, int m, double *b, int l, double *c))
int	intmatmul (char *fname)
void	matmul (a, int n, int m, b, int l, c)

---

Define Documentation

#define A	(	i,
		k		)	a[i + k*n]

Definition at line 47 of file intmatmul.c.

#define B	(	k,
		j		)	b[k + j*m]

Definition at line 48 of file intmatmul.c.

#define C	(	i,
		j		)	c[i + j*n]

Definition at line 49 of file intmatmul.c.

Referenced by C2F(), dgemm(), GetCom(), GetIMATRIX(), GetPOINTER(), GetSPARSE(), intex14c(), intsplin2d(), Real3::main(), Real1::main(), and matmul().

---

Function Documentation

void matmul __PARAMS

(

(double *a, int n, int m, double *b, int l, double *c)

)

int intmatmul

(

char *

fname

)

Definition at line 10 of file intmatmul.c.

References CheckLhs, CheckRhs, CreateVar, Error, GetRhsVar, LhsVar, m1, m2, matmul(), n1, n2, sciprint(), and stk.

{
  static int l1, m1, n1, l2, m2, n2, l3;
  static int minlhs=1, maxlhs=1, minrhs=2, maxrhs=2;

  /* Check number of inputs (rhs=2) and outputs (lhs=1) */
  CheckRhs(minrhs,maxrhs) ;
  CheckLhs(minlhs,maxlhs) ;

  /* Get A (#1) and B (#2) and create C (#3) as double ("d") matrices */
  GetRhsVar(1, "d", &m1, &n1, &l1);   /* m1, n1 (and l1) are output parameters */
  GetRhsVar(2, "d", &m2, &n2, &l2);   /* m1, n1 (and l1) are output parameters */
  CreateVar(3, "d", &m1, &n2, &l3);   /* m1 and n2 are input parameters */

  /* Check dimensions    */
  if (!(n1==m2)) { sciprint("%s: Incompatible inputs \r\n", "matmul");
    Error(999);
    return 0;}

  /* Call multiplication function 
   * inputs:stk(l1)->A, stk(l2)->B 
   *  output:stk(l3)->C               
   */
  matmul(stk(l1), m1, n1, stk(l2), n2, stk(l3));

  /*  Return C (3)  */
  LhsVar(1) = 3;
  return 0;
}

Here is the call graph for this function:

void matmul	(	a	,
		int	n,
		int	m,
		b	,
		int	l,
		c
	)

Definition at line 51 of file intmatmul.c.

References A, B, C, i, j, and s.

Referenced by intmatmul().

{
  int i,j,k;
  double s;
  for( i=0 ; i < n; i++){
    for( j=0; j < l; j++){
      s = 0.;
      for( k=0; k< m; k++){
        s += A(i,k)*B(k,j);
      }
      C(i,j) = s;
    }
  }
}

Here is the caller graph for this function:

---