интерактивный курс: введение в программирование на PARI/GP

[Dmitriy40]

arseniiv
Я бы сказал знак ";" в определённых местах наоборот подавляет вывод результата выражения. Например:

? v=12
%1 = 12
? v=15;
? v=22
%3 = 22
?

Отображено в консоли лишь 12 и 22, но не 15. Хотя результат 15 и запомнен и доступен для ссылок.
Но далеко не везде можно пропускать ";".

-- 03.10.2015, 23:06 --

Begemot82
Спасибо.

[arseniiv]

Dmitriy40 в сообщении #1058928 писал(а):

Я бы сказал знак ";" в определённых местах наоборот подавляет вывод результата выражения.

А, такой эффект в Wolfram тоже есть, потому что там в таком случае a; ; b; считается синонимом a; Null; b; Null, и Null не выводится, если получается в результате чего-нибудь. Может, и здесь снова так?

[Dmitriy40]

arseniiv
Я детально в синтаксис не вникал (вон даже про комментарии не знал), про подавление вывода от наличия ";" заметил где-то в хелпе случайно.

[arseniiv]

(Оффтоп)

Скачал, 2; 3 возвращает 3, хотя print(;2;) почему-то печатает 2. Ладно, любопытство удовлетворил.

[Begemot82]

Обнаружил интересную особенность - большие числа можно записывать через пробелы

a=123 456 789;b=38734 36777;print(a" "b);

[RIP]

Begemot82 в сообщении #1060324 писал(а):

Обнаружил интересную особенность - большие числа можно записывать через пробелы

PARI/GP вообще игнорирует пробелы (за исключением тех, что в кавычках):

gp > s    i    n     (   P         i    /  2  .  0 )    
%1 = 1.0000000000000000000000000000000000000

[VAL]

PARI тормозит
Летом считал на PARI. И он приятно удивил быстродействием. В десятки раз быстрее maple.
Сейчас запустил аналогичную задачу. Считает раз в 10 медленнее maple. Т.е. в сотни раз медленнее себя летнего.
На компе ничего не переставлял.
Для проверки запустил не аналогичную, а ТУ ЖЕ задачу. Тормозит.

Попробовал выделить побольше памяти, увеличить primelimit. Эффекта нет.

В чем может быть дело?

[NT2000]

VAL
А что у вас за ОС?
Посмотрите в процессах, чем занят компьютер?

PS. Бывает, что антивир страшно подвисает комп.
Если у вас Винда, для проверки на время вычислений отключите Интернет и антивирус.

[VAL]

NT2000 в сообщении #1100533 писал(а):

VAL
А что у вас за ОС?
Посмотрите в процессах, чем занят компьютер?

PS. Бывает, что антивир страшно подвисает комп.
Если у вас Винда, для проверки на время вычислений отключите Интернет и антивирус.

ОС - семерка.
Антивирус - Avast. В отличие от Касперского в торможении программ не замечен.
Процессор загружен практически только одним PARI/GP. На 25%.
Кстати!
Может быть в этом и разгадка. Возможно, (что-то такое припоминаю) я летом каким-то образом (а каким, уже не помню ) подключал к работе PARI все (или почти все) ядра. А сейчас используется всего одно.

Тогда вопрос дилетанта: как заставить работать с PARI другие ядра? Процессор - Intel Core 5 i-32

[Dmitriy40]

(Оффтоп)

Очень мне интересно, на какой же это задаче увеличение количества ядер с 1 до 4 даёт 100-кратный выигрыш производительности?!

Команда \d в консоли выдаёт список значений по умолчанию множества параметров , среди них есть два вроде бы относящиеся к многозадачности: nbthreads и threadsize. Можно покопать в их сторону.

[VAL]

Dmitriy40 в сообщении #1100543 писал(а):

Очень мне интересно, на какой же это задаче увеличение количества ядер с 1 до 4 даёт 100-кратный выигрыш производительности?!

Ну... менее эмоциональная оценка, свидетельствует о замедлении в разы, а не в сотни раз

Цитата:

Команда \d в консоли выдаёт список значений по умолчанию множества параметров , среди них есть два вроде бы относящиеся к многозадачности: nbthreads и threadsize. Можно покопать в их сторону.

Попробовал изменить значение nbthreads с 1 до 3 (по аналогии с тем, как изменяется primelimit). Но ничего не изменилось, по-прежнему процессор загружен на 25%

[dmd]

Наткнулся на ошибку, не пойму - у меня или в pari/gp. Простейший код, удаление из списка элемента

tst()=
{

local(H, S):list;

H= List(); S= List();

for(i=1, 100, listput(S, i));

H= S;

print("#H= ", #H);
print("#S= ", #S);

listpop(H, 67);

print("#H= ", #H);
print("#S= ", #S);

};

просто в gp работает правильно, но в gp2c оператор listpop(H, 67) удаляет и из H и из S:

[dmd]

Два варианта посчитать установленные биты:

? bitcount1(n)= sum(i=0, #binary(n), bittest(n, i));
?
? bitcount2(n)={i=0;while(n!=0,i++;n=bitand(n,n-1));return(i)};
?
? for(i=1,10^6,bitcount1(i))
time = 4,660 ms.
?
? for(i=1,10^6,bitcount2(i))
time = 10,266 ms.

Можно ли сделать ещё быстрее?

[Dmitriy40]

Сравнение скорости не совсем корректно, для других диапазонов чисел результат может быть более драматичным, например:

? for(i=10^12, 10^12+10^6, bitcount1(i))
time = 2,555 ms.
? for(i=10^12, 10^12+10^6, bitcount2(i))
time = 1 min, 20,199 ms.

Возможное улучшение:

? bitcount3(n) = vecsum(binary(n));
? for(i=10^12, 10^12+10^6, bitcount3(i))
time = 776 ms.

[dmd]

Здорово, спасибо!


Другой вопрос. Как быстро пробежать по установленным битам в большом числе, допустим размером тысяча бит?

?
?\\ N - длина в битах числа n
?
?bitdel(n, i)= bitnegimply(n, 2^i);
?
?for(i=1, N, if(bittest(n, i), if(<некое условие>, n= bitdel(n, i))));
?

Возможно такое ускорить?