multiply.c

很好的一个约束遗传算法优化程序

#include "genocop.h"

/********************************************************************************/
/*                                                                              */
/*           FUNCTION NAME     :   mmprod()                                     */
/*                                                                              */
/*           SYNOPSIS          :   void mmprod(m,nm,n,mul_cm,mul_am,mul_bm)     */
/*                                                                              */
/*           DESCRIPTION       :   This function finds the product of two       */
/*                                  float matrices                              */
/*                                                                              */
/*           FUNCTIONS CALLED  :   None                                         */
/*                                                                              */
/*           CALLING FUNCITONS :   find_org_in_eq(),                            */
/*                                 main().                                      */
/*                                                                              */
/*           AUTHOR            :   Swarnalatha Swaminathan                      */
/*                                                                              */
/*           DATE              :   1/17/92                                      */
/*                                                                              */
/*                                                                              */
/*           REV            DATE            BY           DESCRIPTION            */
/*           ---            ----            --           -----------            */
/*                                                                              */
/*                                                                              */
/********************************************************************************/




void 
mmprod(m, nm, n, mul_cm, mul_am, mul_bm)
    int             m,			    /* row of mul_am matrix */
                    nm,			    /* column of mul_am and row of
					     * mul_bm matrices */
                    n;			    /* row of mul_bm matrix */
    MATRIX          mul_cm,		    /* the final matrix */
                    mul_am,		    /* the first matrix to be
					     * multiplied */
                    mul_bm;		    /* the second matrix to be
					     * multiplied */
{
    int             i, j, k;		    /* counter variables, where i=m,
					     * j=nm, k=n */

    for (i = 1; i <= m; i++)
	for (j = 1; j <= n; j++)
	{
	    mul_cm[i][j] = 0.0;
	    for (k = 1; k < nm + 1; k++)
		mul_cm[i][j] = mul_am[i][k] * mul_bm[k][j] + mul_cm[i][j];
	}
}



/********************************************************************************/
/*                                                                              */
/*           FUNCTION NAME     :   mvprod()                                     */
/*                                                                              */
/*           SYNOPSIS          :   void mvprod(m,nm,cm,am,bm)                   */
/*                                                                              */
/*           DESCRIPTION       :   This function finds the product of a float   */
/*                                  vector and a float matrix                   */
/*                                                                              */
/*           FUNCTIONS CALLED  :   None                                         */
/*                                                                              */
/*           CALLING FUNCITONS :   find_org_in_eq(),                            */
/*                                 main().                                      */
/*                                                                              */
/*           AUTHOR            :   Swarnalatha Swaminathan                      */
/*                                                                              */
/*           DATE              :   1/17/92                                      */
/*                                                                              */
/*                                                                              */
/*           REV            DATE            BY           DESCRIPTION            */
/*           ---            ----            --           -----------            */
/*                                                                              */
/*                                                                              */
/********************************************************************************/




void 
mvprod(m, nm, cm, am, bm)
    int             m,			    /* row of matrix am and length of
					     * vector cm */
                    nm;			    /* column of matrix am */
    VECTOR          cm,			    /* the final vector */
                    bm;			    /* original vector */
    MATRIX          am;			    /* original matrix */
{
    int             i, k;		    /* counter variables, where i=m,
					     * k=nm */

    for (i = 1; i <= m; i++)
    {
	cm[i] = 0.0;
	for (k = 1; k <= nm; k++)
	    cm[i] = cm[i] + am[i][k] * bm[k];
    }
}