Rendering Text with Direct2D & DirectWrite

#include "stdafx.h"

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

#include "hnum_f2c.h"
namespace harlinn
{
    namespace numerics
    {
        namespace SuperLU
        {
                doublereal ddot_(integer *n, doublereal *dx, integer *incx, doublereal *dy, integer *incy)
                {


                    /* System generated locals */
                    integer i__1;
                    doublereal ret_val;

                    /* Local variables */
                    static integer i, m;
                    static doublereal dtemp;
                    static integer ix, iy, mp1;


                /*     forms the dot product of two vectors.   
                       uses unrolled loops for increments equal to one.   
                       jack dongarra, linpack, 3/11/78.   
                       modified 12/3/93, array(1) declarations changed to array(*)   


    
                   Parameter adjustments   
                       Function Body */
                #define DY(I) dy[(I)-1]
                #define DX(I) dx[(I)-1]


                    ret_val = 0.;
                    dtemp = 0.;
                    if (*n <= 0) {
	                return ret_val;
                    }
                    if (*incx == 1 && *incy == 1) {
	                goto L20;
                    }

                /*        code for unequal increments or equal increments   
                            not equal to 1 */

                    ix = 1;
                    iy = 1;
                    if (*incx < 0) {
	                ix = (-(*n) + 1) * *incx + 1;
                    }
                    if (*incy < 0) {
	                iy = (-(*n) + 1) * *incy + 1;
                    }
                    i__1 = *n;
                    for (i = 1; i <= *n; ++i) {
	                dtemp += DX(ix) * DY(iy);
	                ix += *incx;
	                iy += *incy;
                /* L10: */
                    }
                    ret_val = dtemp;
                    return ret_val;

                /*        code for both increments equal to 1   


                          clean-up loop */

                L20:
                    m = *n % 5;
                    if (m == 0) {
	                goto L40;
                    }
                    i__1 = m;
                    for (i = 1; i <= m; ++i) {
	                dtemp += DX(i) * DY(i);
                /* L30: */
                    }
                    if (*n < 5) {
	                goto L60;
                    }
                L40:
                    mp1 = m + 1;
                    i__1 = *n;
                    for (i = mp1; i <= *n; i += 5) {
	                dtemp = dtemp + DX(i) * DY(i) + DX(i + 1) * DY(i + 1) + DX(i + 2) * 
		                DY(i + 2) + DX(i + 3) * DY(i + 3) + DX(i + 4) * DY(i + 4);
                /* L50: */
                    }
                L60:
                    ret_val = dtemp;
                    return ret_val;
                } /* ddot_ */

        };
    };
};