Октава вместо bc

[vak]

Обычно я для мелких вычислений пользовался калькулятором bc. Командная строка всегда под рукой: сложить пару чисел или процент прикинуть. Но теперь всё больше переключаюсь на Octave. Двадцать первый век на дворе, всё-таки. У неё и оконный интерфейс есть, хотя он мне без надобности. А вот вектора, матрицы и комплексные числа иногда полезно иметь.



Простой алиас в шелле и подсвеченное приглашение превращают Октаву в удобную замену bc:

$ alias octave='/usr/bin/octave -q'
$ cat .octaverc
PS1('\[\e[0;32m\]\#>>\[\e[m\] ')

Comments

[vak]

Это я не пробовал.

[http://users.livejournal.com/_slw/]

а попробуйте?
если я не напортачил, то fftw3 сам дохнет. а тут интересно как?

[vak]

Да где ж мне столько данных набрать.

[winpooh.livejournal.com]

Из dev/random?

[http://users.livejournal.com/_slw/]

из /dev/random не понятно как результат валидировать.
тут вся ссуть эксперимента -- не будет ли мусор на выходе ну и сколько ресурсов будет потреблено.

[http://users.livejournal.com/_slw/]

for(int i=0; i<180M; i++)
printf("%f\n", sin(i*314159.26/60M));

как-нибудь так.

[vak]

Попробовал на массиве в 128*1024*1024 сэмплов - вроде работает, не глючит:

$ ./fourier.m
Creating 1:134217728 data vector...
Computing Fourier transform...
Done.
Time: 21 seconds
Memory footprint: 4074096 kbytes

Вот такой скрипт:

#! /usr/bin/octave -qf

n = 2^27;
step = 2.0 / n;
t = 0:step:2-step;

printf("Creating 1:%d data vector...\n", columns(t));
x = 2*sin(20*pi*t) + sin(100*pi*t);

printf("Computing Fourier transform...\n");
t0 = time();
y = fft(x);
t1 = time();
printf("Done.\n");

printf("Time: %d seconds\n", t1 - t0);
printf("Memory footprint: %d kbytes\n", getrusage().maxrss);

[http://users.livejournal.com/_slw/]

а в результат-то посмотрели?
а fftw3 отработать-то отработал, но такую хрень выдал наружу...
через одну частоты нулевый, максимальная частота -- нулевая и все в таком духе. может конечно и я чего намудрил.

[vak]

Попробовал на меньших размерах.
Результат похож на правду.

[http://users.livejournal.com/_slw/]

на сильно меньших? потому что вопрос-то именно о правде на такого размера данных. на меньших и fftw3 все хорошо делал