Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лабораторная работа №32.doc
Скачиваний:
12
Добавлен:
01.05.2014
Размер:
129.02 Кб
Скачать

Void main ()

{

double res;

double pi = 3.14159265359;

SAMPLE;//51

twopi = 2.0 / sqrt(pi);

lower = 0.0;

upper = 1;

if (upper >= 0.0)

{

res = simps(lower,upper,tol);

erf = twopi * res;

}

SAMPLE;//60

}

Отчет о результатах измерений для программы LABMET.EXE.

Создан программой Sampler ( версия от Feb 15 1999 )

1995-98 (C) сПбГэту, Мойсейчук Леонид

Список обработанных файлов.

----------------------------------------------------------------------

NN Имя обработанного файла

----------------------------------------------------------------------

1. LABMET.CPP

----------------------------------------------------------------------

Таблица с результатами измерений ( используется 2 из 416 записей )

----------------------------------------------------------------------

Исх.Поз. Прием.Поз. Общее время(мкс) Кол-во прох. Среднее время(мкс)

----------------------------------------------------------------------

1 : 51 1 : 60 71.24 1 71.24

---------------------------------------------------------------------

Профилирование оптимизированной программы интегрирования методом Симпсона

(измерение времен выполнения функциональных участков)

#include <iostream.h>

#include <conio.h>

#include <stdio.h>

#include <math.h>

#include "Sampler.h"

const double tol = 1.0E-4;

double sum, upper, lower, erf, twopi;

double fx(double x)

{

return exp(-x * x);

}

double simps(double lower,double upper,double tol)

{

int i;

double x,delta_x,even_sum,odd_sum,end_sum, sum1;

int pieces;

double sum;

pieces = 2;

delta_x = (upper - lower) / pieces;

odd_sum = fx(lower + delta_x);

even_sum = 0.0;

end_sum = fx(lower) + fx(upper);

sum = (end_sum + 4.0 * odd_sum) * delta_x / 3.0;

printf("%d", pieces);

printf("%f", sum);

do

{

pieces = pieces * 2;

sum1 = sum;

delta_x = (upper - lower) / pieces;

even_sum = even_sum + odd_sum;

odd_sum = 0.0;

for (i = 1; i <= (pieces / 2); i++)

{

x = lower + delta_x * (2.0 * i - 1.0);

odd_sum = odd_sum + fx(x);

}

sum = (end_sum + 4.0 * odd_sum + 2.0 * even_sum) * delta_x / 3.0;

}

while (abs(sum - sum1) > abs(tol * sum1));

return sum;

}

Void main ()

{

double res;

double pi = 3.14159265359;

SAMPLE;//51

twopi = 2.0 / sqrt(pi);

SAMPLE;//53

lower = 0.0;

SAMPLE;//55

upper = 1;

SAMPLE;//57

if (upper >= 0.0)

{

SAMPLE;//60

res = simps(lower,upper,tol);

SAMPLE;//62

erf = twopi * res;

SAMPLE;//64

}

SAMPLE;//62

}

Отчет о результатах измерений для программы LABMET.EXE.

Создан программой Sampler ( версия от Feb 15 1999 )

1995-98 (C) сПбГэту, Мойсейчук Леонид

Список обработанных файлов.

----------------------------------------------------------------------

NN Имя обработанного файла

----------------------------------------------------------------------

1. LABMET.CPP

----------------------------------------------------------------------

Таблица с результатами измерений ( используется 8 из 416 записей )

----------------------------------------------------------------------

Исх.Поз. Прием.Поз. Общее время(мкс) Кол-во прох. Среднее время(мкс)

----------------------------------------------------------------------

1 : 51 1 : 53 2.51 1 2.51

----------------------------------------------------------------------

1 : 53 1 : 55 0.84 1 0.84

----------------------------------------------------------------------

1 : 55 1 : 57 0.84 1 0.84

----------------------------------------------------------------------

1 : 57 1 : 60 0.84 1 0.84

----------------------------------------------------------------------

1 : 60 1 : 62 62.86 1 62.86

----------------------------------------------------------------------

1 : 62 1 : 64 1.68 1 1.68

----------------------------------------------------------------------

1 : 64 1 : 66 0.00 1 0.00

----------------------------------------------------------------------

Профилирование программы интегрирования методом Симпсона

(измерение времени полного выполнения программы)

На Паскале

uses sampler;

const szUnit : string = 'Plab1.pas';

const tol = 1.0E-4;

var sum,upper,lower,

erf,twopi : real;

function fx(x: real): real;

begin

fx:=exp(-x*x)

end;

procedure simps(lower,upper,tol : real; var sum : real);

var i : integer;

x,delta_x,even_sum,

odd_sum,end_sum,

sum1 : real;

pieces : integer;

begin

pieces:=2;

delta_x:=(upper-lower)/pieces;

odd_sum:=fx(lower+delta_x);

even_sum:=0.0;

end_sum:=fx(lower)+fx(upper);

sum:=(end_sum+4.0*odd_sum)*delta_x/3.0;

writeln(pieces:5,sum);

repeat

pieces:=pieces*2;

sum1:=sum;

delta_x:=(upper-lower)/pieces;

even_sum:=even_sum+odd_sum;

odd_sum:=0.0;

for i:=1 to pieces div 2 do

begin

x:=lower+delta_x*(2.0*i-1.0);

odd_sum:=odd_sum+fx(x)

end;

sum:=(end_sum+4.0*odd_sum+2.0*even_sum)*delta_x/3.0;

until abs(sum-sum1)<=abs(tol*sum1)

end;

begin

Sample(szUnit,1);

twopi := 2.0/sqrt(pi);

lower := 0.0;

upper := 1;

writeln;

writeln('Erf? ');

if upper >= 0.0 then

begin

simps(lower,upper,tol,sum);

erf:=twopi*sum;

writeln('Erf of ',upper:7:2,', is ',erf:12:8)

end;

Sample(szUnit,2);

end.

Отчет о результатах измерений для программы PLAB1.EXE.

Создан программой Sampler ( версия от Feb 15 1999 )