сестра таланта
Mar. 20th, 2013 01:31 am![[personal profile]](https://www.dreamwidth.org/img/silk/identity/user.png){kind=small}
avvaНаткнулся на этот исторический документ, интерпретатор игрушечного фрагмента языка J, написанный в 89-м году создателем родственного языка K.
Если вы не знаете J/APL (как и я; давно хочу изучить, да все как-то...), то предлагаю вам задачу, которую я поставил перед собой из чистоголюбопытства мазохизма упрямства: прочитайте и поймите эту программу (включая понимание того, как устроен язык, который она интерпретирует) только с листа, не перенося ее в компьютер, не запуская, не переписывая исходники более понятным образом итп.
typedef char C;typedef long I;
typedef struct a{I t,r,d[3],p[2];}*A;
#define P printf
#define R return
#define V1(f) A f(w)A w;
#define V2(f) A f(a,w)A a,w;
#define DO(n,x) {I i=0,_n=(n);for(;i<_n;++i){x;}}
I *ma(n){R(I*)malloc(n*4);}mv(d,s,n)I *d,*s;{DO(n,d[i]=s[i]);}
tr(r,d)I *d;{I z=1;DO(r,z=z*d[i]);R z;}
A ga(t,r,d)I *d;{A z=(A)ma(5+tr(r,d));z->t=t,z->r=r,mv(z->d,d,r);R z;}
V1(iota){I n=*w->p;A z=ga(0,1,&n);DO(n,z->p[i]=i);R z;}
V2(plus){I r=w->r,*d=w->d,n=tr(r,d);A z=ga(0,r,d);
DO(n,z->p[i]=a->p[i]+w->p[i]);R z;}
V2(from){I r=w->r-1,*d=w->d+1,n=tr(r,d);
A z=ga(w->t,r,d);mv(z->p,w->p+(n**a->p),n);R z;}
V1(box){A z=ga(1,0,0);*z->p=(I)w;R z;}
V2(cat){I an=tr(a->r,a->d),wn=tr(w->r,w->d),n=an+wn;
A z=ga(w->t,1,&n);mv(z->p,a->p,an);mv(z->p+an,w->p,wn);R z;}
V2(find){}
V2(rsh){I r=a->r?*a->d:1,n=tr(r,a->p),wn=tr(w->r,w->d);
A z=ga(w->t,r,a->p);mv(z->p,w->p,wn=n>wn?wn:n);
if(n-=wn)mv(z->p+wn,z->p,n);R z;}
V1(sha){A z=ga(0,1,&w->r);mv(z->p,w->d,w->r);R z;}
V1(id){R w;}V1(size){A z=ga(0,0,0);*z->p=w->r?*w->d:1;R z;}
pi(i){P("%d ",i);}nl(){P("\n");}
pr(w)A w;{I r=w->r,*d=w->d,n=tr(r,d);DO(r,pi(d[i]));nl();
if(w->t)DO(n,P("< ");pr(w->p[i]))else DO(n,pi(w->p[i]));nl();}
C vt[]="+{~<#,";
A(*vd[])()={0,plus,from,find,0,rsh,cat},
(*vm[])()={0,id,size,iota,box,sha,0};
I st[26]; qp(a){R a>='a'&&a<='z';}qv(a){R a<'a';}
A ex(e)I *e;{I a=*e;
if(qp(a)){if(e[1]=='=')R st[a-'a']=ex(e+2);a= st[ a-'a'];}
R qv(a)?(*vm[a])(ex(e+1)):e[1]?(*vd[e[1]])(a,ex(e+2)):(A)a;}
noun(c){A z;if(c<'0'||c>'9')R 0;z=ga(0,0,0);*z->p=c-'0';R z;}
verb(c){I i=0;for(;vt[i];)if(vt[i++]==c)R i;R 0;}
I *wd(s)C *s;{I a,n=strlen(s),*e=ma(n+1);C c;
DO(n,e[i]=(a=noun(c=s[i]))?a:(a=verb(c))?a:c);e[n]=0;R e;}
main(){C s[99];while(gets(s))pr(ex(wd(s)));}
Я потратил на это немало времени сегодня вечером, но, кажется, добился.
Но черт побери. Черт побери! (подавляет в себе сложные эмоции и уходит, качая головой)
Если вы не знаете J/APL (как и я; давно хочу изучить, да все как-то...), то предлагаю вам задачу, которую я поставил перед собой из чистого
typedef char C;typedef long I;
typedef struct a{I t,r,d[3],p[2];}*A;
#define P printf
#define R return
#define V1(f) A f(w)A w;
#define V2(f) A f(a,w)A a,w;
#define DO(n,x) {I i=0,_n=(n);for(;i<_n;++i){x;}}
I *ma(n){R(I*)malloc(n*4);}mv(d,s,n)I *d,*s;{DO(n,d[i]=s[i]);}
tr(r,d)I *d;{I z=1;DO(r,z=z*d[i]);R z;}
A ga(t,r,d)I *d;{A z=(A)ma(5+tr(r,d));z->t=t,z->r=r,mv(z->d,d,r);R z;}
V1(iota){I n=*w->p;A z=ga(0,1,&n);DO(n,z->p[i]=i);R z;}
V2(plus){I r=w->r,*d=w->d,n=tr(r,d);A z=ga(0,r,d);
DO(n,z->p[i]=a->p[i]+w->p[i]);R z;}
V2(from){I r=w->r-1,*d=w->d+1,n=tr(r,d);
A z=ga(w->t,r,d);mv(z->p,w->p+(n**a->p),n);R z;}
V1(box){A z=ga(1,0,0);*z->p=(I)w;R z;}
V2(cat){I an=tr(a->r,a->d),wn=tr(w->r,w->d),n=an+wn;
A z=ga(w->t,1,&n);mv(z->p,a->p,an);mv(z->p+an,w->p,wn);R z;}
V2(find){}
V2(rsh){I r=a->r?*a->d:1,n=tr(r,a->p),wn=tr(w->r,w->d);
A z=ga(w->t,r,a->p);mv(z->p,w->p,wn=n>wn?wn:n);
if(n-=wn)mv(z->p+wn,z->p,n);R z;}
V1(sha){A z=ga(0,1,&w->r);mv(z->p,w->d,w->r);R z;}
V1(id){R w;}V1(size){A z=ga(0,0,0);*z->p=w->r?*w->d:1;R z;}
pi(i){P("%d ",i);}nl(){P("\n");}
pr(w)A w;{I r=w->r,*d=w->d,n=tr(r,d);DO(r,pi(d[i]));nl();
if(w->t)DO(n,P("< ");pr(w->p[i]))else DO(n,pi(w->p[i]));nl();}
C vt[]="+{~<#,";
A(*vd[])()={0,plus,from,find,0,rsh,cat},
(*vm[])()={0,id,size,iota,box,sha,0};
I st[26]; qp(a){R a>='a'&&a<='z';}qv(a){R a<'a';}
A ex(e)I *e;{I a=*e;
if(qp(a)){if(e[1]=='=')R st[a-'a']=ex(e+2);a= st[ a-'a'];}
R qv(a)?(*vm[a])(ex(e+1)):e[1]?(*vd[e[1]])(a,ex(e+2)):(A)a;}
noun(c){A z;if(c<'0'||c>'9')R 0;z=ga(0,0,0);*z->p=c-'0';R z;}
verb(c){I i=0;for(;vt[i];)if(vt[i++]==c)R i;R 0;}
I *wd(s)C *s;{I a,n=strlen(s),*e=ma(n+1);C c;
DO(n,e[i]=(a=noun(c=s[i]))?a:(a=verb(c))?a:c);e[n]=0;R e;}
main(){C s[99];while(gets(s))pr(ex(wd(s)));}
Я потратил на это немало времени сегодня вечером, но, кажется, добился.
Но черт побери. Черт побери! (подавляет в себе сложные эмоции и уходит, качая головой)
![[identity profile]](https://www.dreamwidth.org/img/silk/identity/openid.png){kind=small}
no subject
Date: 2013-03-19 11:45 pm (UTC)первая: "Усердие всё превозмогает!"
вторая: "Бывает, что усердие превозмогает и рассудок."
no subject
Date: 2013-03-19 11:56 pm (UTC)Воплощается в трех лицах.
Но все трое, если нужно,
Воедино могут слиться;
Что постигнуть это чудо
И поверить не на шутку
Может только тот, кто бросит
Вызов здравому рассудку;
no subject
Date: 2013-03-20 12:06 am (UTC)Стоить заметить, что программы на J часто выглядят куда менее понятными. Особенно для непосвященных. Тем более, что их всегда можно "свернуть" в одну строчку оператором компиляции. Классика жанра -- программа, реализующая сиплекс-метод:
[-(({~"1<:@{:)-(i.@#@[=<:@{.@]))*/(({~<:@{.)%({~<@<:))>
(отсюда (http://www.jsoftware.com/jwiki/Stories/RichardBrown))
no subject
Date: 2013-03-20 06:43 am (UTC)no subject
Date: 2013-03-20 01:06 am (UTC)Разрешается делать заметки карандашом или нужно все понять, смотря на экран? Если второй вариант, то сразу сдаюсь :)
no subject
Date: 2013-03-20 01:41 am (UTC)no subject
Date: 2013-03-20 02:02 am (UTC)no subject
Date: 2013-03-20 02:49 am (UTC)no subject
Date: 2013-03-20 03:25 am (UTC)Общий исток APL - вряд ли для S (R). Люди разные, места разные, время создания примерно одно и тоже...
Насчет "очень близкий синтаксис"... Вот, например, одна и та же по смыслу программа на R и J (http://progopedia.ru/example/fibonacci/)
R:
fib <- function(n) { if (n < 2) n else fib(n - 1) + fib(n - 2) } print.table(lapply(1:16, fib))J:
no subject
Date: 2013-03-20 01:20 pm (UTC)no subject
Date: 2013-03-20 03:11 pm (UTC)Однако, даже к вашему примеру есть куча вопросов. Векторный язык, а функцию к вектору не применяете? Как так? Вместо этого какой-то lapply, который потом требует какой-то print.table.
То есть, то что "программа на нём выглядит совершенно понятной", как написал migmit - это обманчиво.
fib = function(n){ifelse(n < 2, 1, fib(n-1)+fib(n-2))}
paste(1:16, fib(1:16-1), sep=": ")
no subject
Date: 2013-03-23 11:31 am (UTC)fib <- function(n){ # объявить функцию ifelse(n < 2, 1, fib(n-1)+fib(n-2)) # применить рекурсивно } paste(1:16, # склеить вектор 1-16 fib(1:16-1), # с результатами функции sep=": ") # разделитель :_no subject
Date: 2013-04-04 01:02 am (UTC)Во-первых, искренне, спасибо за ваши посты в R-community и за хороший код.
Во-вторых, спасибо за комментирование/форматирование моего кода. Ах, если б R код, который пишут учьоные, был так оформлен. (Кстатте, "1:16-1" тоже заслуживает комментария)
Но. Не-спасибо за не-комментирование кода kit1980ukr, который даже простой пример написал кое-как.
И отдельное не-спасибо за замену = на <- в моем коде.
no subject
Date: 2013-04-04 06:31 am (UTC)Вы конечно не поверите, но некие флюиды витали :)
Да, логика алгоритма "размазана" между двумя операторами. Я видел это, но это уже не комментарий кода. Это "шапка" описание, или чунк (как там было? :) "всякого ....?).
PS в конце концов именно так написав автор кода делает код более понятным :) R это вполне себе язык, на котором именно описывают что было сделано с данными (а не "портянки кода" пишут ()хотя если хочет человек, то портянку кода он на чем угодно напишет :))
"
1)"учьоных" старательно приучают пользовать for() ("они же начинающие") 2) и не редактировать свой код в emacs (там это практически автоматически получается именно в таком виде)
за исключением когда оператор присваивания не пишут в стандартной для R форме :)
Кроме того a<-b , a->b ; и для того что бы красиво (компактно) оформить код (именно использовав APL-стиля операторы) очень даже может пригодится.
no subject
Date: 2013-03-20 11:41 pm (UTC)fib =: 3 : '{. (}.,+/)^:y i.2'
Смысл (для тех кто не читает J) - взять пару чисел, и заменить их на пару из хвоста и суммы чисел. Проделать это N раз, начиная с (0 1), и в конце взять голову пары.
no subject
Date: 2013-03-24 07:19 pm (UTC)A вы знаете язык лучше R справляющийся с задачами анали
Date: 2013-03-20 02:49 am (UTC)Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-20 03:16 am (UTC)И что значит "лучше справляющийся"? Ну, допустим, я скажу Питон, или не дай б-г Джава лучше справляются. Что вы возразите - что там нет rnorm и ggplot2? И при чем здесь язык?
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-20 09:48 am (UTC)Потому что обычно все что я вижу, гораздо легче и понятнее выражается там (и мне и насколько я понимаю аналитикам), а фишек которые дейтсвительно нужны от мейнстримных языков программирования (типа возможности обработки больших объемов данных) в R все равно нет.
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 11:45 am (UTC)а при чем тут ваши трудности в портировании на питон к матлабу и математике? матлаб и математика на питон легче портируются? :)
что касается больших данных --- все для их обработки в R есть, уже даже в роликах ютюба для ленивых выкладывают.
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 12:02 pm (UTC)> же даже в роликах ютюба для ленивых выкладывают.
> матлаб и математика на питон легче портируются? :)
Да, аналитики не считают себя программистами и пишут, то что надо сделать, а не то как по их мнению это надо сделать.
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 12:22 pm (UTC)> же даже в роликах ютюба для ленивых выкладывают.
> матлаб и математика на питон легче портируются? :)
Да, аналитики не считают себя программистами и пишут, то что надо сделать, а не то как по их мнению это надо сделать.
---------------
Думаю что некая "ненависть" к R действительно обоснована. Язык действительно высоуровневый настолько что позволяет писать крайне компактные и легко развиваемые программы человеку для которого программирование является всего лишь инструментом, а не основной специальностью. Таким образом R сокращает рынок написания программного обеспечения "для ученых". Все сводится фактически к релизам прототипов написанных на нем. Естественно с использованием чего либо индустриального из языков программирования.
Что касается размещения объектов в памяти. Одинаковый массив занимает в памяти ровно столько же сколько и в матлабе и в математике. Да и динамически подгружать данные с диска возможно, и распределено хранить - обрабатывать сверхбольшие объемы возможно. Помимо свободных реализаций для сверх больших данных есть несколько проприетарных реализаций для R.
Оба этих изложенных момента упираются в квалификацию пользователя. Увы компания что ученому надо лишь "уметь тыкать мышкой" в офисных продуктах микрософт и различных родственных этим продуктам "ориджинах" сделала свое черное дело.... всем этим "ельфам" теперь нужны "морлоки" для отправления "естественных обязанностей". надеюсь они эльфов съедят в конце концов :)
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 01:10 pm (UTC)А можно узнать в каких это областях ? Потому что например в кв. химии и других вещах на кластерах и CUDA х до сих пор рулит С, и Fortran изредка с питоновскими API.
У биологов вроде тот же python и изредка, в последнее время, прости господи, Ruby. NLP - это Java и Python точнее NLTK.
Computer vision - С++ и матлаб почему то.
Зоопарк у "ученых" еще тот, и лично у меня лозунги типа "высокоуровневый язык для ученых" доверия не вызывают с тех пор как я прочитал про историю создания Фортрана. И с R то столкнулся впервые когда "ушел из науки", и "R сокращает рынок", это как то вообще громкое заявление.
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 02:14 pm (UTC)Электофизиологи какой то то хренью пользуются записывая свои сигналы... я бы сказал там даже секта целая :)
Биоинформатика явно "деградирует" :) раньше помню был perl во все поля :) так что руби вполне логичен.
В отличии от фортрана в основе появления R лежит свежепройденный курс по сицип от Засмана и Абельсона. Обычно люди кидаются свою систему типов писать, а тут взяли и клонировали S. Получилась "еще одна схема" которая по "заветам Столмана" служит высокоуровневым клеем для кода на си и фортране. Сейчас сишный код можно прямо мешать прямо в репле с R кодом.
Например R дает среду в которой сразу пишешь статью (литературное программирование), код практически декларативно описывает что надо делать и рисовать, синтаксис языка носителя нигде не загромождает смысл программы.
Если просто посмотреть, то все происходит в репле. Это делает незаметными расходы на вычисления (в 90% случаев). Возможность, не зарабатывать грыжу пиша полноценные программы, а просто модифицировать минимальными усилиями созданные объекты данных и кода.
За последний год R реализовал функционала больше чем SAS за все время существования :) R включен как составная часть во все крупные проприетарные системы анализа данных.
PS Машинное зрение это OpenCV + pcl + ROS, и писать там можно хоть на лиспе :)
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 04:21 pm (UTC)> R лежит свежепройденный курс по сицип от Засмана и Абельсона.
SICP я конечно смотрел и читал, но видимо я мало общался с R, да и то что я видел напоминало PHP style,
да и ... тоже сильно подкосило
неподготовленную к такой жести психику.
Код который я портировал чаще всего писался в спешке и т.д. Надо будет посмотреть повнимательнее. Спасибо.
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 04:40 pm (UTC)хотя итераторы написали давно, и в принципе должно хватать и так возможностей писать именно циклы.
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-24 07:29 pm (UTC)Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-24 07:35 pm (UTC)Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-04-04 01:15 am (UTC)легко развиваемые === одноразовые
Вы знаете, зачем вообще портируют Р код на Питон и проч? Для ПРИЛОЖЕНИЙ.
для ученых Р - рай. Использовал пакетик по специальности - хорошо (трибализм). Код нечитаемый - хорошо (элитизм).
Заявим, что все воспроизводимо, т.к. маркдаун/книтор/проч.дерьмо (карго культ).
Плюс пишут этот код далеко не всегда эти ученые.
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-04-04 07:01 am (UTC)никакой "нечитаемости" в нем нет.
большинство "пакетиков" посвящено именно заведению еще одного фрагмента "пазла". и этот фрагмент очень хорошо и удобно вписывается в уже имеющуюся систему и удачно её дополняет.
"легко развиваемые" и "нечитаемые" явно противоположные понятия. все рассуждения притянутые из программинга в случае R "мимо кассы". Язык используется в режиме интерактивной обработки данных. "программа" на нем обрабатывает в 99% процентах случаев конкретную выборку данных (и да, она не обязана быть корректной на другой выборке :) )
что касается "дерьма". побольше бы такого "дерьма" которое повышает производительность труда (конкретно написание статей и презентаций) на _два_ порядка.
PS ...а ученые и диссертации свои не всегда сами пишут :)
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-04-04 08:21 am (UTC)Золотые слова. Я выше это уже написал. Р - замечательный калькулятор, экономит кучу времени как учьоным, так и учёным. И нормальным людям тоже. Но это ужасный язык программирования, который поощряет нечитаемый одноразовый код. Нечитаемый. Одноразовый. 😝
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-04-04 09:36 am (UTC)Только надо помнить что: писать придется "как на схеме", с использованием объектной системы "как в лиспе", осознанно выбрав стиль кодирования в котором будет решена задача. При чем сам язык не накладывает _никаких_ ограничений.
Оба режима одинаково важны и нужны. Умение вовремя "переключиться" от одного подхода к другому и составляет "опытность" пользователя R.
2) Не читаемый и одноразовый код люди пишут на любом языке программирования :)
Код R вставленный (в 90% случаев 1-2 строчками) в текст статьи заменяет целые предложения (а то и страницы :) ) невразумительных объяснений "что авторы имели в виду". Если это "нечитаемость", то я за неё. :)
Re: A вы знаете язык лучше R справляющийся с задачами ана
Date: 2013-03-23 12:10 pm (UTC)no subject
Date: 2013-03-20 02:08 am (UTC)no subject
Date: 2013-03-20 06:25 am (UTC)гораздо проще ее отформатировать по человечески, удалить макросы, переименовать переменные и так разобраться
no subject
Date: 2013-03-20 10:55 am (UTC)#define P printf
#define R return
P - наПечатать (rint)
R - веРнуть (eturn)
no subject
Date: 2013-03-20 06:57 am (UTC)no subject
Date: 2013-03-20 10:08 am (UTC)Да, у Артура, голова работает явно по-другому. Но что-то притягательное в этом есть. Когда начинал писать на q (развитие k), всегда использовал длинные имена функций, которые всё объясняют. А потом стал всё больше писать однострочники - которые все перед глазами и в объяснениях не нуждаются.
Что меня в q/k поражает гораздо больше - это контекстозависимоть. $ - может быть как cast, так и if/else, так и case. Если второе и трерье ещё как-то связаны, то почему тот же самый символ используется для первого - для меня - полная загадка.
no subject
Date: 2013-03-20 12:41 pm (UTC)t - тип 0 или 1 (box?)
r - размерность,
d - размеры r-мерного прямоугольника
p - данные.
tr -- количество "элементарных" элементов
ga - конструктор
iota - создание 1-мерного массива [0, 1, ... n]
plus - сложение массивов одинаковой размерности и размеров
from - доступ к произвольному n-1 мерному ряду/элементу
box - запаковать в один элемент специального типа
cat -- конкатенировать и запаковать в линейный массив
size - количество рядов (первый размер)
rsh - тоже какое-то склеивание
sha - размеры
парсинг и выполнение лень было читать.
no subject
Date: 2013-03-20 12:49 pm (UTC)t - обозначает, что это на самом деле запакованный указатель на другой элемент
rsh - ключевая штука, потому что единственное, что позволяет увеличить размерность
в парсинге оригинально сделана поддержка той запаковки, что указана в t, но в целом ничего очень сложного нет
no subject
Date: 2013-03-20 01:02 pm (UTC)Создает массив с размерами из данных второго параметра (и размерностью = ширина второго параметра), и укладывате/обрезает данные первого параметра.
no subject
Date: 2013-03-23 07:00 pm (UTC)Язык, увы, не Тюринг-полный, так что не получится устроить конкурса на самую интересную программу на этом миниязыке.
Спасибо за публикацию. С удовольствием прочитал код. Совсем с листа не хватило терпения - в паре мест вставил переводы строки, догадки проверял запуском.
no subject
Date: 2013-03-20 04:41 pm (UTC)no subject
Date: 2013-03-20 04:42 pm (UTC)no subject
Date: 2013-03-20 06:58 pm (UTC)no subject
Date: 2013-03-20 07:57 pm (UTC)J - волшебный язык
Date: 2013-03-20 11:46 pm (UTC)Re: J - волшебный язык
Date: 2013-03-23 10:22 pm (UTC)С такого ммммм, оляля, языка программирования - просто недетская радость.
Re: J - волшебный язык
Date: 2013-03-24 06:17 am (UTC)