开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 查看源码
CALC.C
资源名称：MSDN_VC98.zip　[点击查看]
上传用户：bangxh
上传日期：2007-01-31
资源大小：42235k
文件大小：2k
源码类别：

Windows编程

开发平台：

Visual C++

CALC.C：源码内容
  2. /****************************************************************************
  3.                    Microsoft RPC Version 2.0
  4.            Copyright Microsoft Corp. 1992, 1993, 1994- 1996
  5.                         mandel Example
  7.     FILE:       calc.c
  9.     PURPOSE:    Server side of the RPC distributed application Mandel
  11.     FUNCTIONS:  MandelCalc() - Do the calculations for the Windows
  12.                                Mandelbrot Set distributed drawing program.
  14. ****************************************************************************/
  16. #include <stdlib.h>
  17. #include <stdio.h>
  18. #include <windows.h>
  20. #ifdef RPC
  21. #include "mdlrpc.h"
  22. #endif
  23. #include "mandel.h"
  26. short calcmand(double, double, short);
  29. void MandelCalc(PCPOINT    pcptLL,
  30.                 PLONGRECT  prcDraw,
  31.                 double     precision,
  32.                 DWORD      ulThreshold,
  33.                 LINEBUF *  pbBuf)
  34. {
  35.     DWORD   h, height, width;
  36.     double  dreal, dimag, dimag2;
  37.     short   maxit = 0;
  38.     short * pbPtr;
  40.     pbPtr = *pbBuf;   // LINEBUF is an array of shorts
  42.     dreal = pcptLL->real + ((double)prcDraw->xLeft * precision);
  43.     dimag = pcptLL->imag + ((double)prcDraw->yBottom * precision);
  45.     maxit = (short) ulThreshold;
  47.     height = (prcDraw->yTop - prcDraw->yBottom) + 1;
  48.     width = (prcDraw->xRight - prcDraw->xLeft) + 1;
  50.     for ( ; width > 0; --width, dreal += precision) {
  51.         for (dimag2 = dimag, h = height; h > 0; --h, dimag2 += precision) {
  52.             if ((dreal > 4.0) || (dreal < -4.0) ||
  53.                 (dimag2 > 4.0) || (dimag2 < -4.0))
  54.                 *(pbPtr++) = 0;
  55.             else
  56.                 *(pbPtr++) = calcmand(dreal, dimag2, maxit);
  57.         }
  58.     }
  59. }
  62. /* C version of the assembly language program */
  63. short calcmand(double dreal,
  64.                double dimag,
  65.                short  maxit)
  66. {
  67.     double x, y, xsq, ysq;
  68.     short k;
  70.     k = maxit;
  71.     x = dreal;
  72.     y = dimag;
  74.     while (1) {
  75.         xsq = x * x;
  76.         ysq = y * y;
  77.         y = 2.0 * x * y + dimag;
  78.         x = (xsq - ysq) + dreal;
  79.         if (--k == 0)
  80.             return((short) (maxit - k));
  81.         if ((xsq + ysq) > 4.0)
  82.             return((short) (maxit - k));
  83.     }
  84. }