/* SC.f -- translated by f2c (version 19960827).
   You must link the resulting object file with the libraries:
	-lf2c -lm   (in that order)
*/

#include "Data_f2c.h"
#ifdef WNT
#include <MathBase.h>
#else 
#define  __MathBase_API
#endif
__MathBase_API doublereal msc_(ndimen, vecte1, vecte2)
integer *ndimen;
doublereal *vecte1, *vecte2;
{
    /* System generated locals */
    integer i__1;
    doublereal ret_val;

    /* Local variables */
    static integer i__;
    static doublereal x;



/* < */
/* **NOTICE */
/*  THIS SOFTWARE IS THE PROPERTY OF CISIGRAPH. */
/*  THIS CODE MUST NOT BE DISTRIBUTED OR COPIED WITHOUT THE PRIOR */
/*  WRITTEN PERMISSION OF CISIGRAPH AND IS ONLY TO BE USED ON THE */
/*  SITE WHERE IT IS INSTALLED BY CISIGRAPH */
/* **NOTICE */

/************************************************************************
*******/

/*     FONCTION : */
/*     ---------- */
/*        Calcul du produit scalaire de 2 vecteurs dans l' espace */
/*        de dimension NDIMEN. */

/*     MOTS CLES : */
/*     ----------- */
/*        PRODUIT MSCALAIRE. */

/*     ARGUMENTS D'ENTREE : */
/*     ------------------ */
/*        NDIMEN : Dimension de l' espace. */
/*        VECTE1,VECTE2: Les vecteurs. */

/*     ARGUMENTS DE SORTIE : */
/*     ------------------- */

/*     COMMONS UTILISES   : */
/*     ---------------- */

/*     REFERENCES APPELEES   : */
/*     ----------------------- */

/*     DESCRIPTION/REMARQUES/LIMITATIONS : */
/*     ----------------------------------- */

/* $    HISTORIQUE DES MODIFICATIONS   : */
/*     -------------------------------- */
/*        18-07-1988: RBD ; Changement de nom des arguments pour plus */
/*                          de lisibilite. */
/*        16-01-1987: Verification implicite NDIMEN >0 JJM. */

/* > */
/* ***********************************************************************
 */


/*     PRODUIT MSCALAIRE */
    /* Parameter adjustments */
    --vecte2;
    --vecte1;

    /* Function Body */
    x = 0.;

    i__1 = *ndimen;
    for (i__ = 1; i__ <= i__1; ++i__) {
	x += vecte1[i__] * vecte2[i__];
/* L100: */
    }
    ret_val = x;

/* ----------------------------------- THE END -------------------------- 
*/

    return ret_val;
} /* msc_ */

