VC++用MKL产生正态分布随机数、计算样本均值与标准差源程序

// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。
//

#include “stdafx.h“
#include 
#include 
#include 
#include 
// Main program

int _tmain（int argc _TCHAR* argv[]）
{
   double r[1000]; 
   double ss1s2Raw2MomCentral2Mom; 
   VSLStreamStatePtr stream;
   VSLSSTaskPtr task;
   MKL_INT i j errstatusxstoragestatusDIMN;

   s = 0.0;
   vslNewStream（ &stream VSL_BRNG_MCG31 777 ）;
  
      vdRngGaussian（ VSL_RNG_METHOD_GAUSSIAN_ICDF stream 1000 r 5.0 2.0 ）;
      for （ j=0; j<1000; j++ ） {
         s += r[j];
      }
   s /= 1000.0;
   errstatus = vslSaveStreamF（ stream “d:\\abc“ ）;
   if（（errstatus==VSL_ERROR_OK）||（errstatus==VSL_STATUS_OK））
   	   printf（“\n stored successd\n“）;
   else
   	   printf（“\n err no. : %d stored no successd\n“errstatus）;
    
  
   vslDeleteStream（ &stream ）;
  
    printf（ “1.Sample mean of normal distribution = %f\n“ s ）;

	DIM=1;
	N=1000;
	xstorage=VSL_SS_MATRIX_STORAGE_ROWS;

	/* Create a task */     
	status = vsldSSNewTask（ &task &DIM &N &xstorage r 0 0 ）;       

	/* Initialize the task parameters */    

	status = vsldSSEditTask（ task VSL_SS_ED_VARIATION &s2 ）;        

	status = vsldSSEditMoments（ task &s1 &Raw2Mom 0 0&Central2Mom 0 0 ）; 
	status = vsldSSCompute（task  VSL_SS_MEAN | VSL_SS_VARIATION  VSL_SS_METHOD_FAST）;

	printf（“status3=%d\n“status）;
    printf（ “3.Sample mean and stdev of normal distribution = （%lf %lf）\n“ ss2*s1*sqrt（N/（N-1.0）） ）;

    status=vslSSDeleteTask （&task）;
	return 0;
}


/*

//*******************************************************************************
!  Copyright（C） 2003-2014 Intel Corporation. All Rights Reserved.
!
!  The source code information  and  material （“Material“） contained herein is
!  owned  by Intel Corporation or its suppliers or licensors and title to such
!  Material remains  with Intel Corporation  or its suppliers or licensors. The
!  Material  contains proprietary information  of  Intel or  its  suppliers and
!  licensors. The  Material is protected by worldwide copyright laws and treaty
!  provisions. No  part  of  the  Material  may  be  used  copied reproduced
!  modified published uploaded posted transmitted distributed or disclosed
!  in any way  without Intel‘s  prior  express written  permission. No  license
!  under  any patent copyright  or  other intellectual property rights  in the
!  Material  is  granted  to  or  conferred  upon  you  either  expressly  by
!  implication inducement  estoppel or  otherwise.  Any  license  under  such
!  intellectual  property  rights must  be express  and  approved  by  Intel in
!  writing.
!
!  *Third Party trademarks are the property of their respective owners.
!
!  Unless otherwise  agreed  by Intel  in writing you may not remove  or alter
!  this  notice or  any other notice embedded  in Materials by Intel or Intel‘s
!  suppliers or licensors