FreePascal.ru

Sergei I. Gorelkin писал(а):если бы все extended и превращались в double, первоначальная проблема вроде не должна возникать...

это не происходит и проблема возникает...

Table 3.4: Supported Real types

Type Range Significant digits Size
Real platform dependant ??? 4 or 8
Single 1.5E-45 .. 3.4E38 7-8 4
Double 5.0E-324 .. 1.7E308 15-16 8
Extended 1.9E-4932 .. 1.1E4932 19-20 10

Спасибо, за интерес!

Выяснил, что проблема в несогласованности типа переменных. Действительно, function ln(d: double) : double, a function Random : extended. Поэтому в ln(1-Random) под ln может оказаться ноль в «double измерении».

Как один из вариантов выхода я поставил проверку на ноль. Может ...

Как известно, функция ln не определена в нуле.

При достаточно большом количестве итераций в строке ln(1-Random) или ln(Random) появляется «Runtime error 200». Из чего можно сделать вывод, что Random включает как 0 так и 1.

Однако, дело может быть и в типе переменных. Например, Random – Extended, а ...

Никто не замечал, что в fpc-2.4.0.x86_64-linux функция Random генерирует случайное число из (0;1], тогда как обычно принято из [0;1).

сайт разработчиков не отвечает...

Что бы это значило?

shade писал(а):в фортране двуменый массив храниться столбцами, а не строками

Совершенно верно

shade писал(а):Так что в FPC нужно оптимизировать адресную арифметику...

вот и задача для программистов появилась

программа:
Код:

uses Dos;
var
i,j,k,l,n : integer;
a,b,c : array [1..2000,1..2000] of double;
col,rol : array [1..2000] of double;
op,mf,s : double;
tim1,tim2,tim : real;

function second :real;
var
h1,m1,s1,ds1 : word;
begin
GetTime(h1,m1,s1,ds1);
second := h1*3600+m1*60+s1+0.01*ds1 ...

Рождённый_в_СССР писал(а):не вычисляя всё, а вычисляя только ту часть, которую потом нужно выводить

А такое бывает? Я сильно сомневаюсь....

Инетесно, заменил
DO 2 I = 1,N
DO L = 1,N
ROW(L) = A(I,L)
END DO
DO 2 J = 1,N
C(I,J) =0.0d0
DO 3 K = 1,N
3 C(I,J) = C(I,J) + ROW(K)*B(K,J)
2 CONTINUE
на
DO 2 I = 1,N
DO 2 J = 1,N
C(I,J) =0.0d0
DO 3 K = 1,N
3 C(I,J ...

volvo877 писал(а):И не надо говорить того, о чем не знаешь.

А Вы что ли знаете

Turbo Pascal Version 8.0
Borland International

ev писал(а):потому что фортран уже медленно вымирает

Я бы так не сказал, а вот действительно мощная контора Borland почему то прекратила поддержку Паскаля.

вряд ли найдёшь сейчас даже несколько школ где ему учат...

А зря, я сам был испорчен в школе и универе Паскалем, хотя может это и логично: ведь Паскаль разрабатывался в своё время наряду с Бейсиком как учебный язык.
Был сегодня в книжном всё - завалено Паскалем, Си, Бейсиком и др. по Фортрану ...

то любой компилятор для вас также является черным ящиком

Согласен, но этот ящик, по краней мере, на уровень ниже, чем пакеты MatLAB, MatCAD, Mathemetica и прочее

я знаю много людей (занимающихся подобными задачами) и ни одного кто бы это делал на фортране - паскаль, си или маткад/подобные ...

Brainenjii писал(а):G95 несколько уступает fpc

Возможно это связано со спецификой задачи.