jk.c 164 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798991001011021031041051061071081091101111121131141151161171181191201211221231241251261271281291301311321331341351361371381391401411421431441451461471481491501511521531541551561571581591601611621631641651661671681691701711721731741751761771781791801811821831841851861871881891901911921931941951961971981992002012022032042052062072082092102112122132142152162172182192202212222232242252262272282292302312322332342352362372382392402412422432442452462472482492502512522532542552562572582592602612622632642652662672682692702712722732742752762772782792802812822832842852862872882892902912922932942952962972982993003013023033043053063073083093103113123133143153163173183193203213223233243253263273283293303313323333343353363373383393403413423433443453463473483493503513523533543553563573583593603613623633643653663673683693703713723733743753763773783793803813823833843853863873883893903913923933943953963973983994004014024034044054064074084094104114124134144154164174184194204214224234244254264274284294304314324334344354364374384394404414424434444454464474484494504514524534544554564574584594604614624634644654664674684694704714724734744754764774784794804814824834844854864874884894904914924934944954964974984995005015025035045055065075085095105115125135145155165175185195205215225235245255265275285295305315325335345355365375385395405415425435445455465475485495505515525535545555565575585595605615625635645655665675685695705715725735745755765775785795805815825835845855865875885895905915925935945955965975985996006016026036046056066076086096106116126136146156166176186196206216226236246256266276286296306316326336346356366376386396406416426436446456466476486496506516526536546556566576586596606616626636646656666676686696706716726736746756766776786796806816826836846856866876886896906916926936946956966976986997007017027037047057067077087097107117127137147157167177187197207217227237247257267277287297307317327337347357367377387397407417427437447457467477487497507517527537547557567577587597607617627637647657667677687697707717727737747757767777787797807817827837847857867877887897907917927937947957967977987998008018028038048058068078088098108118128138148158168178188198208218228238248258268278288298308318328338348358368378388398408418428438448458468478488498508518528538548558568578588598608618628638648658668678688698708718728738748758768778788798808818828838848858868878888898908918928938948958968978988999009019029039049059069079089099109119129139149159169179189199209219229239249259269279289299309319329339349359369379389399409419429439449459469479489499509519529539549559569579589599609619629639649659669679689699709719729739749759769779789799809819829839849859869879889899909919929939949959969979989991000100110021003100410051006100710081009101010111012101310141015101610171018101910201021102210231024102510261027102810291030103110321033103410351036103710381039104010411042104310441045104610471048104910501051105210531054105510561057105810591060106110621063106410651066106710681069107010711072107310741075107610771078107910801081108210831084108510861087108810891090109110921093109410951096109710981099110011011102110311041105110611071108110911101111111211131114111511161117111811191120112111221123112411251126112711281129113011311132113311341135113611371138113911401141114211431144114511461147114811491150115111521153115411551156115711581159116011611162116311641165116611671168116911701171117211731174117511761177117811791180118111821183118411851186118711881189119011911192119311941195119611971198119912001201120212031204120512061207120812091210121112121213121412151216121712181219122012211222122312241225122612271228122912301231123212331234123512361237123812391240124112421243124412451246124712481249125012511252125312541255125612571258125912601261126212631264126512661267126812691270127112721273127412751276127712781279128012811282128312841285128612871288128912901291129212931294129512961297129812991300130113021303130413051306130713081309131013111312131313141315131613171318131913201321132213231324132513261327132813291330133113321333133413351336133713381339134013411342134313441345134613471348134913501351135213531354135513561357135813591360136113621363136413651366136713681369137013711372137313741375137613771378137913801381138213831384138513861387138813891390139113921393139413951396139713981399140014011402140314041405140614071408140914101411141214131414141514161417141814191420142114221423142414251426142714281429143014311432143314341435143614371438143914401441144214431444144514461447144814491450145114521453145414551456145714581459146014611462146314641465146614671468146914701471147214731474147514761477147814791480148114821483148414851486148714881489149014911492149314941495149614971498149915001501150215031504150515061507150815091510151115121513151415151516151715181519152015211522152315241525152615271528152915301531153215331534153515361537153815391540154115421543154415451546154715481549155015511552155315541555155615571558155915601561156215631564156515661567156815691570157115721573157415751576157715781579158015811582158315841585158615871588158915901591159215931594159515961597159815991600160116021603160416051606160716081609161016111612161316141615161616171618161916201621162216231624162516261627162816291630163116321633163416351636163716381639164016411642164316441645164616471648164916501651165216531654165516561657165816591660166116621663166416651666166716681669167016711672167316741675167616771678167916801681168216831684168516861687168816891690169116921693169416951696169716981699170017011702170317041705170617071708170917101711171217131714171517161717171817191720172117221723172417251726172717281729173017311732173317341735173617371738173917401741174217431744174517461747174817491750175117521753175417551756175717581759176017611762176317641765176617671768176917701771177217731774177517761777177817791780178117821783178417851786178717881789179017911792179317941795179617971798179918001801180218031804180518061807180818091810181118121813181418151816181718181819182018211822182318241825182618271828182918301831183218331834183518361837183818391840184118421843184418451846184718481849185018511852185318541855185618571858185918601861186218631864186518661867186818691870187118721873187418751876187718781879188018811882188318841885188618871888188918901891189218931894189518961897189818991900190119021903190419051906190719081909191019111912191319141915191619171918191919201921192219231924192519261927192819291930193119321933193419351936193719381939194019411942194319441945194619471948194919501951195219531954195519561957195819591960196119621963196419651966196719681969197019711972197319741975197619771978197919801981198219831984198519861987198819891990199119921993199419951996199719981999200020012002200320042005200620072008200920102011201220132014201520162017201820192020202120222023202420252026202720282029203020312032203320342035203620372038203920402041204220432044204520462047204820492050205120522053205420552056205720582059206020612062206320642065206620672068206920702071207220732074207520762077207820792080208120822083208420852086208720882089209020912092209320942095209620972098209921002101210221032104210521062107210821092110211121122113211421152116211721182119212021212122212321242125212621272128212921302131213221332134213521362137213821392140214121422143214421452146214721482149215021512152215321542155215621572158215921602161216221632164216521662167216821692170217121722173217421752176217721782179218021812182218321842185218621872188218921902191219221932194219521962197219821992200220122022203220422052206220722082209221022112212221322142215221622172218221922202221222222232224222522262227222822292230223122322233223422352236223722382239224022412242224322442245224622472248224922502251225222532254225522562257225822592260226122622263226422652266226722682269227022712272227322742275227622772278227922802281228222832284228522862287228822892290229122922293229422952296229722982299230023012302230323042305230623072308230923102311231223132314231523162317231823192320232123222323232423252326232723282329233023312332233323342335233623372338233923402341234223432344234523462347234823492350235123522353235423552356235723582359236023612362236323642365236623672368236923702371237223732374237523762377237823792380238123822383238423852386238723882389239023912392239323942395239623972398239924002401240224032404240524062407240824092410241124122413241424152416241724182419242024212422242324242425242624272428242924302431243224332434243524362437243824392440244124422443244424452446244724482449245024512452245324542455245624572458245924602461246224632464246524662467246824692470247124722473247424752476247724782479248024812482248324842485248624872488248924902491249224932494249524962497249824992500250125022503250425052506250725082509251025112512251325142515251625172518251925202521252225232524252525262527252825292530253125322533253425352536253725382539254025412542254325442545254625472548254925502551255225532554255525562557255825592560256125622563256425652566256725682569257025712572257325742575257625772578257925802581258225832584258525862587258825892590259125922593259425952596259725982599260026012602260326042605260626072608260926102611261226132614261526162617261826192620262126222623262426252626262726282629263026312632263326342635263626372638263926402641264226432644264526462647264826492650265126522653265426552656265726582659266026612662266326642665266626672668266926702671267226732674267526762677267826792680268126822683268426852686268726882689269026912692269326942695269626972698269927002701270227032704270527062707270827092710271127122713271427152716271727182719272027212722272327242725272627272728272927302731273227332734273527362737273827392740274127422743274427452746274727482749275027512752275327542755275627572758275927602761276227632764276527662767276827692770277127722773277427752776277727782779278027812782278327842785278627872788278927902791279227932794279527962797279827992800280128022803280428052806280728082809281028112812281328142815281628172818281928202821282228232824282528262827282828292830283128322833283428352836283728382839284028412842284328442845284628472848284928502851285228532854285528562857285828592860286128622863286428652866286728682869287028712872287328742875287628772878287928802881288228832884288528862887288828892890289128922893289428952896289728982899290029012902290329042905290629072908290929102911291229132914291529162917291829192920292129222923292429252926292729282929293029312932293329342935293629372938293929402941294229432944294529462947294829492950295129522953295429552956295729582959296029612962296329642965296629672968296929702971297229732974297529762977297829792980298129822983298429852986298729882989299029912992299329942995299629972998299930003001300230033004300530063007300830093010301130123013301430153016301730183019302030213022302330243025302630273028302930303031303230333034303530363037303830393040304130423043304430453046304730483049305030513052305330543055305630573058305930603061306230633064306530663067306830693070307130723073307430753076307730783079308030813082308330843085308630873088308930903091309230933094309530963097309830993100310131023103310431053106310731083109311031113112311331143115311631173118311931203121312231233124312531263127312831293130313131323133313431353136313731383139314031413142314331443145314631473148314931503151315231533154315531563157315831593160316131623163316431653166316731683169317031713172317331743175317631773178317931803181318231833184318531863187318831893190319131923193319431953196319731983199320032013202320332043205320632073208320932103211321232133214321532163217321832193220322132223223322432253226322732283229323032313232323332343235323632373238323932403241324232433244324532463247324832493250325132523253325432553256325732583259326032613262326332643265326632673268326932703271327232733274327532763277327832793280328132823283328432853286328732883289329032913292329332943295329632973298329933003301330233033304330533063307330833093310331133123313331433153316331733183319332033213322332333243325332633273328332933303331333233333334333533363337333833393340334133423343334433453346334733483349335033513352335333543355335633573358335933603361336233633364336533663367336833693370337133723373337433753376337733783379338033813382338333843385338633873388338933903391339233933394339533963397339833993400340134023403340434053406340734083409341034113412341334143415341634173418341934203421342234233424342534263427342834293430343134323433343434353436343734383439344034413442344334443445344634473448344934503451345234533454345534563457345834593460346134623463346434653466346734683469347034713472347334743475347634773478347934803481348234833484348534863487348834893490349134923493349434953496349734983499350035013502350335043505350635073508350935103511351235133514351535163517351835193520352135223523352435253526352735283529353035313532353335343535353635373538353935403541354235433544354535463547354835493550355135523553355435553556355735583559356035613562356335643565356635673568356935703571357235733574357535763577357835793580358135823583358435853586358735883589359035913592359335943595359635973598359936003601360236033604360536063607360836093610361136123613361436153616361736183619362036213622362336243625362636273628362936303631363236333634363536363637363836393640364136423643364436453646364736483649365036513652365336543655365636573658365936603661366236633664366536663667366836693670367136723673367436753676367736783679368036813682368336843685368636873688368936903691369236933694369536963697369836993700370137023703370437053706370737083709371037113712371337143715371637173718371937203721372237233724372537263727372837293730373137323733373437353736373737383739374037413742374337443745374637473748374937503751375237533754375537563757375837593760376137623763376437653766376737683769377037713772377337743775377637773778377937803781378237833784378537863787378837893790379137923793379437953796379737983799380038013802380338043805380638073808380938103811381238133814381538163817381838193820382138223823382438253826382738283829383038313832383338343835383638373838383938403841384238433844384538463847384838493850385138523853385438553856385738583859386038613862386338643865386638673868386938703871387238733874387538763877387838793880388138823883388438853886388738883889389038913892389338943895389638973898389939003901390239033904390539063907390839093910391139123913391439153916391739183919392039213922392339243925392639273928392939303931393239333934393539363937393839393940394139423943394439453946394739483949395039513952395339543955395639573958395939603961396239633964396539663967396839693970397139723973397439753976397739783979398039813982398339843985398639873988398939903991399239933994399539963997399839994000400140024003400440054006400740084009401040114012401340144015401640174018401940204021402240234024402540264027402840294030403140324033403440354036403740384039404040414042404340444045404640474048404940504051405240534054405540564057405840594060406140624063406440654066406740684069407040714072407340744075407640774078407940804081408240834084408540864087408840894090409140924093409440954096409740984099410041014102410341044105410641074108410941104111411241134114411541164117411841194120412141224123412441254126412741284129413041314132413341344135413641374138413941404141414241434144414541464147414841494150415141524153415441554156415741584159416041614162416341644165416641674168416941704171417241734174417541764177417841794180418141824183418441854186418741884189419041914192419341944195419641974198419942004201420242034204420542064207420842094210421142124213421442154216421742184219422042214222422342244225422642274228422942304231423242334234423542364237423842394240424142424243424442454246424742484249425042514252425342544255425642574258425942604261426242634264426542664267426842694270427142724273427442754276427742784279428042814282428342844285428642874288428942904291429242934294429542964297429842994300430143024303430443054306430743084309431043114312431343144315431643174318431943204321432243234324432543264327432843294330433143324333433443354336433743384339434043414342434343444345434643474348434943504351435243534354435543564357435843594360436143624363436443654366436743684369437043714372437343744375437643774378437943804381438243834384438543864387438843894390439143924393439443954396439743984399440044014402440344044405440644074408440944104411441244134414441544164417441844194420442144224423442444254426442744284429443044314432443344344435443644374438443944404441444244434444444544464447444844494450445144524453445444554456445744584459446044614462446344644465446644674468446944704471447244734474447544764477447844794480448144824483448444854486448744884489449044914492449344944495449644974498449945004501450245034504450545064507450845094510451145124513451445154516451745184519452045214522452345244525452645274528452945304531453245334534453545364537453845394540454145424543454445454546454745484549455045514552455345544555455645574558455945604561456245634564456545664567456845694570457145724573457445754576457745784579458045814582458345844585458645874588458945904591459245934594459545964597459845994600460146024603460446054606460746084609461046114612461346144615461646174618461946204621462246234624462546264627462846294630463146324633463446354636463746384639464046414642464346444645464646474648464946504651465246534654465546564657465846594660466146624663466446654666466746684669467046714672467346744675467646774678467946804681468246834684468546864687468846894690469146924693469446954696469746984699470047014702470347044705470647074708470947104711471247134714471547164717471847194720472147224723472447254726472747284729473047314732473347344735473647374738473947404741474247434744474547464747474847494750475147524753475447554756475747584759476047614762476347644765476647674768476947704771477247734774477547764777477847794780478147824783478447854786478747884789479047914792479347944795479647974798479948004801480248034804480548064807480848094810481148124813481448154816481748184819482048214822482348244825482648274828482948304831483248334834483548364837483848394840484148424843484448454846484748484849485048514852485348544855485648574858485948604861486248634864486548664867486848694870487148724873487448754876487748784879488048814882488348844885488648874888488948904891489248934894489548964897489848994900490149024903490449054906490749084909491049114912491349144915491649174918491949204921492249234924492549264927492849294930493149324933493449354936493749384939494049414942494349444945494649474948494949504951495249534954495549564957495849594960496149624963496449654966496749684969497049714972497349744975497649774978497949804981498249834984498549864987498849894990499149924993499449954996499749984999500050015002500350045005500650075008500950105011501250135014501550165017501850195020502150225023502450255026502750285029503050315032503350345035503650375038503950405041504250435044504550465047504850495050505150525053505450555056505750585059506050615062506350645065506650675068506950705071507250735074507550765077507850795080508150825083508450855086508750885089509050915092509350945095509650975098509951005101510251035104510551065107510851095110511151125113511451155116511751185119512051215122512351245125512651275128512951305131513251335134513551365137513851395140514151425143514451455146514751485149515051515152515351545155515651575158515951605161516251635164516551665167516851695170517151725173517451755176517751785179518051815182518351845185518651875188518951905191519251935194519551965197519851995200520152025203520452055206520752085209521052115212521352145215521652175218521952205221522252235224522552265227522852295230523152325233523452355236523752385239524052415242524352445245524652475248524952505251525252535254525552565257525852595260526152625263526452655266526752685269527052715272527352745275527652775278527952805281528252835284528552865287528852895290529152925293529452955296529752985299530053015302530353045305530653075308530953105311531253135314531553165317531853195320532153225323532453255326532753285329533053315332533353345335533653375338533953405341534253435344534553465347534853495350535153525353535453555356535753585359536053615362536353645365536653675368536953705371537253735374537553765377537853795380538153825383538453855386538753885389539053915392539353945395539653975398539954005401540254035404540554065407540854095410541154125413541454155416541754185419542054215422542354245425542654275428542954305431543254335434543554365437543854395440544154425443544454455446544754485449545054515452545354545455545654575458545954605461546254635464546554665467546854695470547154725473547454755476547754785479548054815482548354845485548654875488548954905491549254935494549554965497549854995500550155025503550455055506550755085509551055115512551355145515551655175518551955205521552255235524552555265527552855295530553155325533553455355536553755385539554055415542554355445545554655475548554955505551555255535554555555565557555855595560556155625563556455655566556755685569557055715572557355745575557655775578557955805581558255835584558555865587558855895590559155925593559455955596559755985599560056015602560356045605560656075608560956105611561256135614561556165617561856195620562156225623562456255626562756285629563056315632563356345635563656375638563956405641564256435644564556465647564856495650565156525653565456555656565756585659566056615662566356645665566656675668566956705671567256735674567556765677567856795680568156825683568456855686568756885689569056915692569356945695569656975698569957005701570257035704570557065707570857095710571157125713571457155716571757185719572057215722572357245725572657275728572957305731573257335734573557365737573857395740574157425743574457455746574757485749575057515752575357545755575657575758575957605761576257635764576557665767576857695770577157725773577457755776577757785779578057815782578357845785578657875788578957905791579257935794579557965797579857995800580158025803580458055806580758085809581058115812581358145815581658175818581958205821582258235824582558265827582858295830583158325833583458355836583758385839584058415842584358445845584658475848584958505851585258535854585558565857585858595860586158625863586458655866586758685869587058715872587358745875587658775878587958805881588258835884588558865887588858895890589158925893589458955896589758985899590059015902590359045905590659075908590959105911591259135914591559165917591859195920592159225923592459255926592759285929593059315932593359345935593659375938593959405941594259435944594559465947594859495950595159525953595459555956595759585959596059615962596359645965596659675968596959705971597259735974597559765977597859795980598159825983598459855986598759885989599059915992599359945995599659975998599960006001600260036004600560066007600860096010601160126013601460156016601760186019602060216022602360246025602660276028602960306031603260336034603560366037603860396040604160426043604460456046604760486049605060516052605360546055605660576058605960606061606260636064606560666067606860696070607160726073607460756076607760786079608060816082608360846085608660876088608960906091609260936094609560966097609860996100610161026103610461056106610761086109611061116112611361146115611661176118611961206121612261236124612561266127612861296130613161326133613461356136613761386139614061416142614361446145614661476148614961506151615261536154615561566157615861596160616161626163616461656166616761686169617061716172617361746175617661776178617961806181618261836184618561866187618861896190619161926193619461956196619761986199620062016202620362046205620662076208620962106211621262136214621562166217621862196220622162226223622462256226622762286229623062316232623362346235623662376238623962406241624262436244624562466247624862496250625162526253625462556256625762586259626062616262626362646265626662676268626962706271627262736274627562766277627862796280628162826283628462856286628762886289629062916292629362946295629662976298629963006301630263036304630563066307630863096310631163126313631463156316631763186319632063216322632363246325632663276328632963306331633263336334633563366337633863396340634163426343634463456346634763486349635063516352635363546355635663576358635963606361636263636364636563666367636863696370637163726373637463756376637763786379638063816382638363846385638663876388638963906391639263936394639563966397639863996400640164026403640464056406640764086409641064116412641364146415641664176418641964206421642264236424642564266427642864296430643164326433643464356436643764386439644064416442644364446445644664476448644964506451645264536454645564566457645864596460646164626463646464656466646764686469647064716472647364746475647664776478647964806481648264836484648564866487648864896490649164926493649464956496649764986499650065016502650365046505650665076508650965106511651265136514651565166517651865196520652165226523652465256526652765286529653065316532653365346535653665376538653965406541654265436544654565466547654865496550655165526553655465556556655765586559656065616562656365646565656665676568656965706571657265736574657565766577657865796580658165826583658465856586658765886589659065916592659365946595659665976598659966006601660266036604660566066607660866096610661166126613661466156616661766186619662066216622662366246625662666276628662966306631663266336634663566366637663866396640664166426643664466456646664766486649665066516652665366546655665666576658665966606661666266636664666566666667666866696670667166726673667466756676667766786679668066816682668366846685668666876688668966906691669266936694669566966697669866996700670167026703670467056706670767086709671067116712671367146715671667176718671967206721672267236724672567266727672867296730673167326733673467356736673767386739674067416742674367446745674667476748674967506751675267536754675567566757675867596760676167626763676467656766676767686769677067716772677367746775677667776778677967806781678267836784678567866787678867896790679167926793679467956796679767986799680068016802680368046805680668076808680968106811681268136814681568166817681868196820682168226823682468256826682768286829683068316832683368346835683668376838683968406841684268436844684568466847684868496850685168526853685468556856685768586859686068616862686368646865686668676868686968706871687268736874687568766877687868796880688168826883688468856886688768886889689068916892689368946895689668976898689969006901690269036904690569066907690869096910691169126913691469156916691769186919692069216922692369246925692669276928692969306931693269336934693569366937693869396940694169426943694469456946694769486949695069516952695369546955695669576958695969606961696269636964696569666967696869696970697169726973697469756976697769786979698069816982698369846985698669876988698969906991699269936994699569966997699869997000700170027003700470057006700770087009701070117012701370147015701670177018701970207021702270237024702570267027702870297030703170327033703470357036703770387039704070417042704370447045704670477048704970507051705270537054705570567057705870597060706170627063706470657066706770687069707070717072707370747075707670777078707970807081708270837084708570867087708870897090709170927093709470957096709770987099710071017102710371047105710671077108710971107111711271137114711571167117711871197120712171227123712471257126712771287129713071317132713371347135713671377138713971407141714271437144714571467147714871497150715171527153715471557156715771587159716071617162716371647165716671677168716971707171717271737174717571767177717871797180718171827183718471857186718771887189719071917192719371947195719671977198719972007201720272037204720572067207720872097210721172127213721472157216721772187219722072217222722372247225722672277228722972307231723272337234723572367237723872397240724172427243724472457246724772487249725072517252725372547255725672577258725972607261726272637264726572667267726872697270727172727273727472757276727772787279728072817282728372847285728672877288728972907291729272937294729572967297729872997300730173027303730473057306730773087309731073117312731373147315731673177318731973207321732273237324732573267327732873297330733173327333733473357336733773387339734073417342734373447345734673477348734973507351735273537354735573567357735873597360736173627363736473657366736773687369737073717372737373747375737673777378737973807381738273837384738573867387738873897390739173927393739473957396739773987399740074017402740374047405740674077408740974107411741274137414741574167417741874197420742174227423742474257426742774287429743074317432743374347435743674377438743974407441744274437444744574467447744874497450745174527453745474557456745774587459746074617462746374647465746674677468746974707471
  1. #include "version.h"
  2. #include <ctype.h>
  3. #include <dlfcn.h>
  4. #include <errno.h>
  5. #include <ffi.h>
  6. #include <gc.h>
  7. #include <limits.h>
  8. #include <math.h>
  9. #include <pcre.h>
  10. #include <setjmp.h>
  11. #include <stdbool.h>
  12. #include <stddef.h>
  13. #include <stdint.h>
  14. #include <stdio.h>
  15. #include <stdlib.h>
  16. #include <string.h>
  17. #include <sys/wait.h>
  18. #include <time.h>
  19. #include <unistd.h>
  20. jmp_buf interactive_checkpoint;
  21. bool is_interactive;
  22. size_t max_rec_depth = 1000;
  23. size_t rec_depth = 0;
  24. void *malloc_checked(size_t size) {
  25. void *p;
  26. if (!(p = GC_MALLOC(size)))
  27. abort();
  28. memset(p, 0, size);
  29. return p;
  30. }
  31. void *malloc_checked_atomic(size_t size) {
  32. void *p;
  33. if (!(p = GC_malloc_atomic(size)))
  34. abort();
  35. memset(p, 0, size);
  36. return p;
  37. }
  38. void *malloc_checked_uncollectable(size_t size) {
  39. void *p;
  40. if (!(p = GC_malloc_uncollectable(size)))
  41. abort();
  42. memset(p, 0, size);
  43. return p;
  44. }
  45. void *realloc_checked(void *p, size_t size) {
  46. if (!(p = GC_REALLOC(p, size)))
  47. abort();
  48. return p;
  49. }
  50. char *strdup_checked(char *s) {
  51. char *p = GC_strdup(s);
  52. if (!p)
  53. abort();
  54. return p;
  55. }
  56. typedef struct {
  57. void **data;
  58. size_t length;
  59. } list_t;
  60. list_t *list_new(void) {
  61. list_t *list = malloc_checked(sizeof(list_t));
  62. list->data = NULL;
  63. list->length = 0;
  64. return list;
  65. }
  66. list_t *list_newk(size_t k) {
  67. list_t *list = malloc_checked(sizeof(list_t));
  68. list->data = malloc_checked(k * sizeof(void *));
  69. list->length = k;
  70. return list;
  71. }
  72. list_t *list_copy(list_t *l) {
  73. list_t *list = list_newk(l->length);
  74. for (size_t i = 0; i < l->length; i++)
  75. list->data[i] = l->data[i];
  76. return list;
  77. }
  78. void list_push(list_t *l, void *v) {
  79. size_t i = l->length++;
  80. l->data = realloc_checked(l->data, l->length * sizeof(void *));
  81. l->data[i] = v;
  82. }
  83. void *list_pop(list_t *l) {
  84. if (!l->data)
  85. return NULL;
  86. size_t i = --l->length;
  87. void *v = l->data[i];
  88. l->data[i] = NULL;
  89. if (!l->length) {
  90. GC_FREE(l->data);
  91. l->data = NULL;
  92. } else
  93. l->data = realloc_checked(l->data, l->length * sizeof(void *));
  94. return v;
  95. }
  96. void *list_index(list_t *l, ssize_t index) {
  97. if (!l->data)
  98. return NULL;
  99. if (index < 0)
  100. index += ((ssize_t)l->length);
  101. if (index < 0 || index >= l->length)
  102. return NULL;
  103. return l->data[index];
  104. }
  105. void list_set(list_t *l, ssize_t index, void *v) {
  106. if (!l->data)
  107. return;
  108. if (index < 0)
  109. index += ((ssize_t)l->length);
  110. if (index < 0 || index >= l->length)
  111. return;
  112. l->data[index] = v;
  113. }
  114. typedef struct {
  115. char *str;
  116. size_t size;
  117. } buffer_t;
  118. buffer_t *buffer_new(void) {
  119. buffer_t *buf = malloc_checked(sizeof(buffer_t));
  120. buf->str = NULL;
  121. buf->size = 0;
  122. return buf;
  123. }
  124. void buffer_append(buffer_t *buf, char c) {
  125. buf->size++;
  126. void *p = malloc_checked_atomic(sizeof(char) * buf->size);
  127. if (buf->str) {
  128. memcpy(p, buf->str, buf->size - 1);
  129. GC_FREE(buf->str);
  130. }
  131. buf->str = p;
  132. buf->str[buf->size - 1] = c;
  133. }
  134. char *buffer_read(buffer_t *buf) {
  135. if (buf->size == 0 || buf->str[buf->size - 1])
  136. buffer_append(buf, 0);
  137. char *str = buf->str;
  138. GC_FREE(buf);
  139. return str;
  140. }
  141. void buffer_append_str(buffer_t *buf, char *s) {
  142. for (size_t i = 0; i < strlen(s); i++)
  143. buffer_append(buf, s[i]);
  144. }
  145. typedef struct {
  146. enum token_tag_t {
  147. T_PUNCT,
  148. T_LPAR,
  149. T_RPAR,
  150. T_NAME,
  151. T_NUMBER,
  152. T_BNUMBER,
  153. T_QUOTE
  154. } tag;
  155. char *text;
  156. } token_t;
  157. typedef struct {
  158. char *source;
  159. size_t len;
  160. size_t pos;
  161. list_t *tokens;
  162. } lexer_t;
  163. lexer_t *lexer_new(void) {
  164. lexer_t *lexer = malloc_checked(sizeof(lexer_t));
  165. return lexer;
  166. }
  167. char lexer_lookahead(lexer_t *lexer, size_t offset) {
  168. size_t pos = lexer->pos + offset;
  169. if (pos >= lexer->len)
  170. return 0;
  171. return lexer->source[pos];
  172. }
  173. char lexer_eat(lexer_t *lexer) {
  174. if (lexer->pos >= lexer->len)
  175. return 0;
  176. return lexer->source[lexer->pos++];
  177. }
  178. void lexer_push_token(lexer_t *lexer, enum token_tag_t tag, char *text) {
  179. token_t *token = malloc_checked(sizeof(token_t));
  180. token->tag = tag;
  181. token->text = text;
  182. list_push(lexer->tokens, token);
  183. }
  184. list_t *guards;
  185. typedef struct {
  186. jmp_buf lb;
  187. size_t rec_depth;
  188. } guard_t;
  189. guard_t *guard() {
  190. guard_t *g = malloc_checked_atomic(sizeof(guard_t));
  191. g->rec_depth = rec_depth;
  192. list_push(guards, g);
  193. return g;
  194. }
  195. guard_t *guarding() { return list_index(guards, -1); }
  196. void unguard() { GC_FREE(list_pop(guards)); }
  197. void fatal(char *s) {
  198. guard_t *g = guarding();
  199. if (g) {
  200. rec_depth = g->rec_depth;
  201. longjmp(g->lb, 1);
  202. }
  203. fprintf(stderr, "|%s error\n", s);
  204. if (is_interactive) {
  205. rec_depth = 0;
  206. longjmp(interactive_checkpoint, 1 );
  207. }
  208. exit(1);
  209. }
  210. void lexer_error(lexer_t *lexer, char *s) { fatal(s); }
  211. void lexer_lex_number(lexer_t *lexer, bool is_negative) {
  212. buffer_t *buf = buffer_new();
  213. if (is_negative)
  214. buffer_append(buf, '-');
  215. if (lexer_lookahead(lexer, 0) == '.') {
  216. buffer_append(buf, lexer_eat(lexer));
  217. if (!(isdigit(lexer_lookahead(lexer, 0))))
  218. lexer_error(lexer, "trailing-dot");
  219. }
  220. do {
  221. buffer_append(buf, lexer_eat(lexer));
  222. if (lexer_lookahead(lexer, 0) == '`' && isdigit(lexer_lookahead(lexer, 1)))
  223. lexer_eat(lexer);
  224. } while (isdigit(lexer_lookahead(lexer, 0)));
  225. if (lexer_lookahead(lexer, 0) == '.') {
  226. buffer_append(buf, lexer_eat(lexer));
  227. if (!(isdigit(lexer_lookahead(lexer, 0))))
  228. lexer_error(lexer, "trailing-dot");
  229. do {
  230. buffer_append(buf, lexer_eat(lexer));
  231. } while (isdigit(lexer_lookahead(lexer, 0)));
  232. }
  233. lexer_push_token(lexer, T_NUMBER, buffer_read(buf));
  234. }
  235. void lexer_lex(lexer_t *lexer, char *s) {
  236. lexer->source = s;
  237. lexer->len = strlen(s);
  238. lexer->pos = 0;
  239. lexer->tokens = list_new();
  240. while (lexer->pos < lexer->len) {
  241. char c = lexer_lookahead(lexer, 0);
  242. if (c == '/' && !lexer->tokens->data)
  243. break;
  244. if (isspace(c)) {
  245. lexer_eat(lexer);
  246. if (lexer_lookahead(lexer, 0) == '/')
  247. break;
  248. } else if (c == '0' && (lexer_lookahead(lexer, 1) == 'x' ||
  249. lexer_lookahead(lexer, 1) == 'b' ||
  250. lexer_lookahead(lexer, 1) == 'o')) {
  251. lexer_eat(lexer);
  252. buffer_t *buf = buffer_new();
  253. char b = lexer_eat(lexer);
  254. buffer_append(buf, b);
  255. const char *base = b == 'x' ? "0123456789abcdefABCDEF"
  256. : b == 'b' ? "01"
  257. : "01234567";
  258. while (strchr(base, lexer_lookahead(lexer, 0)) != NULL)
  259. buffer_append(buf, lexer_eat(lexer));
  260. lexer_push_token(lexer, T_BNUMBER, buffer_read(buf));
  261. } else if (isdigit(c) || c == '.') {
  262. lexer_lex_number(lexer, false);
  263. } else if (isalpha(c)) {
  264. buffer_t *buf = buffer_new();
  265. do {
  266. buffer_append(buf, lexer_eat(lexer));
  267. } while (isalpha(lexer_lookahead(lexer, 0)));
  268. if (buf->size == 1 && lexer_lookahead(lexer, 0) == '.') {
  269. buffer_append(buf, lexer_eat(lexer));
  270. lexer_push_token(lexer, T_PUNCT, buffer_read(buf));
  271. } else
  272. lexer_push_token(lexer, T_NAME, buffer_read(buf));
  273. } else if (c == '(' || c == ')') {
  274. lexer_eat(lexer);
  275. lexer_push_token(lexer, c == '(' ? T_LPAR : T_RPAR, NULL);
  276. } else if (c == '\'') {
  277. buffer_t *buf = buffer_new();
  278. lexer_eat(lexer);
  279. for (;;) {
  280. if (lexer->pos >= lexer->len)
  281. lexer_error(lexer, "unmatched-quote");
  282. if (lexer_lookahead(lexer, 0) == '\'') {
  283. if (lexer_lookahead(lexer, 1) == '\'') {
  284. buffer_append(buf, lexer_eat(lexer));
  285. lexer_eat(lexer);
  286. continue;
  287. }
  288. lexer_eat(lexer);
  289. break;
  290. }
  291. buffer_append(buf, lexer_eat(lexer));
  292. }
  293. lexer_push_token(lexer, T_QUOTE, buffer_read(buf));
  294. } else if (ispunct(c)) {
  295. char buf[3];
  296. buf[0] = lexer_eat(lexer);
  297. buf[1] = 0;
  298. if (lexer_lookahead(lexer, 0) == '.' ||
  299. lexer_lookahead(lexer, 0) == ':') {
  300. buf[1] = lexer_eat(lexer);
  301. buf[2] = 0;
  302. }
  303. if (strcmp(buf, "-") == 0 && isdigit(lexer_lookahead(lexer, 0))) {
  304. lexer_lex_number(lexer, true);
  305. continue;
  306. }
  307. lexer_push_token(lexer, T_PUNCT, strdup_checked(buf));
  308. } else
  309. lexer_error(lexer, "lex");
  310. }
  311. }
  312. typedef struct _table_t table_t;
  313. typedef struct _table_entry_t table_entry_t;
  314. struct _table_entry_t {
  315. char *key;
  316. void *value;
  317. bool is_deleted;
  318. };
  319. struct _table_t {
  320. table_entry_t *entries;
  321. size_t used;
  322. size_t capacity;
  323. };
  324. #define TABLE_MIN_SIZE 32
  325. table_t *table_new(void) {
  326. table_t *table = malloc_checked(sizeof(table_t));
  327. table->used = 0;
  328. table->capacity = TABLE_MIN_SIZE;
  329. table->entries = malloc_checked(table->capacity * sizeof(table_entry_t));
  330. return table;
  331. }
  332. size_t table_length(table_t *table) { return table->used; }
  333. bool table_empty(table_t *table) { return table->used == 0; }
  334. static uint64_t MM86128(void *key, const int len, uint32_t seed) {
  335. #define ROTL32(x, r) ((x << r) | (x >> (32 - r)))
  336. #define FMIX32(h) \
  337. h ^= h >> 16; \
  338. h *= 0x85ebca6b; \
  339. h ^= h >> 13; \
  340. h *= 0xc2b2ae35; \
  341. h ^= h >> 16;
  342. const uint8_t *data = (const uint8_t *)key;
  343. const int nblocks = len / 16;
  344. uint32_t h1 = seed;
  345. uint32_t h2 = seed;
  346. uint32_t h3 = seed;
  347. uint32_t h4 = seed;
  348. uint32_t c1 = 0x239b961b;
  349. uint32_t c2 = 0xab0e9789;
  350. uint32_t c3 = 0x38b34ae5;
  351. uint32_t c4 = 0xa1e38b93;
  352. const uint32_t *blocks = (const uint32_t *)(data + nblocks * 16);
  353. for (int i = -nblocks; i; i++) {
  354. uint32_t k1 = blocks[i * 4 + 0];
  355. uint32_t k2 = blocks[i * 4 + 1];
  356. uint32_t k3 = blocks[i * 4 + 2];
  357. uint32_t k4 = blocks[i * 4 + 3];
  358. k1 *= c1;
  359. k1 = ROTL32(k1, 15);
  360. k1 *= c2;
  361. h1 ^= k1;
  362. h1 = ROTL32(h1, 19);
  363. h1 += h2;
  364. h1 = h1 * 5 + 0x561ccd1b;
  365. k2 *= c2;
  366. k2 = ROTL32(k2, 16);
  367. k2 *= c3;
  368. h2 ^= k2;
  369. h2 = ROTL32(h2, 17);
  370. h2 += h3;
  371. h2 = h2 * 5 + 0x0bcaa747;
  372. k3 *= c3;
  373. k3 = ROTL32(k3, 17);
  374. k3 *= c4;
  375. h3 ^= k3;
  376. h3 = ROTL32(h3, 15);
  377. h3 += h4;
  378. h3 = h3 * 5 + 0x96cd1c35;
  379. k4 *= c4;
  380. k4 = ROTL32(k4, 18);
  381. k4 *= c1;
  382. h4 ^= k4;
  383. h4 = ROTL32(h4, 13);
  384. h4 += h1;
  385. h4 = h4 * 5 + 0x32ac3b17;
  386. }
  387. const uint8_t *tail = (const uint8_t *)(data + nblocks * 16);
  388. uint32_t k1 = 0;
  389. uint32_t k2 = 0;
  390. uint32_t k3 = 0;
  391. uint32_t k4 = 0;
  392. switch (len & 15) {
  393. case 15:
  394. k4 ^= tail[14] << 16;
  395. case 14:
  396. k4 ^= tail[13] << 8;
  397. case 13:
  398. k4 ^= tail[12] << 0;
  399. k4 *= c4;
  400. k4 = ROTL32(k4, 18);
  401. k4 *= c1;
  402. h4 ^= k4;
  403. case 12:
  404. k3 ^= tail[11] << 24;
  405. case 11:
  406. k3 ^= tail[10] << 16;
  407. case 10:
  408. k3 ^= tail[9] << 8;
  409. case 9:
  410. k3 ^= tail[8] << 0;
  411. k3 *= c3;
  412. k3 = ROTL32(k3, 17);
  413. k3 *= c4;
  414. h3 ^= k3;
  415. case 8:
  416. k2 ^= tail[7] << 24;
  417. case 7:
  418. k2 ^= tail[6] << 16;
  419. case 6:
  420. k2 ^= tail[5] << 8;
  421. case 5:
  422. k2 ^= tail[4] << 0;
  423. k2 *= c2;
  424. k2 = ROTL32(k2, 16);
  425. k2 *= c3;
  426. h2 ^= k2;
  427. case 4:
  428. k1 ^= tail[3] << 24;
  429. case 3:
  430. k1 ^= tail[2] << 16;
  431. case 2:
  432. k1 ^= tail[1] << 8;
  433. case 1:
  434. k1 ^= tail[0] << 0;
  435. k1 *= c1;
  436. k1 = ROTL32(k1, 15);
  437. k1 *= c2;
  438. h1 ^= k1;
  439. }
  440. h1 ^= len;
  441. h2 ^= len;
  442. h3 ^= len;
  443. h4 ^= len;
  444. h1 += h2;
  445. h1 += h3;
  446. h1 += h4;
  447. h2 += h1;
  448. h3 += h1;
  449. h4 += h1;
  450. FMIX32(h1);
  451. FMIX32(h2);
  452. FMIX32(h3);
  453. FMIX32(h4);
  454. h1 += h2;
  455. h1 += h3;
  456. h1 += h4;
  457. h2 += h1;
  458. h3 += h1;
  459. h4 += h1;
  460. return (((uint64_t)h2) << 32) | h1;
  461. }
  462. static uint32_t HASH_SEED = 0;
  463. void *table_get(table_t *table, char *key) {
  464. if (table_empty(table))
  465. return NULL;
  466. uint64_t hash = MM86128(key, strlen(key), HASH_SEED);
  467. size_t index = hash % table->capacity;
  468. size_t i = index;
  469. while (table->entries[i].key) {
  470. if (!table->entries[i].is_deleted &&
  471. strcmp(table->entries[i].key, key) == 0)
  472. return table->entries[i].value;
  473. i++;
  474. if (i >= table->capacity)
  475. i = 0;
  476. if (i == index)
  477. break;
  478. }
  479. return NULL;
  480. }
  481. bool table_has(table_t *table, char *key) {
  482. if (table_empty(table))
  483. return false;
  484. uint64_t hash = MM86128(key, strlen(key), HASH_SEED);
  485. size_t index = hash % table->capacity;
  486. size_t i = index;
  487. while (table->entries[i].key) {
  488. if (!table->entries[i].is_deleted &&
  489. strcmp(table->entries[i].key, key) == 0)
  490. return true;
  491. i++;
  492. if (i >= table->capacity)
  493. i = 0;
  494. if (i == index)
  495. break;
  496. }
  497. return false;
  498. }
  499. static void table_entry_set(table_entry_t *entries, char *key, void *value,
  500. size_t capacity, size_t *used) {
  501. uint64_t hash = MM86128(key, strlen(key), HASH_SEED);
  502. size_t index = hash % capacity;
  503. size_t i = index;
  504. while (entries[i].key) {
  505. if (strcmp(entries[i].key, key) == 0) {
  506. entries[i].value = value;
  507. if (entries[i].is_deleted) {
  508. if (used)
  509. (*used)++;
  510. entries[i].is_deleted = false;
  511. }
  512. return;
  513. } else if (entries[i].is_deleted)
  514. break;
  515. i++;
  516. if (i >= capacity)
  517. i = 0;
  518. if (i == index)
  519. break;
  520. }
  521. if (used)
  522. (*used)++;
  523. entries[i].key = key;
  524. entries[i].value = value;
  525. entries[i].is_deleted = false;
  526. }
  527. table_t *table_set(table_t *table, char *key, void *value) {
  528. if (table->used >= table->capacity) {
  529. size_t capacity = table->capacity + TABLE_MIN_SIZE;
  530. table_entry_t *entries = malloc_checked(capacity * sizeof(table_entry_t));
  531. for (size_t i = 0; i < table->capacity; i++) {
  532. table_entry_t entry = table->entries[i];
  533. if (entry.key && !entry.is_deleted)
  534. table_entry_set(entries, entry.key, entry.value, capacity, NULL);
  535. }
  536. GC_FREE(table->entries);
  537. table->entries = entries;
  538. table->capacity = capacity;
  539. }
  540. table_entry_set(table->entries, key, value, table->capacity, &table->used);
  541. return table;
  542. }
  543. bool table_delete(table_t *table, char *key) {
  544. uint64_t hash = MM86128(key, strlen(key), HASH_SEED);
  545. size_t index = hash % table->capacity;
  546. size_t i = index;
  547. while (table->entries[i].key) {
  548. if (!table->entries[i].is_deleted &&
  549. strcmp(table->entries[i].key, key) == 0) {
  550. table->entries[i].value = NULL;
  551. table->entries[i].is_deleted = true;
  552. table->used--;
  553. if (table->capacity > TABLE_MIN_SIZE &&
  554. table->used <= table->capacity - TABLE_MIN_SIZE) {
  555. size_t capacity = table->capacity - TABLE_MIN_SIZE;
  556. table_entry_t *entries =
  557. malloc_checked(capacity * sizeof(table_entry_t));
  558. for (size_t i = 0; i < table->capacity; i++) {
  559. table_entry_t entry = table->entries[i];
  560. if (entry.key && !entry.is_deleted)
  561. table_entry_set(entries, entry.key, entry.value, capacity, NULL);
  562. }
  563. GC_FREE(table->entries);
  564. table->entries = entries;
  565. table->capacity = capacity;
  566. }
  567. return true;
  568. }
  569. i++;
  570. if (i >= table->capacity)
  571. i = 0;
  572. if (i == index)
  573. break;
  574. }
  575. return false;
  576. }
  577. typedef struct _value_t value_t;
  578. typedef struct _interpreter_t interpreter_t;
  579. typedef struct _verb_t verb_t;
  580. struct _interpreter_t {
  581. table_t *env;
  582. list_t *args;
  583. list_t *selfrefs;
  584. value_t *nil;
  585. value_t *udf;
  586. value_t *unit;
  587. verb_t *at;
  588. bool bn;
  589. };
  590. struct _verb_t {
  591. char *name;
  592. unsigned int rank[3];
  593. list_t *bonds;
  594. bool mark;
  595. bool is_fun;
  596. value_t *(*monad)(interpreter_t *, verb_t *, value_t *);
  597. value_t *(*dyad)(interpreter_t *, verb_t *, value_t *, value_t *);
  598. };
  599. typedef struct {
  600. char *name;
  601. verb_t *(*adverb)(interpreter_t *, value_t *);
  602. verb_t *(*conjunction)(interpreter_t *, value_t *, value_t *);
  603. } adverb_t;
  604. struct _value_t {
  605. enum value_tag_t { ARRAY, VERB, SYMBOL, NUMBER, CHAR, NIL, UDF } tag;
  606. union {
  607. list_t *array;
  608. verb_t *verb;
  609. char *symbol;
  610. double number;
  611. unsigned char _char;
  612. } val;
  613. };
  614. verb_t *verb_new() {
  615. verb_t *verb = malloc_checked(sizeof(verb_t));
  616. return verb;
  617. }
  618. value_t *value_new(enum value_tag_t tag) {
  619. value_t *val;
  620. if (tag > SYMBOL)
  621. val = malloc_checked_atomic(sizeof(value_t));
  622. else
  623. val = malloc_checked(sizeof(value_t));
  624. val->tag = tag;
  625. return val;
  626. }
  627. value_t *value_new_const(enum value_tag_t tag) {
  628. value_t *val = malloc_checked_uncollectable(sizeof(value_t));
  629. val->tag = tag;
  630. return val;
  631. }
  632. value_t *_UNIT;
  633. value_t *value_new_array(list_t *array) {
  634. if (!array->data) {
  635. GC_FREE(array);
  636. return _UNIT;
  637. }
  638. value_t *val = value_new(ARRAY);
  639. val->val.array = array;
  640. return val;
  641. }
  642. table_t *VCACHE;
  643. value_t *value_new_verb(verb_t *verb) {
  644. value_t *val;
  645. if ((val = table_get(VCACHE, verb->name)))
  646. return val;
  647. val = value_new(VERB);
  648. val->val.verb = verb;
  649. return val;
  650. }
  651. table_t *SCACHE;
  652. value_t *value_new_symbol(char *symbol) {
  653. value_t *val;
  654. if ((val = table_get(SCACHE, symbol)))
  655. return val;
  656. val = value_new_const(SYMBOL);
  657. val->val.symbol = symbol;
  658. table_set(SCACHE, symbol, val);
  659. return val;
  660. }
  661. value_t *_NAN, *INF, *NINF;
  662. value_t *NNUMS[8];
  663. value_t *NUMS[256];
  664. value_t *CHARS[256];
  665. value_t *value_new_number(double number) {
  666. if (isnan(number))
  667. return _NAN;
  668. else if (number == INFINITY)
  669. return INF;
  670. else if (number == -INFINITY)
  671. return NINF;
  672. else if (number >= 0 && number < 256 && number == (double)((size_t)number))
  673. return NUMS[(size_t)number];
  674. else if (number < 0 && number >= -8 &&
  675. fabs(number) == (double)((size_t)fabs(number)))
  676. return NNUMS[((size_t)fabs(number)) - 1];
  677. value_t *val = value_new(NUMBER);
  678. val->val.number = number;
  679. return val;
  680. }
  681. value_t *value_new_char(unsigned char _char) { return CHARS[_char]; }
  682. bool value_equals(value_t *x, value_t *y) {
  683. if (x->tag != y->tag)
  684. return false;
  685. switch (x->tag) {
  686. case ARRAY: {
  687. list_t *tx = x->val.array;
  688. list_t *ty = y->val.array;
  689. if (tx->length == 0 && ty->length == 0)
  690. break;
  691. if (tx->length != ty->length)
  692. return false;
  693. for (size_t i = 0; i < tx->length; i++)
  694. if (!value_equals(tx->data[i], ty->data[i]))
  695. return false;
  696. }
  697. case VERB:
  698. return strcmp(x->val.verb->name, x->val.verb->name) == 0;
  699. case SYMBOL:
  700. return strcmp(x->val.symbol, y->val.symbol) == 0;
  701. case NUMBER:
  702. if (isnan(x->val.number) && isnan(y->val.number))
  703. break;
  704. return x->val.number == y->val.number;
  705. case CHAR:
  706. return x == y;
  707. case NIL:
  708. case UDF:
  709. break;
  710. }
  711. return true;
  712. }
  713. bool is_char_array(list_t *a) {
  714. for (size_t i = 0; i < a->length; i++) {
  715. value_t *v = a->data[i];
  716. if (v->tag != CHAR || !isprint(v->val._char))
  717. return false;
  718. }
  719. return true;
  720. }
  721. bool is_bytes_array(list_t *a) {
  722. for (size_t i = 0; i < a->length; i++) {
  723. value_t *v = a->data[i];
  724. if (v->tag != CHAR)
  725. return false;
  726. }
  727. return true;
  728. }
  729. bool is_arrays_array(list_t *a) {
  730. for (size_t i = 0; i < a->length; i++) {
  731. value_t *v = a->data[i];
  732. if (v->tag != ARRAY)
  733. return false;
  734. }
  735. return true;
  736. }
  737. bool is_not_arrays_array(list_t *a) {
  738. if (!a->data)
  739. return true;
  740. for (size_t i = 1; i < a->length; i++) {
  741. value_t *v = a->data[i];
  742. if (v->tag == ARRAY)
  743. return false;
  744. }
  745. return true;
  746. }
  747. bool is_matrix(list_t *a) {
  748. if (a->length < 2)
  749. return false;
  750. size_t rwl = ((value_t *)a->data[0])->val.array->length;
  751. if (rwl < 1)
  752. return false;
  753. for (size_t i = 0; i < a->length; i++) {
  754. value_t *v = a->data[i];
  755. if (v->tag != ARRAY || v->val.array->length != rwl ||
  756. !is_not_arrays_array(v->val.array) || is_char_array(v->val.array))
  757. return false;
  758. }
  759. return true;
  760. }
  761. char *value_show(value_t *v);
  762. char *show_array(value_t *v) {
  763. if (v->tag != ARRAY)
  764. return value_show(v);
  765. list_t *t = v->val.array;
  766. if (!t->data)
  767. return strdup_checked("()");
  768. buffer_t *buf = buffer_new();
  769. if (t->length == 1) {
  770. buffer_append(buf, ',');
  771. char *ts = value_show(t->data[0]);
  772. buffer_append_str(buf, ts);
  773. GC_FREE(ts);
  774. return buffer_read(buf);
  775. }
  776. if (is_char_array(t)) {
  777. for (size_t i = 0; i < t->length; i++)
  778. buffer_append(buf, ((value_t *)t->data[i])->val._char);
  779. return buffer_read(buf);
  780. }
  781. if (!is_arrays_array(t))
  782. for (size_t i = 0; i < t->length; i++) {
  783. char *ts = value_show(t->data[i]);
  784. buffer_append_str(buf, ts);
  785. GC_FREE(ts);
  786. if (i != t->length - 1)
  787. buffer_append(buf, ' ');
  788. }
  789. else if (is_matrix(t)) {
  790. size_t rwl = 0;
  791. size_t pad = 0;
  792. size_t padl = 0;
  793. list_t *ss = list_new();
  794. for (size_t i = 0; i < t->length; i++) {
  795. value_t *rw = t->data[i];
  796. list_t *rwt = rw->val.array;
  797. if (rwl < 1)
  798. rwl = rwt->length;
  799. for (size_t j = 0; j < rwt->length; j++) {
  800. char *s = value_show(rwt->data[j]);
  801. size_t z = strlen(s);
  802. if (z > pad)
  803. pad = z;
  804. if (j == 0 && z > padl)
  805. padl = z;
  806. list_push(ss, s);
  807. }
  808. }
  809. size_t k = 0;
  810. for (size_t i = 0; i < ss->length; i++) {
  811. char *s = ss->data[i];
  812. size_t mp = (k == 0 ? padl : pad) - strlen(s);
  813. while (mp--)
  814. buffer_append(buf, ' ');
  815. buffer_append_str(buf, s);
  816. GC_FREE(s);
  817. if (i != ss->length - 1) {
  818. if (k == rwl - 1) {
  819. k = 0;
  820. buffer_append(buf, '\n');
  821. } else {
  822. buffer_append(buf, ' ');
  823. k++;
  824. }
  825. }
  826. }
  827. GC_FREE(ss->data);
  828. GC_FREE(ss);
  829. } else
  830. for (size_t i = 0; i < t->length; i++) {
  831. value_t *rw = t->data[i];
  832. char *ts = show_array(rw);
  833. buffer_append_str(buf, ts);
  834. GC_FREE(ts);
  835. if (i != t->length - 1)
  836. buffer_append(buf, '\n');
  837. }
  838. return buffer_read(buf);
  839. }
  840. char *value_show(value_t *v) {
  841. switch (v->tag) {
  842. case ARRAY:
  843. return show_array(v);
  844. case VERB:
  845. return strdup_checked(v->val.verb->name);
  846. case SYMBOL:
  847. return strdup_checked(v->val.symbol);
  848. case NUMBER: {
  849. char buf[128];
  850. snprintf(buf, sizeof(buf), "%.15g", v->val.number);
  851. return strdup_checked(buf);
  852. }
  853. case CHAR: {
  854. if (!isprint(v->val._char)) {
  855. char buf[16];
  856. snprintf(buf, sizeof(buf), "4t.%d", v->val._char);
  857. return strdup_checked(buf);
  858. }
  859. char buf[2];
  860. buf[0] = v->val._char;
  861. buf[1] = 0;
  862. return strdup_checked(buf);
  863. }
  864. case NIL:
  865. return strdup_checked("nil");
  866. case UDF:
  867. return strdup_checked("udf");
  868. }
  869. return strdup_checked("<?>");
  870. }
  871. char *value_str(value_t *v) {
  872. if (v->tag == ARRAY && v->val.array->length == 1 &&
  873. ((value_t *)v->val.array->data[0])->tag == CHAR)
  874. return value_show(v->val.array->data[0]);
  875. else if (v->tag == ARRAY && is_bytes_array(v->val.array)) {
  876. buffer_t *buf = buffer_new();
  877. for (size_t i = 0; i < v->val.array->length; i++)
  878. buffer_append(buf, ((value_t *)v->val.array->data[i])->val._char);
  879. return buffer_read(buf);
  880. }
  881. return value_show(v);
  882. }
  883. double get_numeric(value_t *v) {
  884. if (v->tag == CHAR)
  885. return v->val._char;
  886. return v->val.number;
  887. }
  888. bool value_is_truthy(value_t *x) {
  889. switch (x->tag) {
  890. case ARRAY:
  891. return x->val.array->length != 0;
  892. case NUMBER:
  893. case CHAR:
  894. return get_numeric(x) != 0;
  895. case NIL:
  896. case UDF:
  897. return false;
  898. default:
  899. return true;
  900. }
  901. }
  902. verb_t *find_verb(char *s);
  903. interpreter_t *interpreter_new(void) {
  904. interpreter_t *state = malloc_checked(sizeof(interpreter_t));
  905. state->env = table_new();
  906. state->args = list_new();
  907. state->selfrefs = list_new();
  908. state->nil = value_new(NIL);
  909. state->udf = value_new(UDF);
  910. state->unit = _UNIT;
  911. state->at = find_verb("@");
  912. return state;
  913. }
  914. value_t *each_rank(interpreter_t *state, verb_t *f, value_t *x, unsigned int d,
  915. unsigned int rm) {
  916. if (!f->monad)
  917. return state->udf;
  918. if (rec_depth >= max_rec_depth)
  919. fatal("recursion-depth");
  920. rec_depth++;
  921. if (d >= rm || x->tag != ARRAY) {
  922. if (f->mark)
  923. list_push(state->selfrefs, f);
  924. value_t *r = f->monad(state, f, x);
  925. if (f->mark)
  926. list_pop(state->selfrefs);
  927. if (rec_depth > 0) rec_depth--;
  928. return r;
  929. }
  930. list_t *t = x->val.array;
  931. if (!t->data) {
  932. if (rec_depth > 0) rec_depth--;
  933. return x;
  934. }
  935. list_t *l = list_newk(t->length);
  936. for (size_t i = 0; i < t->length; i++)
  937. l->data[i] = each_rank(state, f, t->data[i], d + 1, rm);
  938. if (rec_depth > 0) rec_depth--;
  939. return value_new_array(l);
  940. }
  941. value_t *verb_at(interpreter_t *state, verb_t *self, value_t *x, value_t *y);
  942. value_t *apply_monad(interpreter_t *state, value_t *f, value_t *x) {
  943. if (f->tag == ARRAY)
  944. return verb_at(state, NULL, f, x);
  945. if (f->tag != VERB)
  946. return state->udf;
  947. if (!f->val.verb->monad)
  948. return state->udf;
  949. return each_rank(state, f->val.verb, x, 0, f->val.verb->rank[0]);
  950. }
  951. value_t *together(interpreter_t *state, verb_t *f, value_t *x, value_t *y,
  952. unsigned int dl, unsigned int dr, unsigned int rl,
  953. unsigned int rr) {
  954. if (!f->dyad)
  955. return state->udf;
  956. if (rec_depth >= max_rec_depth)
  957. fatal("recursion-depth");
  958. rec_depth++;
  959. if (dl >= rl && dr >= rr) {
  960. if (f->mark)
  961. list_push(state->selfrefs, f);
  962. value_t *r = f->dyad(state, f, x, y);
  963. if (f->mark)
  964. list_pop(state->selfrefs);
  965. if (rec_depth > 0) rec_depth--;
  966. return r;
  967. }
  968. if (dl < rl && dr < rr && x->tag == ARRAY && y->tag == ARRAY) {
  969. list_t *tx = x->val.array;
  970. list_t *ty = y->val.array;
  971. if (!tx->data || !ty->data) {
  972. if (rec_depth > 0) rec_depth--;
  973. return !tx->data ? x : y;
  974. }
  975. list_t *t = list_newk(ty->length < tx->length ? ty->length : tx->length);
  976. for (size_t i = 0; i < tx->length; i++) {
  977. if (i >= ty->length)
  978. break;
  979. t->data[i] =
  980. together(state, f, tx->data[i], ty->data[i], dl + 1, dr + 1, rl, rr);
  981. }
  982. if (rec_depth > 0) rec_depth--;
  983. return value_new_array(t);
  984. } else if ((x->tag != ARRAY || dl >= rl) && y->tag == ARRAY && dr < rr) {
  985. list_t *ty = y->val.array;
  986. if (!ty->data) {
  987. if (rec_depth > 0) rec_depth--;
  988. return y;
  989. }
  990. list_t *t = list_newk(ty->length);
  991. for (size_t i = 0; i < ty->length; i++)
  992. t->data[i] = together(state, f, x, ty->data[i], dl, dr + 1, rl, rr);
  993. return value_new_array(t);
  994. } else if ((y->tag != ARRAY || dr >= rr) && x->tag == ARRAY && dl < rl) {
  995. list_t *tx = x->val.array;
  996. if (!tx->data) {
  997. if (rec_depth > 0) rec_depth--;
  998. return x;
  999. }
  1000. list_t *t = list_newk(tx->length);
  1001. for (size_t i = 0; i < tx->length; i++)
  1002. t->data[i] = together(state, f, tx->data[i], y, dl + 1, dr, rl, rr);
  1003. if (rec_depth > 0) rec_depth--;
  1004. return value_new_array(t);
  1005. }
  1006. if (f->mark)
  1007. list_push(state->selfrefs, f);
  1008. value_t *r = f->dyad(state, f, x, y);
  1009. if (f->mark)
  1010. list_pop(state->selfrefs);
  1011. if (rec_depth > 0) rec_depth--;
  1012. return r;
  1013. }
  1014. value_t *apply_dyad(interpreter_t *state, value_t *f, value_t *x, value_t *y) {
  1015. if (f->tag == ARRAY)
  1016. return verb_at(state, NULL, verb_at(state, NULL, f, x), y);
  1017. if (f->tag != VERB)
  1018. return state->nil;
  1019. return together(state, f->val.verb, x, y, 0, 0, f->val.verb->rank[1],
  1020. f->val.verb->rank[2]);
  1021. }
  1022. typedef struct _node_t node_t;
  1023. struct _node_t {
  1024. enum node_tag_t {
  1025. N_STRAND,
  1026. N_LITERAL,
  1027. N_INDEX1,
  1028. N_INDEX2,
  1029. N_FUN,
  1030. N_MONAD,
  1031. N_DYAD,
  1032. N_ADV,
  1033. N_CONJ,
  1034. N_PARTIAL_CONJ,
  1035. N_FORK,
  1036. N_HOOK,
  1037. N_BOND,
  1038. N_OVER,
  1039. N_BIND
  1040. } tag;
  1041. adverb_t *av;
  1042. value_t *v;
  1043. list_t *l;
  1044. node_t *a;
  1045. node_t *b;
  1046. node_t *c;
  1047. size_t dp;
  1048. };
  1049. char *node_show(node_t *n) {
  1050. switch (n->tag) {
  1051. case N_STRAND: {
  1052. buffer_t *buf = buffer_new();
  1053. for (size_t i = 0; i < n->l->length; i++) {
  1054. if (i != 0)
  1055. buffer_append_str(buf, ",:");
  1056. char *s = node_show(n->l->data[i]);
  1057. buffer_append_str(buf, s);
  1058. GC_FREE(s);
  1059. }
  1060. return buffer_read(buf);
  1061. }
  1062. case N_LITERAL:
  1063. return value_show(n->v);
  1064. case N_INDEX1: {
  1065. char *s;
  1066. buffer_t *buf = buffer_new();
  1067. s = node_show(n->a);
  1068. buffer_append_str(buf, s);
  1069. GC_FREE(s);
  1070. buffer_append(buf, ' ');
  1071. s = node_show(n->b);
  1072. buffer_append_str(buf, s);
  1073. GC_FREE(s);
  1074. return buffer_read(buf);
  1075. }
  1076. case N_INDEX2: {
  1077. char *s;
  1078. buffer_t *buf = buffer_new();
  1079. s = node_show(n->a);
  1080. buffer_append_str(buf, s);
  1081. GC_FREE(s);
  1082. buffer_append(buf, ' ');
  1083. s = node_show(n->b);
  1084. buffer_append_str(buf, s);
  1085. GC_FREE(s);
  1086. buffer_append(buf, ' ');
  1087. s = node_show(n->c);
  1088. buffer_append_str(buf, s);
  1089. GC_FREE(s);
  1090. return buffer_read(buf);
  1091. }
  1092. case N_FUN: {
  1093. buffer_t *buf = buffer_new();
  1094. buffer_append(buf, ':');
  1095. char *s = node_show(n->a);
  1096. buffer_append_str(buf, s);
  1097. GC_FREE(s);
  1098. return buffer_read(buf);
  1099. }
  1100. case N_MONAD:
  1101. case N_HOOK:
  1102. case N_BOND:
  1103. case N_OVER: {
  1104. char *s;
  1105. buffer_t *buf = buffer_new();
  1106. s = node_show(n->a);
  1107. buffer_append_str(buf, s);
  1108. GC_FREE(s);
  1109. s = node_show(n->b);
  1110. buffer_append_str(buf, s);
  1111. GC_FREE(s);
  1112. return buffer_read(buf);
  1113. }
  1114. case N_DYAD: {
  1115. char *s;
  1116. buffer_t *buf = buffer_new();
  1117. s = node_show(n->b);
  1118. buffer_append_str(buf, s);
  1119. GC_FREE(s);
  1120. s = node_show(n->a);
  1121. buffer_append_str(buf, s);
  1122. GC_FREE(s);
  1123. s = node_show(n->c);
  1124. buffer_append_str(buf, s);
  1125. GC_FREE(s);
  1126. return buffer_read(buf);
  1127. }
  1128. case N_ADV:
  1129. case N_PARTIAL_CONJ: {
  1130. buffer_t *buf = buffer_new();
  1131. char *s = node_show(n->a);
  1132. buffer_append_str(buf, s);
  1133. GC_FREE(s);
  1134. buffer_append_str(buf, n->av->name);
  1135. return buffer_read(buf);
  1136. }
  1137. case N_CONJ: {
  1138. char *s;
  1139. buffer_t *buf = buffer_new();
  1140. s = node_show(n->a);
  1141. buffer_append_str(buf, s);
  1142. GC_FREE(s);
  1143. buffer_append_str(buf, n->av->name);
  1144. s = node_show(n->b);
  1145. buffer_append_str(buf, s);
  1146. GC_FREE(s);
  1147. return buffer_read(buf);
  1148. }
  1149. case N_FORK: {
  1150. char *s;
  1151. buffer_t *buf = buffer_new();
  1152. s = node_show(n->a);
  1153. buffer_append_str(buf, s);
  1154. GC_FREE(s);
  1155. s = node_show(n->b);
  1156. buffer_append_str(buf, s);
  1157. GC_FREE(s);
  1158. s = node_show(n->c);
  1159. buffer_append_str(buf, s);
  1160. GC_FREE(s);
  1161. return buffer_read(buf);
  1162. }
  1163. case N_BIND: {
  1164. char *s;
  1165. buffer_t *buf = buffer_new();
  1166. s = node_show(n->a);
  1167. buffer_append_str(buf, s);
  1168. GC_FREE(s);
  1169. buffer_append(buf, ':');
  1170. s = node_show(n->b);
  1171. buffer_append_str(buf, s);
  1172. GC_FREE(s);
  1173. return buffer_read(buf);
  1174. }
  1175. }
  1176. return strdup_checked("<?>");
  1177. }
  1178. value_t *_fork_monad(interpreter_t *state, verb_t *self, value_t *x) {
  1179. verb_t *f = list_index(self->bonds, 0);
  1180. verb_t *g = list_index(self->bonds, 1);
  1181. verb_t *h = list_index(self->bonds, 2);
  1182. value_t *l = each_rank(state, f, x, 0, f->rank[0]);
  1183. value_t *r = each_rank(state, h, x, 0, f->rank[0]);
  1184. return together(state, g, l, r, 0, 0, g->rank[1], g->rank[2]);
  1185. }
  1186. value_t *_fork_dyad(interpreter_t *state, verb_t *self, value_t *x,
  1187. value_t *y) {
  1188. verb_t *f = list_index(self->bonds, 0);
  1189. verb_t *g = list_index(self->bonds, 1);
  1190. verb_t *h = list_index(self->bonds, 2);
  1191. value_t *l = each_rank(state, f, x, 0, f->rank[0]);
  1192. value_t *r = each_rank(state, h, y, 0, f->rank[0]);
  1193. return together(state, g, l, r, 0, 0, g->rank[1], g->rank[2]);
  1194. }
  1195. value_t *_hook_monad(interpreter_t *state, verb_t *self, value_t *x) {
  1196. verb_t *f = list_index(self->bonds, 0);
  1197. verb_t *g = list_index(self->bonds, 1);
  1198. value_t *r = each_rank(state, g, x, 0, g->rank[0]);
  1199. return each_rank(state, f, r, 0, f->rank[0]);
  1200. }
  1201. value_t *_hook_dyad(interpreter_t *state, verb_t *self, value_t *x,
  1202. value_t *y) {
  1203. verb_t *f = list_index(self->bonds, 0);
  1204. verb_t *g = list_index(self->bonds, 1);
  1205. value_t *r = together(state, g, x, y, 0, 0, g->rank[1], g->rank[2]);
  1206. return each_rank(state, f, r, 0, f->rank[0]);
  1207. }
  1208. value_t *_bond_monad(interpreter_t *state, verb_t *self, value_t *x) {
  1209. verb_t *f = list_index(self->bonds, 0);
  1210. value_t *g = list_index(self->bonds, 1);
  1211. return together(state, f, g, x, 0, 0, f->rank[1], f->rank[2]);
  1212. }
  1213. value_t *_bond_dyad(interpreter_t *state, verb_t *self, value_t *x,
  1214. value_t *y) {
  1215. verb_t *f = list_index(self->bonds, 0);
  1216. value_t *g = list_index(self->bonds, 1);
  1217. value_t *r = together(state, f, x, y, 0, 0, f->rank[1], f->rank[2]);
  1218. return together(state, f, g, r, 0, 0, f->rank[1], f->rank[2]);
  1219. }
  1220. value_t *_over_monad(interpreter_t *state, verb_t *self, value_t *x) {
  1221. value_t *f = list_index(self->bonds, 0);
  1222. verb_t *g = list_index(self->bonds, 1);
  1223. verb_t *h = list_index(self->bonds, 2);
  1224. value_t *l = each_rank(state, h, x, 0, h->rank[0]);
  1225. return together(state, g, f, l, 0, 0, g->rank[1], g->rank[2]);
  1226. }
  1227. value_t *_over_dyad(interpreter_t *state, verb_t *self, value_t *x,
  1228. value_t *y) {
  1229. value_t *f = list_index(self->bonds, 0);
  1230. verb_t *g = list_index(self->bonds, 1);
  1231. verb_t *h = list_index(self->bonds, 2);
  1232. value_t *l = together(state, h, x, y, 0, 0, h->rank[1], h->rank[2]);
  1233. return together(state, g, f, l, 0, 0, g->rank[1], g->rank[2]);
  1234. }
  1235. bool function_collect_args(node_t *node, unsigned int *argc) {
  1236. if (!node)
  1237. return false;
  1238. if (node->tag == N_LITERAL && node->v->tag == SYMBOL &&
  1239. strcmp(node->v->val.symbol, "y") == 0) {
  1240. *argc = 2;
  1241. return true;
  1242. } else if (node->tag == N_LITERAL && node->v->tag == SYMBOL &&
  1243. strcmp(node->v->val.symbol, "x") == 0) {
  1244. if (*argc < 2)
  1245. *argc = 1;
  1246. } else if (node->tag == N_MONAD || node->tag == N_CONJ ||
  1247. node->tag == N_HOOK || node->tag == N_BOND ||
  1248. node->tag == N_INDEX1) {
  1249. if (function_collect_args(node->a, argc))
  1250. return true;
  1251. if (function_collect_args(node->b, argc))
  1252. return true;
  1253. } else if (node->tag == N_DYAD || node->tag == N_FORK ||
  1254. node->tag == N_OVER || node->tag == N_INDEX2) {
  1255. if (function_collect_args(node->a, argc))
  1256. return true;
  1257. if (function_collect_args(node->b, argc))
  1258. return true;
  1259. if (function_collect_args(node->c, argc))
  1260. return true;
  1261. } else if (node->tag == N_ADV) {
  1262. if (function_collect_args(node->a, argc))
  1263. return true;
  1264. } else if (node->tag == N_STRAND) {
  1265. list_t *t = node->l;
  1266. for (size_t i = 0; i < t->length; i++)
  1267. if (function_collect_args(t->data[i], argc))
  1268. return true;
  1269. } else if (node->tag == N_BIND) {
  1270. if (function_collect_args(node->b, argc))
  1271. return true;
  1272. }
  1273. return false;
  1274. }
  1275. value_t *interpreter_walk(interpreter_t *state, node_t *node);
  1276. value_t *_const_monad(interpreter_t *state, verb_t *self, value_t *x) {
  1277. return interpreter_walk(state, self->bonds->data[0]);
  1278. }
  1279. value_t *_const_dyad(interpreter_t *state, verb_t *self, value_t *x,
  1280. value_t *y) {
  1281. return interpreter_walk(state, self->bonds->data[0]);
  1282. }
  1283. value_t *_constv_monad(interpreter_t *state, verb_t *self, value_t *x) {
  1284. return self->bonds->data[0];
  1285. }
  1286. value_t *_constv_dyad(interpreter_t *state, verb_t *self, value_t *x,
  1287. value_t *y) {
  1288. return self->bonds->data[0];
  1289. }
  1290. value_t *_fun_monad(interpreter_t *state, verb_t *self, value_t *x) {
  1291. list_t *args = list_new();
  1292. list_push(args, x);
  1293. list_push(args, self);
  1294. list_push(state->args, args);
  1295. value_t *r = interpreter_walk(state, self->bonds->data[0]);
  1296. list_pop(state->args);
  1297. GC_FREE(args);
  1298. return r;
  1299. }
  1300. value_t *_fun_dyad(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1301. list_t *args = list_new();
  1302. list_push(args, x);
  1303. list_push(args, y);
  1304. list_push(args, self);
  1305. list_push(state->args, args);
  1306. value_t *r = interpreter_walk(state, self->bonds->data[1]);
  1307. list_pop(state->args);
  1308. GC_FREE(args);
  1309. return r;
  1310. }
  1311. value_t *_partial_conjunction(interpreter_t *state, verb_t *self, value_t *x) {
  1312. adverb_t *av = self->bonds->data[0];
  1313. value_t *a = self->bonds->data[1];
  1314. return value_new_verb(av->conjunction(state, a, x));
  1315. }
  1316. node_t *node_new1(enum node_tag_t tag, node_t *a);
  1317. value_t *interpreter_walk(interpreter_t *state, node_t *node) {
  1318. if (!node)
  1319. return state->nil;
  1320. switch (node->tag) {
  1321. case N_STRAND: {
  1322. list_t *t = list_copy(node->l);
  1323. for (size_t i = 0; i < t->length; i++)
  1324. t->data[i] = interpreter_walk(state, t->data[i]);
  1325. return value_new_array(t);
  1326. }
  1327. case N_LITERAL: {
  1328. value_t *v = node->v;
  1329. value_t *t = NULL;
  1330. if (v->tag == SYMBOL) {
  1331. char *n = v->val.symbol;
  1332. if (state->args->data) {
  1333. list_t *args = list_index(state->args, -1);
  1334. size_t argc = args->length - 1;
  1335. if (argc == 2 && strcmp(n, "y") == 0)
  1336. return args->data[1];
  1337. else if (strcmp(n, "x") == 0)
  1338. return args->data[0];
  1339. }
  1340. if ((t = table_get(state->env, n)))
  1341. return t;
  1342. if (strcmp(n, "T") == 0)
  1343. return value_new_number(time(NULL));
  1344. }
  1345. return v;
  1346. }
  1347. case N_INDEX1:
  1348. return together(state, state->at, interpreter_walk(state, node->a),
  1349. interpreter_walk(state, node->b), 0, 0, state->at->rank[1],
  1350. state->at->rank[2]);
  1351. case N_INDEX2:
  1352. return together(state, state->at,
  1353. together(state, state->at, interpreter_walk(state, node->a),
  1354. interpreter_walk(state, node->b), 0, 0,
  1355. state->at->rank[1], state->at->rank[2]),
  1356. interpreter_walk(state, node->c), 0, 0, state->at->rank[1],
  1357. state->at->rank[2]);
  1358. case N_FUN: {
  1359. unsigned int argc = 0;
  1360. function_collect_args(node->a, &argc);
  1361. verb_t *nv = verb_new();
  1362. if (argc > 0)
  1363. nv->is_fun = true;
  1364. nv->bonds = list_new();
  1365. char *s = node_show(node->a);
  1366. size_t z = strlen(s) + 2;
  1367. nv->name = malloc_checked_atomic(z);
  1368. snprintf(nv->name, z, ":%s", s);
  1369. GC_FREE(s);
  1370. nv->rank[0] = 0;
  1371. nv->rank[1] = 0;
  1372. nv->rank[2] = 0;
  1373. if (argc == 0) {
  1374. list_push(nv->bonds, node->a);
  1375. nv->monad = _const_monad;
  1376. nv->dyad = _const_dyad;
  1377. } else if (argc == 1) {
  1378. list_push(nv->bonds, node->a);
  1379. nv->monad = _fun_monad;
  1380. nv->dyad = NULL;
  1381. } else {
  1382. nv->monad = NULL;
  1383. list_push(nv->bonds, state->udf);
  1384. list_push(nv->bonds, node->a);
  1385. nv->dyad = _fun_dyad;
  1386. }
  1387. return value_new_verb(nv);
  1388. }
  1389. case N_MONAD:
  1390. return apply_monad(state, interpreter_walk(state, node->a),
  1391. interpreter_walk(state, node->b));
  1392. case N_DYAD:
  1393. return apply_dyad(state, interpreter_walk(state, node->a),
  1394. interpreter_walk(state, node->b),
  1395. interpreter_walk(state, node->c));
  1396. case N_ADV: {
  1397. value_t *v = interpreter_walk(state, node->a);
  1398. verb_t *nv = node->av->adverb(state, v);
  1399. if (node->dp < 2)
  1400. nv->mark = true;
  1401. return value_new_verb(nv);
  1402. }
  1403. case N_CONJ: {
  1404. value_t *v1 = interpreter_walk(state, node->a);
  1405. value_t *v2 = interpreter_walk(state, node->b);
  1406. verb_t *nv = node->av->conjunction(state, v1, v2);
  1407. if (node->dp < 2)
  1408. nv->mark = true;
  1409. return value_new_verb(nv);
  1410. }
  1411. case N_PARTIAL_CONJ: {
  1412. verb_t *nv = verb_new();
  1413. value_t *a = interpreter_walk(state, node->a);
  1414. char *r = value_show(a);
  1415. size_t l = strlen(r) + strlen(node->av->name) + 1;
  1416. nv->name = malloc_checked_atomic(l);
  1417. snprintf(nv->name, l, "%s%s", r, node->av->name);
  1418. GC_FREE(r);
  1419. nv->bonds = list_new();
  1420. list_push(nv->bonds, node->av);
  1421. list_push(nv->bonds, a);
  1422. nv->rank[0] = 0;
  1423. nv->rank[1] = 0;
  1424. nv->rank[2] = 0;
  1425. nv->monad = _partial_conjunction;
  1426. nv->dyad = NULL;
  1427. if (node->dp < 2)
  1428. nv->mark = true;
  1429. return value_new_verb(nv);
  1430. }
  1431. case N_FORK: {
  1432. value_t *_f = interpreter_walk(state, node->a);
  1433. if (_f->tag != VERB)
  1434. return state->udf;
  1435. value_t *_g = interpreter_walk(state, node->b);
  1436. if (_g->tag != VERB)
  1437. return state->udf;
  1438. value_t *_h = interpreter_walk(state, node->c);
  1439. if (_h->tag != VERB)
  1440. return state->udf;
  1441. verb_t *f = _f->val.verb;
  1442. verb_t *g = _g->val.verb;
  1443. verb_t *h = _h->val.verb;
  1444. verb_t *nv = verb_new();
  1445. nv->bonds = list_newk(3);
  1446. nv->bonds->data[0] = f;
  1447. nv->bonds->data[1] = g;
  1448. nv->bonds->data[2] = h;
  1449. size_t l = strlen(f->name) + strlen(g->name) + strlen(h->name) + 1;
  1450. nv->name = malloc_checked_atomic(l);
  1451. snprintf(nv->name, l, "%s%s%s", f->name, g->name, h->name);
  1452. nv->rank[0] = 0;
  1453. nv->rank[1] = 0;
  1454. nv->rank[2] = 0;
  1455. nv->monad = _fork_monad;
  1456. nv->dyad = _fork_dyad;
  1457. if (node->dp < 2)
  1458. nv->mark = true;
  1459. return value_new_verb(nv);
  1460. }
  1461. case N_HOOK: {
  1462. value_t *_f = interpreter_walk(state, node->a);
  1463. if (_f->tag != VERB)
  1464. return state->udf;
  1465. value_t *_g = interpreter_walk(state, node->b);
  1466. if (_g->tag != VERB)
  1467. return state->udf;
  1468. verb_t *f = _f->val.verb;
  1469. verb_t *g = _g->val.verb;
  1470. verb_t *nv = verb_new();
  1471. nv->bonds = list_newk(2);
  1472. nv->bonds->data[0] = f;
  1473. nv->bonds->data[1] = g;
  1474. size_t l = strlen(f->name) + strlen(g->name) + 1;
  1475. nv->name = malloc_checked_atomic(l);
  1476. snprintf(nv->name, l, "%s%s", f->name, g->name);
  1477. nv->rank[0] = 0;
  1478. nv->rank[1] = 0;
  1479. nv->rank[2] = 0;
  1480. nv->monad = _hook_monad;
  1481. nv->dyad = _hook_dyad;
  1482. if (node->dp < 2)
  1483. nv->mark = true;
  1484. return value_new_verb(nv);
  1485. }
  1486. case N_BOND: {
  1487. value_t *_f = interpreter_walk(state, node->a);
  1488. if (_f->tag != VERB)
  1489. return state->udf;
  1490. value_t *g = interpreter_walk(state, node->b);
  1491. verb_t *f = _f->val.verb;
  1492. verb_t *nv = verb_new();
  1493. nv->bonds = list_newk(2);
  1494. nv->bonds->data[0] = f;
  1495. nv->bonds->data[1] = g;
  1496. char *r = value_show(g);
  1497. size_t l = strlen(r) + strlen(f->name) + 1;
  1498. nv->name = malloc_checked_atomic(l);
  1499. snprintf(nv->name, l, "%s%s", r, f->name);
  1500. GC_FREE(r);
  1501. nv->rank[0] = 0;
  1502. nv->rank[1] = 0;
  1503. nv->rank[2] = 0;
  1504. nv->monad = _bond_monad;
  1505. nv->dyad = _bond_dyad;
  1506. if (node->dp < 2)
  1507. nv->mark = true;
  1508. return value_new_verb(nv);
  1509. }
  1510. case N_OVER: {
  1511. value_t *f = interpreter_walk(state, node->a);
  1512. value_t *_g = interpreter_walk(state, node->b);
  1513. if (_g->tag != VERB)
  1514. return state->udf;
  1515. value_t *_h = interpreter_walk(state, node->c);
  1516. if (_h->tag != VERB)
  1517. return state->udf;
  1518. verb_t *g = _g->val.verb;
  1519. verb_t *h = _h->val.verb;
  1520. verb_t *nv = verb_new();
  1521. nv->bonds = list_newk(3);
  1522. nv->bonds->data[0] = f;
  1523. nv->bonds->data[1] = g;
  1524. nv->bonds->data[2] = h;
  1525. char *r = value_show(f);
  1526. size_t l = strlen(r) + strlen(g->name) + strlen(h->name) + 1;
  1527. nv->name = malloc_checked_atomic(l);
  1528. snprintf(nv->name, l, "%s%s%s", r, g->name, h->name);
  1529. GC_FREE(r);
  1530. nv->rank[0] = 0;
  1531. nv->rank[1] = 0;
  1532. nv->rank[2] = 0;
  1533. nv->monad = _over_monad;
  1534. nv->dyad = _over_dyad;
  1535. if (node->dp < 2)
  1536. nv->mark = true;
  1537. return value_new_verb(nv);
  1538. }
  1539. case N_BIND: {
  1540. value_t *l = node->a->v;
  1541. node_t *b = node->b;
  1542. if (state->bn || state->args->data || node->dp != 0) {
  1543. table_set(state->env, l->val.symbol, interpreter_walk(state, b));
  1544. break;
  1545. }
  1546. unsigned int argc = 0;
  1547. function_collect_args(b, &argc);
  1548. if (argc != 0)
  1549. b = node_new1(N_FUN, b);
  1550. bool t = state->bn;
  1551. state->bn = true;
  1552. value_t *r = interpreter_walk(state, b);
  1553. state->bn = t;
  1554. if (argc != 0) {
  1555. GC_FREE(r->val.verb->name);
  1556. r->val.verb->name = l->val.symbol;
  1557. }
  1558. if (r->tag == VERB && argc == 0)
  1559. r->val.verb->mark = true;
  1560. value_t *ov = table_get(state->env, l->val.symbol);
  1561. if (ov && ov->tag == VERB && ov->val.verb->is_fun && r->tag == VERB &&
  1562. r->val.verb->is_fun) {
  1563. if (!ov->val.verb->monad && r->val.verb->monad) {
  1564. list_set(ov->val.verb->bonds, 0, r->val.verb->bonds->data[0]);
  1565. ov->val.verb->monad = r->val.verb->monad;
  1566. break;
  1567. }
  1568. if (!ov->val.verb->dyad && r->val.verb->dyad) {
  1569. list_push(ov->val.verb->bonds, r->val.verb->bonds->data[1]);
  1570. ov->val.verb->dyad = r->val.verb->dyad;
  1571. break;
  1572. }
  1573. }
  1574. table_set(state->env, l->val.symbol, r);
  1575. } break;
  1576. }
  1577. return state->nil;
  1578. }
  1579. value_t *verb_const(interpreter_t *state, verb_t *self, value_t *x) {
  1580. verb_t *nv = verb_new();
  1581. nv->bonds = list_newk(1);
  1582. nv->bonds->data[0] = x;
  1583. char *r = value_show(x);
  1584. size_t l = strlen(r) + 2;
  1585. nv->name = malloc_checked_atomic(l);
  1586. snprintf(nv->name, l, ":%s", r);
  1587. nv->rank[0] = 0;
  1588. nv->rank[1] = 0;
  1589. nv->rank[2] = 0;
  1590. nv->monad = _constv_monad;
  1591. nv->dyad = _constv_dyad;
  1592. return value_new_verb(nv);
  1593. }
  1594. value_t *verb_bind(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1595. if (x->tag == SYMBOL) {
  1596. if (y->tag == VERB)
  1597. y->val.verb->mark = true;
  1598. table_set(state->env, x->val.symbol, y);
  1599. }
  1600. return state->udf;
  1601. }
  1602. table_t *Inverses;
  1603. value_t *verb_unbind(interpreter_t *state, verb_t *self, value_t *x) {
  1604. if (x->tag == SYMBOL) {
  1605. table_delete(state->env, x->val.symbol);
  1606. return state->nil;
  1607. }
  1608. return state->udf;
  1609. }
  1610. value_t *verb_obverse(interpreter_t *state, verb_t *self, value_t *x,
  1611. value_t *y) {
  1612. if (x->tag == VERB && y->tag == VERB) {
  1613. verb_t *vx = x->val.verb;
  1614. if (!y->val.verb->monad)
  1615. return state->udf;
  1616. if (vx->is_fun)
  1617. return state->udf;
  1618. if (table_has(Inverses, vx->name))
  1619. return state->udf;
  1620. table_set(Inverses, vx->name, y->val.verb);
  1621. return state->nil;
  1622. }
  1623. return state->udf;
  1624. }
  1625. value_t *verb_flip(interpreter_t *state, verb_t *self, value_t *x) {
  1626. if (x->tag != ARRAY || !x->val.array->data)
  1627. return state->udf;
  1628. if (!is_arrays_array(x->val.array))
  1629. return state->udf;
  1630. list_t *r = list_new();
  1631. value_t *c0 = x->val.array->data[0];
  1632. list_t *c0t = c0->val.array;
  1633. size_t c0l = c0t->length;
  1634. for (size_t i = 0; i < c0l; i++) {
  1635. list_t *nc = list_new();
  1636. for (size_t j = 0; j < x->val.array->length; j++) {
  1637. value_t *rw = x->val.array->data[j];
  1638. list_t *rwt = rw->val.array;
  1639. if (!rwt->data)
  1640. return state->udf;
  1641. value_t *v = list_index(rwt, i);
  1642. if (!v)
  1643. v = rwt->data[0];
  1644. list_push(nc, v);
  1645. }
  1646. list_push(r, value_new_array(nc));
  1647. }
  1648. return value_new_array(r);
  1649. }
  1650. value_t *verb_plus(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1651. if ((x->tag == NUMBER || x->tag == CHAR) &&
  1652. (y->tag == NUMBER || y->tag == CHAR)) {
  1653. if (x->tag == CHAR || y->tag == CHAR)
  1654. return value_new_char(get_numeric(x) + get_numeric(y));
  1655. return value_new_number(get_numeric(x) + get_numeric(y));
  1656. }
  1657. return _NAN;
  1658. }
  1659. value_t *verb_sign(interpreter_t *state, verb_t *self, value_t *x) {
  1660. if (x->tag == NUMBER)
  1661. return x->val.number < 0 ? NNUMS[0] : x->val.number > 0 ? NUMS[1] : NUMS[0];
  1662. return _NAN;
  1663. }
  1664. double gcd(double a, double b) {
  1665. if (b != 0)
  1666. return gcd(b, fmod(a, b));
  1667. else
  1668. return fabs(a);
  1669. }
  1670. value_t *verb_gcd(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1671. if (x->tag == NUMBER && y->tag == NUMBER)
  1672. return value_new_number(gcd(x->val.number, y->val.number));
  1673. return _NAN;
  1674. }
  1675. value_t *verb_sin(interpreter_t *state, verb_t *self, value_t *x) {
  1676. if (x->tag == NUMBER)
  1677. return value_new_number(sin(x->val.number));
  1678. return _NAN;
  1679. }
  1680. value_t *verb_square(interpreter_t *state, verb_t *self, value_t *x) {
  1681. if (x->tag == NUMBER)
  1682. return value_new_number(x->val.number * x->val.number);
  1683. return _NAN;
  1684. }
  1685. value_t *verb_negate(interpreter_t *state, verb_t *self, value_t *x) {
  1686. if (x->tag == NUMBER)
  1687. return value_new_number(-x->val.number);
  1688. return _NAN;
  1689. }
  1690. value_t *verb_minus(interpreter_t *state, verb_t *self, value_t *x,
  1691. value_t *y) {
  1692. if ((x->tag == NUMBER || x->tag == CHAR) &&
  1693. (y->tag == NUMBER || y->tag == CHAR)) {
  1694. if (x->tag == CHAR || y->tag == CHAR)
  1695. return value_new_char(get_numeric(x) - get_numeric(y));
  1696. return value_new_number(get_numeric(x) - get_numeric(y));
  1697. }
  1698. return _NAN;
  1699. }
  1700. value_t *verb_atan(interpreter_t *state, verb_t *self, value_t *x) {
  1701. if (x->tag == NUMBER)
  1702. return value_new_number(atan(x->val.number));
  1703. return _NAN;
  1704. }
  1705. value_t *verb_atan2(interpreter_t *state, verb_t *self, value_t *x,
  1706. value_t *y) {
  1707. if (x->tag == NUMBER && y->tag == NUMBER)
  1708. return value_new_number(atan2(x->val.number, y->val.number));
  1709. return _NAN;
  1710. }
  1711. value_t *verb_first(interpreter_t *state, verb_t *self, value_t *x) {
  1712. if (x->tag != ARRAY)
  1713. return x;
  1714. if (!x->val.array->data)
  1715. return state->udf;
  1716. return x->val.array->data[0];
  1717. }
  1718. value_t *verb_times(interpreter_t *state, verb_t *self, value_t *x,
  1719. value_t *y) {
  1720. if ((x->tag == NUMBER || x->tag == CHAR) &&
  1721. (y->tag == NUMBER || y->tag == CHAR)) {
  1722. if (x->tag == CHAR || y->tag == CHAR)
  1723. return value_new_char(get_numeric(x) * get_numeric(y));
  1724. return value_new_number(get_numeric(x) * get_numeric(y));
  1725. }
  1726. return _NAN;
  1727. }
  1728. double lcm(double a, double b) { return (a * b) / gcd(a, b); }
  1729. uint64_t factorial(uint64_t n) {
  1730. uint64_t r = 1;
  1731. while (n > 0)
  1732. r *= n--;
  1733. return r;
  1734. }
  1735. value_t *verb_factorial(interpreter_t *state, verb_t *self, value_t *x) {
  1736. if (x->tag == NUMBER)
  1737. return value_new_number(factorial((uint64_t)fabs(x->val.number)));
  1738. return _NAN;
  1739. }
  1740. value_t *verb_lcm(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1741. if (x->tag == NUMBER && y->tag == NUMBER)
  1742. return value_new_number(lcm(x->val.number, y->val.number));
  1743. return _NAN;
  1744. }
  1745. value_t *verb_double(interpreter_t *state, verb_t *self, value_t *x) {
  1746. if (x->tag == NUMBER)
  1747. return value_new_number(x->val.number * 2);
  1748. return _NAN;
  1749. }
  1750. value_t *verb_replicate(interpreter_t *state, verb_t *self, value_t *x,
  1751. value_t *y) {
  1752. if (x->tag == NUMBER) {
  1753. size_t k = fabs(x->val.number);
  1754. list_t *r = list_new();
  1755. while (k--)
  1756. list_push(r, y);
  1757. return value_new_array(r);
  1758. }
  1759. return state->udf;
  1760. }
  1761. value_t *verb_reciprocal(interpreter_t *state, verb_t *self, value_t *x) {
  1762. if (x->tag == NUMBER)
  1763. return value_new_number(1 / x->val.number);
  1764. return _NAN;
  1765. }
  1766. value_t *verb_divide(interpreter_t *state, verb_t *self, value_t *x,
  1767. value_t *y) {
  1768. if (x->tag == NUMBER && y->tag == NUMBER) {
  1769. double ny = y->val.number;
  1770. if (ny == 0)
  1771. return INF;
  1772. return value_new_number(x->val.number / ny);
  1773. }
  1774. return _NAN;
  1775. }
  1776. double npower(double base, int n) {
  1777. if (n < 0)
  1778. return npower(1 / base, -n);
  1779. else if (n == 0)
  1780. return 1.0;
  1781. else if (n == 1)
  1782. return base;
  1783. else if (n % 2)
  1784. return base * npower(base * base, n / 2);
  1785. else
  1786. return npower(base * base, n / 2);
  1787. }
  1788. double nroot(double base, int n) {
  1789. if (n == 1)
  1790. return base;
  1791. else if (n <= 0 || base < 0)
  1792. return NAN;
  1793. else {
  1794. double delta, x = base / n;
  1795. do {
  1796. delta = (base / npower(x, n - 1) - x) / n;
  1797. x += delta;
  1798. } while (fabs(delta) >= 1e-8);
  1799. return x;
  1800. }
  1801. }
  1802. value_t *verb_sqrt(interpreter_t *state, verb_t *self, value_t *x) {
  1803. if (x->tag == NUMBER)
  1804. return value_new_number(sqrt(x->val.number));
  1805. return _NAN;
  1806. }
  1807. value_t *verb_root(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1808. if (x->tag == NUMBER && y->tag == NUMBER)
  1809. return value_new_number(nroot(y->val.number, x->val.number));
  1810. return _NAN;
  1811. }
  1812. value_t *verb_halve(interpreter_t *state, verb_t *self, value_t *x) {
  1813. if (x->tag == NUMBER)
  1814. return value_new_number(x->val.number / 2);
  1815. return _NAN;
  1816. }
  1817. value_t *verb_idivide(interpreter_t *state, verb_t *self, value_t *x,
  1818. value_t *y) {
  1819. if ((x->tag == NUMBER || x->tag == CHAR) &&
  1820. (y->tag == NUMBER || y->tag == CHAR)) {
  1821. double a = get_numeric(x);
  1822. double b = get_numeric(y);
  1823. if (x->tag == CHAR || y->tag == CHAR)
  1824. return b == 0 ? state->udf : value_new_char(fabs(trunc(a / b)));
  1825. if (b == 0)
  1826. return INF;
  1827. return value_new_number(trunc(a / b));
  1828. }
  1829. return _NAN;
  1830. }
  1831. value_t *verb_enlist(interpreter_t *state, verb_t *self, value_t *x);
  1832. value_t *verb_pred(interpreter_t *state, verb_t *self, value_t *x);
  1833. value_t *verb_range(interpreter_t *state, verb_t *self, value_t *x, value_t *y);
  1834. value_t *verb_enum(interpreter_t *state, verb_t *self, value_t *x) {
  1835. if (value_equals(x, NUMS[1]))
  1836. return verb_enlist(state, NULL, NUMS[0]);
  1837. else if (value_equals(x, NUMS[0]))
  1838. return state->unit;
  1839. return verb_range(state, self, NUMS[0], verb_pred(state, self, x));
  1840. }
  1841. value_t *verb_mod(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1842. if (x->tag == NUMBER && y->tag == NUMBER) {
  1843. double ny = y->val.number;
  1844. if (ny == 0)
  1845. return _NAN;
  1846. return value_new_number(fmod(x->val.number, ny));
  1847. }
  1848. return _NAN;
  1849. }
  1850. value_t *verb_take(interpreter_t *state, verb_t *self, value_t *x, value_t *y);
  1851. value_t *verb_drop(interpreter_t *state, verb_t *self, value_t *x, value_t *y);
  1852. bool is_bad_num(double v) {
  1853. return isnan(v) || v == INFINITY || v == -INFINITY;
  1854. }
  1855. value_t *verb_odometer(interpreter_t *state, verb_t *self, value_t *x) {
  1856. if (x->tag != ARRAY)
  1857. x = verb_enlist(state, NULL, x);
  1858. else if (x->val.array->length < 2)
  1859. return state->udf;
  1860. size_t p = 1;
  1861. size_t xl = x->val.array->length;
  1862. for (size_t i = 0; i < xl; i++) {
  1863. value_t *it = x->val.array->data[i];
  1864. if (it->tag != NUMBER || is_bad_num(it->val.number))
  1865. return state->udf;
  1866. p *= (size_t)(it->val.number);
  1867. }
  1868. if (p < 1)
  1869. return state->unit;
  1870. uint64_t *lims = malloc_checked_atomic(sizeof(uint64_t) * xl);
  1871. for (size_t i = 0; i < xl; i++)
  1872. lims[i] = (size_t)(((value_t *)x->val.array->data[i])->val.number);
  1873. uint64_t **z = malloc_checked(sizeof(uint64_t *) * p);
  1874. for (size_t i = 0; i < p; i++)
  1875. z[i] = malloc_checked_atomic(sizeof(uint64_t) * xl);
  1876. for (size_t i = 0; i < p - 1; i++) {
  1877. uint64_t *r = z[i];
  1878. uint64_t *s = z[i + 1];
  1879. bool carry = true;
  1880. for (size_t j = 0; j < xl; j++) {
  1881. uint64_t a = xl - 1 - j;
  1882. s[a] = r[a];
  1883. if (carry) {
  1884. s[a]++;
  1885. carry = false;
  1886. }
  1887. if (s[a] >= lims[a]) {
  1888. s[a] = 0;
  1889. carry = true;
  1890. }
  1891. }
  1892. }
  1893. GC_FREE(lims);
  1894. list_t *r = list_newk(p);
  1895. for (size_t i = 0; i < p; i++) {
  1896. list_t *rw = list_newk(xl);
  1897. for (size_t j = 0; j < xl; j++)
  1898. rw->data[j] = value_new_number(z[i][j]);
  1899. r->data[i] = value_new_array(rw);
  1900. GC_FREE(z[i]);
  1901. }
  1902. GC_FREE(z);
  1903. return value_new_array(r);
  1904. }
  1905. value_t *verb_chunks(interpreter_t *state, verb_t *self, value_t *x,
  1906. value_t *y) {
  1907. if (x->tag != NUMBER)
  1908. return state->udf;
  1909. if (y->tag != ARRAY)
  1910. y = verb_enlist(state, NULL, y);
  1911. else if (!y->val.array->data)
  1912. return y;
  1913. list_t *r = list_new();
  1914. size_t cl = fabs(x->val.number);
  1915. for (size_t i = 0; i < y->val.array->length; i += cl)
  1916. list_push(r, verb_take(state, NULL, value_new_number(cl),
  1917. verb_drop(state, NULL, value_new_number(i), y)));
  1918. return value_new_array(r);
  1919. }
  1920. value_t *verb_exp(interpreter_t *state, verb_t *self, value_t *x) {
  1921. if (x->tag == NUMBER)
  1922. return value_new_number(exp(x->val.number));
  1923. return _NAN;
  1924. }
  1925. value_t *verb_power(interpreter_t *state, verb_t *self, value_t *x,
  1926. value_t *y) {
  1927. if (x->tag == NUMBER && y->tag == NUMBER)
  1928. return value_new_number(pow(x->val.number, y->val.number));
  1929. return _NAN;
  1930. }
  1931. value_t *verb_nlog(interpreter_t *state, verb_t *self, value_t *x) {
  1932. if (x->tag == NUMBER)
  1933. return value_new_number(log(x->val.number));
  1934. return _NAN;
  1935. }
  1936. value_t *verb_log(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1937. if (x->tag == NUMBER && y->tag == NUMBER)
  1938. return value_new_number(log(y->val.number) / log(x->val.number));
  1939. return _NAN;
  1940. }
  1941. int bits_needed(uint32_t value) {
  1942. int bits = 0;
  1943. for (int bit_test = 16; bit_test > 0; bit_test >>= 1) {
  1944. if (value >> bit_test != 0) {
  1945. bits += bit_test;
  1946. value >>= bit_test;
  1947. }
  1948. }
  1949. return bits + value;
  1950. }
  1951. value_t *verb_bits(interpreter_t *state, verb_t *self, value_t *x) {
  1952. if (x->tag == NUMBER) {
  1953. int n = x->val.number;
  1954. int bk = bits_needed(n);
  1955. list_t *r = list_newk(bk);
  1956. for (int i = 0; i < bk; i++)
  1957. if ((n & (1 << i)) >> i)
  1958. r->data[i] = NUMS[1];
  1959. else
  1960. r->data[i] = NUMS[0];
  1961. return value_new_array(r);
  1962. }
  1963. return state->udf;
  1964. }
  1965. value_t *verb_reverse(interpreter_t *state, verb_t *self, value_t *x);
  1966. value_t *verb_base(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  1967. if (x->tag == NUMBER && y->tag == NUMBER) {
  1968. size_t v = fabs(y->val.number);
  1969. size_t b = fabs(x->val.number);
  1970. if (b < 2)
  1971. return state->udf;
  1972. list_t *r = list_new();
  1973. while (v > 0) {
  1974. list_push(r, value_new_number(v % b));
  1975. v /= b;
  1976. }
  1977. return verb_reverse(state, NULL, value_new_array(r));
  1978. }
  1979. return state->udf;
  1980. }
  1981. ssize_t indexOf(list_t *l, value_t *x) {
  1982. if (!l->data)
  1983. return -1;
  1984. for (size_t i = 0; i < l->length; i++)
  1985. if (value_equals(l->data[i], x))
  1986. return i;
  1987. return -1;
  1988. }
  1989. value_t *verb_group(interpreter_t *state, verb_t *self, value_t *x) {
  1990. if (x->tag != ARRAY)
  1991. x = verb_enlist(state, NULL, x);
  1992. else if (!x->val.array->data)
  1993. return x;
  1994. list_t *r = list_new();
  1995. list_t *is = list_new();
  1996. for (size_t i = 0; i < x->val.array->length; i++) {
  1997. value_t *v = x->val.array->data[i];
  1998. ssize_t n = indexOf(is, v);
  1999. if (n < 0) {
  2000. list_push(r, verb_enlist(state, NULL, value_new_number(i)));
  2001. list_push(is, v);
  2002. } else {
  2003. value_t *tmp = list_index(r, n);
  2004. list_push(tmp->val.array, value_new_number(i));
  2005. }
  2006. }
  2007. GC_FREE(is->data);
  2008. GC_FREE(is);
  2009. return value_new_array(r);
  2010. }
  2011. value_t *verb_buckets(interpreter_t *state, verb_t *self, value_t *x,
  2012. value_t *y) {
  2013. if (x->tag != ARRAY)
  2014. x = verb_enlist(state, NULL, x);
  2015. else if (!x->val.array->data)
  2016. return y;
  2017. if (y->tag != ARRAY)
  2018. y = verb_enlist(state, NULL, x);
  2019. else if (!y->val.array->data)
  2020. return y;
  2021. list_t *r = list_new();
  2022. size_t mx = 0;
  2023. for (size_t i = 0; i < x->val.array->length; i++) {
  2024. value_t *v = x->val.array->data[i];
  2025. if (v->tag != NUMBER)
  2026. break;
  2027. ssize_t j = v->val.number;
  2028. if (j >= 0 && j > mx)
  2029. mx = j;
  2030. }
  2031. for (size_t i = 0; i < mx + 1; i++)
  2032. list_push(r, list_new());
  2033. if (!r->data) {
  2034. GC_FREE(r);
  2035. return state->unit;
  2036. }
  2037. for (size_t i = 0; i < x->val.array->length; i++) {
  2038. if (i >= y->val.array->length)
  2039. break;
  2040. value_t *v = x->val.array->data[i];
  2041. if (v->tag != NUMBER)
  2042. break;
  2043. ssize_t j = v->val.number;
  2044. if (j >= 0) {
  2045. list_t *b = list_index(r, j);
  2046. if (b)
  2047. list_push(b, y->val.array->data[i]);
  2048. }
  2049. }
  2050. if (x->val.array->length < y->val.array->length) {
  2051. list_t *lb = list_new();
  2052. for (size_t i = x->val.array->length; i < y->val.array->length; i++)
  2053. list_push(lb, y->val.array->data[i]);
  2054. list_push(r, lb);
  2055. }
  2056. for (size_t i = 0; i < r->length; i++)
  2057. r->data[i] = value_new_array(r->data[i]);
  2058. return value_new_array(r);
  2059. }
  2060. value_t *verb_equals(interpreter_t *state, verb_t *self, value_t *x,
  2061. value_t *y) {
  2062. return value_equals(x, y) ? NUMS[1] : NUMS[0];
  2063. }
  2064. value_t *verb_permute(interpreter_t *state, verb_t *self, value_t *x) {
  2065. if (x->tag != ARRAY || x->val.array->length < 2)
  2066. return x;
  2067. list_t *permutation = list_copy(x->val.array);
  2068. size_t length = permutation->length;
  2069. list_t *result = list_new();
  2070. list_push(result, list_copy(permutation));
  2071. list_t *c = list_new();
  2072. for (size_t i = 0; i < length; i++) {
  2073. size_t *n = malloc_checked_atomic(sizeof(size_t));
  2074. list_push(c, n);
  2075. }
  2076. size_t k;
  2077. list_t *p;
  2078. size_t i = 0;
  2079. while (i < length) {
  2080. size_t *n = list_index(c, i);
  2081. if ((*n) < i) {
  2082. k = i % 2 && (*n);
  2083. p = list_index(permutation, i);
  2084. list_set(permutation, i, list_index(permutation, k));
  2085. list_set(permutation, k, p);
  2086. *n = (*n) + 1;
  2087. i = 1;
  2088. list_push(result, list_copy(permutation));
  2089. } else {
  2090. *n = 0;
  2091. i++;
  2092. }
  2093. }
  2094. for (size_t i = 0; i < c->length; i++)
  2095. GC_FREE(c->data[i]);
  2096. GC_FREE(c->data);
  2097. GC_FREE(c);
  2098. GC_FREE(permutation->data);
  2099. GC_FREE(permutation);
  2100. for (size_t i = 0; i < result->length; i++)
  2101. result->data[i] = value_new_array(result->data[i]);
  2102. return value_new_array(result);
  2103. }
  2104. value_t *verb_occurences(interpreter_t *state, verb_t *self, value_t *x) {
  2105. if (x->tag != ARRAY)
  2106. x = verb_enlist(state, NULL, x);
  2107. else if (!x->val.array->data)
  2108. return x;
  2109. list_t *table = list_new();
  2110. list_t *r = list_new();
  2111. for (size_t i = 0; i < x->val.array->length; i++) {
  2112. bool f = false;
  2113. value_t *it = x->val.array->data[i];
  2114. for (size_t j = 0; j < table->length; j++) {
  2115. list_t *p = table->data[j];
  2116. if (value_equals(p->data[0], it)) {
  2117. size_t *n = p->data[1];
  2118. *n = (*n) + 1;
  2119. list_push(r, value_new_number(*n));
  2120. f = true;
  2121. break;
  2122. }
  2123. }
  2124. if (!f) {
  2125. list_t *p = list_newk(2);
  2126. p->data[0] = it;
  2127. size_t *n = malloc_checked_atomic(sizeof(size_t));
  2128. p->data[1] = n;
  2129. list_push(table, p);
  2130. list_push(r, NUMS[0]);
  2131. }
  2132. }
  2133. for (size_t i = 0; i < table->length; i++) {
  2134. list_t *p = table->data[i];
  2135. GC_FREE(p->data[1]);
  2136. GC_FREE(p->data);
  2137. GC_FREE(p);
  2138. }
  2139. GC_FREE(table->data);
  2140. GC_FREE(table);
  2141. return value_new_array(r);
  2142. }
  2143. value_t *verb_mask(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2144. if (x->tag != ARRAY)
  2145. x = verb_enlist(state, NULL, x);
  2146. else if (!x->val.array->data)
  2147. return x;
  2148. if (y->tag != ARRAY)
  2149. y = verb_enlist(state, NULL, y);
  2150. list_t *r = list_new();
  2151. value_t *l = value_new_number(y->val.array->length);
  2152. size_t n = 0;
  2153. size_t k = x->val.array->length;
  2154. for (size_t i = 0; i < k; i++) {
  2155. value_t *s = verb_take(state, NULL, l,
  2156. verb_drop(state, NULL, value_new_number(i), x));
  2157. if (value_equals(s, y)) {
  2158. n++;
  2159. for (size_t j = 0; j < l->val.number; j++, i++)
  2160. list_push(r, value_new_number(n));
  2161. i--;
  2162. } else
  2163. list_push(r, NUMS[0]);
  2164. }
  2165. return value_new_array(r);
  2166. }
  2167. value_t *verb_classify(interpreter_t *state, verb_t *self, value_t *x) {
  2168. if (x->tag != ARRAY)
  2169. x = verb_enlist(state, NULL, x);
  2170. else if (!x->val.array->data)
  2171. return x;
  2172. list_t *table = list_new();
  2173. list_t *r = list_new();
  2174. for (size_t i = 0; i < x->val.array->length; i++) {
  2175. bool f = false;
  2176. value_t *it = x->val.array->data[i];
  2177. for (size_t j = 0; j < table->length; j++) {
  2178. list_t *p = table->data[j];
  2179. if (value_equals(p->data[0], it)) {
  2180. size_t *n = p->data[1];
  2181. list_push(r, value_new_number(*n));
  2182. f = true;
  2183. break;
  2184. }
  2185. }
  2186. if (!f) {
  2187. list_t *p = list_newk(2);
  2188. p->data[0] = it;
  2189. size_t *n = malloc_checked_atomic(sizeof(size_t));
  2190. *n = i++;
  2191. p->data[1] = n;
  2192. list_push(table, p);
  2193. list_push(r, value_new_number(*n));
  2194. }
  2195. }
  2196. for (size_t i = 0; i < table->length; i++) {
  2197. list_t *p = table->data[i];
  2198. GC_FREE(p->data[1]);
  2199. GC_FREE(p->data);
  2200. GC_FREE(p);
  2201. }
  2202. GC_FREE(table->data);
  2203. GC_FREE(table);
  2204. return value_new_array(r);
  2205. }
  2206. value_t *verb_unbits(interpreter_t *state, verb_t *self, value_t *x) {
  2207. if (x->tag != ARRAY)
  2208. x = verb_enlist(state, NULL, x);
  2209. int n = 0;
  2210. for (size_t i = 0; i < x->val.array->length; i++) {
  2211. if (value_is_truthy(x->val.array->data[i]))
  2212. n |= (int)1 << (int)i;
  2213. else
  2214. n &= ~((int)1 << (int)i);
  2215. }
  2216. return value_new_number(n);
  2217. }
  2218. value_t *verb_unbase(interpreter_t *state, verb_t *self, value_t *x,
  2219. value_t *y) {
  2220. if (x->tag == NUMBER) {
  2221. size_t b = fabs(x->val.number);
  2222. if (b < 2)
  2223. return state->udf;
  2224. if (y->tag != ARRAY)
  2225. y = verb_enlist(state, NULL, y);
  2226. size_t n = 0;
  2227. if (!y->val.array->data)
  2228. return state->udf;
  2229. for (size_t i = 0; i < y->val.array->length; i++) {
  2230. value_t *v = y->val.array->data[i];
  2231. if (v->tag != NUMBER)
  2232. break;
  2233. size_t k = fabs(v->val.number);
  2234. n = n * b + k;
  2235. }
  2236. return value_new_number(n);
  2237. }
  2238. return state->udf;
  2239. }
  2240. value_t *verb_not(interpreter_t *state, verb_t *self, value_t *x) {
  2241. return value_is_truthy(x) ? NUMS[0] : NUMS[1];
  2242. }
  2243. value_t *verb_not_equals(interpreter_t *state, verb_t *self, value_t *x,
  2244. value_t *y) {
  2245. return !value_equals(x, y) ? NUMS[1] : NUMS[0];
  2246. }
  2247. value_t *verb_pred(interpreter_t *state, verb_t *self, value_t *x) {
  2248. if (x->tag == NUMBER)
  2249. return value_new_number(x->val.number - 1);
  2250. else if (x->tag == CHAR)
  2251. return value_new_char(x->val._char - 1);
  2252. return _NAN;
  2253. }
  2254. value_t *verb_less(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2255. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2256. (y->tag == NUMBER || y->tag == CHAR)) {
  2257. if (get_numeric(x) < get_numeric(y))
  2258. return NUMS[1];
  2259. return NUMS[0];
  2260. }
  2261. return _NAN;
  2262. }
  2263. value_t *verb_floor(interpreter_t *state, verb_t *self, value_t *x) {
  2264. if (x->tag == NUMBER)
  2265. return value_new_number(floor(x->val.number));
  2266. return _NAN;
  2267. }
  2268. int _compare_up(const void *a, const void *b) {
  2269. value_t *x = (*(list_t **)a)->data[0];
  2270. value_t *y = (*(list_t **)b)->data[0];
  2271. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2272. (y->tag == NUMBER || y->tag == CHAR)) {
  2273. if (get_numeric(x) > get_numeric(y))
  2274. return 1;
  2275. else if (get_numeric(x) < get_numeric(y))
  2276. return -1;
  2277. return 0;
  2278. }
  2279. return 0;
  2280. }
  2281. int _compare_down(const void *a, const void *b) {
  2282. value_t *x = (*(list_t **)a)->data[0];
  2283. value_t *y = (*(list_t **)b)->data[0];
  2284. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2285. (y->tag == NUMBER || y->tag == CHAR)) {
  2286. if (get_numeric(x) > get_numeric(y))
  2287. return -1;
  2288. else if (get_numeric(x) < get_numeric(y))
  2289. return 1;
  2290. return 0;
  2291. }
  2292. return 0;
  2293. }
  2294. value_t *_grade(value_t *x, bool down) {
  2295. if (x->tag != ARRAY || x->val.array->length < 2)
  2296. return x;
  2297. list_t *ps = list_newk(x->val.array->length);
  2298. for (size_t i = 0; i < x->val.array->length; i++) {
  2299. list_t *p = list_newk(2);
  2300. p->data[0] = x->val.array->data[i];
  2301. p->data[1] = value_new_number(i);
  2302. ps->data[i] = p;
  2303. }
  2304. qsort(ps->data, ps->length, sizeof(void *),
  2305. down ? _compare_down : _compare_up);
  2306. for (size_t i = 0; i < ps->length; i++) {
  2307. list_t *p = ps->data[i];
  2308. ps->data[i] = p->data[1];
  2309. GC_FREE(p->data);
  2310. GC_FREE(p);
  2311. }
  2312. return value_new_array(ps);
  2313. }
  2314. value_t *verb_gradedown(interpreter_t *state, verb_t *self, value_t *x) {
  2315. return _grade(x, true);
  2316. }
  2317. value_t *verb_nudge_left(interpreter_t *state, verb_t *self, value_t *x,
  2318. value_t *y) {
  2319. if (y->tag != ARRAY)
  2320. return verb_enlist(state, NULL, x);
  2321. else if (!y->val.array->data)
  2322. return y;
  2323. else if (y->val.array->length < 2)
  2324. return verb_enlist(state, NULL, x);
  2325. list_t *r = list_new();
  2326. for (size_t i = 1; i < y->val.array->length; i++)
  2327. list_push(r, y->val.array->data[i]);
  2328. list_push(r, x);
  2329. return value_new_array(r);
  2330. }
  2331. value_t *verb_lesseq(interpreter_t *state, verb_t *self, value_t *x,
  2332. value_t *y) {
  2333. if (value_equals(x, y))
  2334. return NUMS[1];
  2335. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2336. (y->tag == NUMBER || y->tag == CHAR)) {
  2337. if (get_numeric(x) < get_numeric(y))
  2338. return NUMS[1];
  2339. return NUMS[0];
  2340. }
  2341. return _NAN;
  2342. }
  2343. value_t *verb_succ(interpreter_t *state, verb_t *self, value_t *x) {
  2344. if (x->tag == NUMBER)
  2345. return value_new_number(x->val.number + 1);
  2346. else if (x->tag == CHAR)
  2347. return value_new_char(x->val._char + 1);
  2348. return _NAN;
  2349. }
  2350. value_t *verb_ceil(interpreter_t *state, verb_t *self, value_t *x) {
  2351. if (x->tag == NUMBER)
  2352. return value_new_number(ceil(x->val.number));
  2353. return _NAN;
  2354. }
  2355. value_t *verb_greater(interpreter_t *state, verb_t *self, value_t *x,
  2356. value_t *y) {
  2357. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2358. (y->tag == NUMBER || y->tag == CHAR)) {
  2359. if (get_numeric(x) > get_numeric(y))
  2360. return NUMS[1];
  2361. return NUMS[0];
  2362. }
  2363. return _NAN;
  2364. }
  2365. value_t *verb_greatereq(interpreter_t *state, verb_t *self, value_t *x,
  2366. value_t *y) {
  2367. if (value_equals(x, y))
  2368. return NUMS[1];
  2369. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2370. (y->tag == NUMBER || y->tag == CHAR)) {
  2371. if (get_numeric(x) > get_numeric(y))
  2372. return NUMS[1];
  2373. return NUMS[0];
  2374. }
  2375. return _NAN;
  2376. }
  2377. value_t *verb_gradeup(interpreter_t *state, verb_t *self, value_t *x) {
  2378. return _grade(x, false);
  2379. }
  2380. value_t *verb_nudge_right(interpreter_t *state, verb_t *self, value_t *x,
  2381. value_t *y) {
  2382. if (y->tag != ARRAY)
  2383. return verb_enlist(state, NULL, x);
  2384. else if (!y->val.array->data)
  2385. return y;
  2386. else if (y->val.array->length < 2)
  2387. return verb_enlist(state, NULL, x);
  2388. list_t *r = list_new();
  2389. list_push(r, x);
  2390. for (size_t i = 0; i < y->val.array->length - 1; i++)
  2391. list_push(r, y->val.array->data[i]);
  2392. return value_new_array(r);
  2393. }
  2394. value_t *verb_enlist(interpreter_t *state, verb_t *self, value_t *x) {
  2395. list_t *l = list_new();
  2396. list_push(l, x);
  2397. return value_new_array(l);
  2398. }
  2399. value_t *verb_enpair(interpreter_t *state, verb_t *self, value_t *x,
  2400. value_t *y);
  2401. value_t *verb_join(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2402. list_t *l;
  2403. if (x->tag == ARRAY && y->tag == ARRAY) {
  2404. if (!x->val.array->data && !y->val.array->data)
  2405. return state->unit;
  2406. else if (!x->val.array->data)
  2407. return y;
  2408. else if (!y->val.array->data)
  2409. return x;
  2410. l = list_newk(x->val.array->length + y->val.array->length);
  2411. size_t lp = 0;
  2412. for (size_t i = 0; i < x->val.array->length; i++)
  2413. l->data[lp++] = x->val.array->data[i];
  2414. for (size_t i = 0; i < y->val.array->length; i++)
  2415. l->data[lp++] = y->val.array->data[i];
  2416. } else if (x->tag == ARRAY && y->tag != ARRAY) {
  2417. if (!x->val.array->data)
  2418. return verb_enlist(state, NULL, y);
  2419. l = list_newk(x->val.array->length + 1);
  2420. size_t lp = 0;
  2421. for (size_t i = 0; i < x->val.array->length; i++)
  2422. l->data[lp++] = x->val.array->data[i];
  2423. l->data[lp++] = y;
  2424. } else if (x->tag != ARRAY && y->tag == ARRAY) {
  2425. if (!y->val.array->data)
  2426. return verb_enlist(state, NULL, x);
  2427. l = list_newk(y->val.array->length + 1);
  2428. size_t lp = 0;
  2429. l->data[lp++] = x;
  2430. for (size_t i = 0; i < y->val.array->length; i++)
  2431. l->data[lp++] = y->val.array->data[i];
  2432. } else
  2433. return verb_enpair(state, NULL, x, y);
  2434. return value_new_array(l);
  2435. }
  2436. value_t *verb_enpair(interpreter_t *state, verb_t *self, value_t *x,
  2437. value_t *y) {
  2438. list_t *l = list_newk(2);
  2439. l->data[0] = x;
  2440. l->data[1] = y;
  2441. return value_new_array(l);
  2442. }
  2443. value_t *verb_selfref1(interpreter_t *state, verb_t *self, value_t *x) {
  2444. verb_t *v;
  2445. if (state->args->data)
  2446. v = list_index(list_index(state->args, -1), -1);
  2447. else if (state->selfrefs->data)
  2448. v = list_index(state->selfrefs, -1);
  2449. else
  2450. return state->udf;
  2451. return each_rank(state, v, x, 0, v->rank[0]);
  2452. }
  2453. value_t *verb_selfref2(interpreter_t *state, verb_t *self, value_t *x,
  2454. value_t *y) {
  2455. verb_t *v;
  2456. if (state->args->data)
  2457. v = list_index(list_index(state->args, -1), -1);
  2458. else if (state->selfrefs->data)
  2459. v = list_index(state->selfrefs, -1);
  2460. else
  2461. return state->udf;
  2462. return together(state, v, x, y, 0, 0, v->rank[1], v->rank[2]);
  2463. }
  2464. value_t *verb_take(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2465. if (x->tag == NUMBER) {
  2466. if (y->tag != ARRAY) {
  2467. if (x->val.number == 0)
  2468. return state->unit;
  2469. else
  2470. return x;
  2471. }
  2472. if (x->val.number == 0 || !y->val.array->data)
  2473. return state->unit;
  2474. bool rev = x->val.number < 0;
  2475. size_t k = (size_t)fabs(x->val.number);
  2476. list_t *r = list_newk(y->val.array->length < k ? y->val.array->length : k);
  2477. size_t p = 0;
  2478. if (rev)
  2479. for (ssize_t i = k; i > 0; i--) {
  2480. value_t *v = list_index(y->val.array, -i);
  2481. if (!v)
  2482. continue;
  2483. r->data[p++] = v;
  2484. }
  2485. else
  2486. for (size_t i = 0; i < y->val.array->length && k; i++, k--)
  2487. r->data[p++] = y->val.array->data[i];
  2488. return value_new_array(r);
  2489. }
  2490. return state->udf;
  2491. }
  2492. value_t *verb_where(interpreter_t *state, verb_t *self, value_t *x) {
  2493. if (x->tag != ARRAY)
  2494. x = verb_enlist(state, NULL, x);
  2495. else if (!x->val.array->data)
  2496. return x;
  2497. list_t *r = list_new();
  2498. for (size_t i = 0; i < x->val.array->length; i++) {
  2499. value_t *a = x->val.array->data[i];
  2500. if (a->tag != NUMBER)
  2501. break;
  2502. size_t k = fabs(a->val.number);
  2503. for (size_t j = 0; j < k; j++)
  2504. list_push(r, value_new_number(i));
  2505. }
  2506. return value_new_array(r);
  2507. }
  2508. value_t *verb_copy(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2509. if (x->tag != ARRAY)
  2510. x = verb_enlist(state, NULL, x);
  2511. if (y->tag != ARRAY)
  2512. y = verb_enlist(state, NULL, y);
  2513. list_t *tx = x->val.array;
  2514. list_t *ty = y->val.array;
  2515. if (!tx->data || !ty->data)
  2516. return state->unit;
  2517. list_t *r = list_new();
  2518. for (size_t i = 0; i < tx->length; i++) {
  2519. value_t *a = tx->data[i];
  2520. value_t *b = ty->data[i >= ty->length ? ty->length - 1 : i];
  2521. if (b->tag != NUMBER)
  2522. break;
  2523. size_t k = fabs(b->val.number);
  2524. for (size_t i = 0; i < k; i++)
  2525. list_push(r, a);
  2526. }
  2527. return value_new_array(r);
  2528. }
  2529. value_t *verb_nub(interpreter_t *state, verb_t *self, value_t *x) {
  2530. if (x->tag != ARRAY || !x->val.array->data)
  2531. return x;
  2532. list_t *n = list_newk(x->val.array->length);
  2533. list_t *r = list_new();
  2534. for (size_t i = 0; i < x->val.array->length; i++) {
  2535. bool u = true;
  2536. for (size_t j = 0; j < r->length; j++)
  2537. if (value_equals(x->val.array->data[i], r->data[j])) {
  2538. u = false;
  2539. break;
  2540. }
  2541. if (u)
  2542. list_push(r, x->val.array->data[i]);
  2543. n->data[i] = u ? NUMS[1] : NUMS[0];
  2544. }
  2545. GC_FREE(r->data);
  2546. GC_FREE(r);
  2547. return value_new_array(n);
  2548. }
  2549. value_t *verb_drop(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2550. if (x->tag == NUMBER) {
  2551. if (y->tag != ARRAY) {
  2552. if (x->val.number == 0)
  2553. return y;
  2554. else
  2555. return state->unit;
  2556. }
  2557. if (x->val.number == 0)
  2558. return y;
  2559. if (!y->val.array->data)
  2560. return state->unit;
  2561. bool rev = x->val.number < 0;
  2562. size_t k = (size_t)fabs(x->val.number);
  2563. if (k >= y->val.array->length)
  2564. return state->unit;
  2565. if (rev) {
  2566. size_t l = y->val.array->length;
  2567. if (k >= l)
  2568. return state->unit;
  2569. return verb_take(state, NULL, value_new_number(l - k), y);
  2570. }
  2571. list_t *r = list_newk(y->val.array->length - k);
  2572. size_t rp = 0;
  2573. for (size_t i = k; i < y->val.array->length; i++)
  2574. r->data[rp++] = y->val.array->data[i];
  2575. return value_new_array(r);
  2576. }
  2577. return state->udf;
  2578. }
  2579. value_t *verb_unique(interpreter_t *state, verb_t *self, value_t *x) {
  2580. if (x->tag != ARRAY || !x->val.array->data)
  2581. return x;
  2582. list_t *r = list_new();
  2583. for (size_t i = 0; i < x->val.array->length; i++) {
  2584. bool u = true;
  2585. for (size_t j = 0; j < r->length; j++)
  2586. if (value_equals(x->val.array->data[i], r->data[j])) {
  2587. u = false;
  2588. break;
  2589. }
  2590. if (u)
  2591. list_push(r, x->val.array->data[i]);
  2592. }
  2593. return value_new_array(r);
  2594. }
  2595. value_t *verb_find(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2596. if (y->tag != ARRAY)
  2597. y = verb_enlist(state, self, y);
  2598. else if (!y->val.array->data)
  2599. return state->unit;
  2600. list_t *r = list_new();
  2601. for (size_t i = 0; i < y->val.array->length; i++)
  2602. if (value_equals(y->val.array->data[i], x))
  2603. list_push(r, value_new_number(i));
  2604. return value_new_array(r);
  2605. }
  2606. value_t *verb_indexof(interpreter_t *state, verb_t *self, value_t *x,
  2607. value_t *y) {
  2608. if (y->tag != ARRAY)
  2609. y = verb_enlist(state, self, y);
  2610. else if (!y->val.array->data)
  2611. return state->unit;
  2612. ssize_t n = indexOf(y->val.array, x);
  2613. if (n < 0)
  2614. n = y->val.array->length;
  2615. return value_new_number(n);
  2616. }
  2617. value_t *verb_count(interpreter_t *state, verb_t *self, value_t *x) {
  2618. if (x->tag != ARRAY)
  2619. return NUMS[1];
  2620. return value_new_number(x->val.array->length);
  2621. }
  2622. void flatten(value_t *x, list_t *r) {
  2623. if (x->tag == ARRAY)
  2624. for (size_t i = 0; i < x->val.array->length; i++)
  2625. flatten(x->val.array->data[i], r);
  2626. else
  2627. list_push(r, x);
  2628. }
  2629. value_t *verb_flatten(interpreter_t *state, verb_t *self, value_t *x) {
  2630. if (x->tag != ARRAY || !x->val.array->data)
  2631. return x;
  2632. list_t *r = list_new();
  2633. flatten(x, r);
  2634. return value_new_array(r);
  2635. }
  2636. value_t *verb_minand(interpreter_t *state, verb_t *self, value_t *x,
  2637. value_t *y) {
  2638. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2639. (y->tag == NUMBER || y->tag == CHAR)) {
  2640. if (get_numeric(x) < get_numeric(y))
  2641. return x;
  2642. return y;
  2643. }
  2644. return _NAN;
  2645. }
  2646. value_t *verb_reverse(interpreter_t *state, verb_t *self, value_t *x) {
  2647. if (x->tag != ARRAY)
  2648. return x;
  2649. if (x->val.array->length < 2)
  2650. return x;
  2651. list_t *r = list_newk(x->val.array->length);
  2652. size_t rp = 0;
  2653. for (ssize_t i = x->val.array->length - 1; i >= 0; i--)
  2654. r->data[rp++] = x->val.array->data[i];
  2655. return value_new_array(r);
  2656. }
  2657. value_t *verb_maxor(interpreter_t *state, verb_t *self, value_t *x,
  2658. value_t *y) {
  2659. if ((x->tag == NUMBER || x->tag == CHAR) &&
  2660. (y->tag == NUMBER || y->tag == CHAR)) {
  2661. if (get_numeric(x) > get_numeric(y))
  2662. return x;
  2663. return y;
  2664. }
  2665. return _NAN;
  2666. }
  2667. value_t *verb_rotate(interpreter_t *state, verb_t *self, value_t *x,
  2668. value_t *y) {
  2669. if (y->tag != ARRAY || y->val.array->length < 2)
  2670. return x;
  2671. if (x->tag != NUMBER)
  2672. return state->udf;
  2673. bool rev = x->val.number < 0;
  2674. size_t k = fabs(x->val.number);
  2675. list_t *r = list_copy(y->val.array);
  2676. for (size_t i = 0; i < k; i++) {
  2677. value_t *v;
  2678. if (rev) {
  2679. v = r->data[0];
  2680. for (size_t j = 0; j < r->length - 1; j++)
  2681. r->data[j] = r->data[j + 1];
  2682. r->data[r->length - 1] = v;
  2683. } else {
  2684. v = r->data[r->length - 1];
  2685. for (size_t j = r->length - 1; j > 0; j--)
  2686. r->data[j] = r->data[j - 1];
  2687. r->data[0] = v;
  2688. }
  2689. }
  2690. return value_new_array(r);
  2691. }
  2692. value_t *verb_windows(interpreter_t *state, verb_t *self, value_t *x,
  2693. value_t *y) {
  2694. if (y->tag != ARRAY)
  2695. y = verb_enlist(state, NULL, y);
  2696. else if (!y->val.array->data)
  2697. return y;
  2698. size_t k = fabs(x->val.number);
  2699. size_t l = y->val.array->length;
  2700. list_t *r = list_new();
  2701. for (size_t i = 0; i < l; i++) {
  2702. if (i + k > l)
  2703. break;
  2704. list_push(r, verb_take(state, NULL, value_new_number(k),
  2705. verb_drop(state, NULL, value_new_number(i), y)));
  2706. }
  2707. return value_new_array(r);
  2708. }
  2709. size_t depthOf(value_t *x, size_t d) {
  2710. if (x->tag == ARRAY) {
  2711. if (!x->val.array->data)
  2712. return 0;
  2713. for (size_t i = 0; i < x->val.array->length; i++) {
  2714. size_t d2 = depthOf(x->val.array->data[i], d + 1);
  2715. if (d2 > d)
  2716. d = d2;
  2717. }
  2718. return d;
  2719. }
  2720. return 0;
  2721. }
  2722. value_t *verb_depth(interpreter_t *state, verb_t *self, value_t *x) {
  2723. return value_new_number(depthOf(x, 1));
  2724. }
  2725. value_t *verb_round(interpreter_t *state, verb_t *self, value_t *x) {
  2726. if (x->tag == NUMBER)
  2727. return value_new_number(round(x->val.number));
  2728. return _NAN;
  2729. }
  2730. value_t *verb_abs(interpreter_t *state, verb_t *self, value_t *x) {
  2731. if (x->tag == NUMBER)
  2732. return value_new_number(fabs(x->val.number));
  2733. return _NAN;
  2734. }
  2735. value_t *verb_tail(interpreter_t *state, verb_t *self, value_t *x);
  2736. value_t *verb_at(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2737. if (y->tag != NUMBER)
  2738. return state->udf;
  2739. if (x->tag != ARRAY) {
  2740. if (y->val.number > -1 && y->val.number < 1)
  2741. return x;
  2742. else
  2743. return state->udf;
  2744. }
  2745. if (!x->val.array->data)
  2746. return state->nil;
  2747. value_t *v = list_index(x->val.array, (ssize_t)y->val.number);
  2748. if (!v)
  2749. return state->udf;
  2750. return v;
  2751. }
  2752. value_t *verb_member(interpreter_t *state, verb_t *self, value_t *x,
  2753. value_t *y) {
  2754. if (y->tag != ARRAY)
  2755. y = verb_enlist(state, self, y);
  2756. else if (!y->val.array->data)
  2757. return NUMS[0];
  2758. for (size_t i = 0; i < y->val.array->length; i++)
  2759. if (value_equals(y->val.array->data[i], x))
  2760. return NUMS[1];
  2761. return NUMS[0];
  2762. }
  2763. value_t *verb_shuffle(interpreter_t *state, verb_t *self, value_t *x) {
  2764. if (x->tag != ARRAY)
  2765. x = verb_enlist(state, self, x);
  2766. else if (!x->val.array->data)
  2767. return x;
  2768. list_t *r = list_copy(x->val.array);
  2769. for (size_t i = 0; i < r->length; i++) {
  2770. size_t j = rand() % r->length;
  2771. value_t *tmp = r->data[i];
  2772. r->data[i] = r->data[j];
  2773. r->data[j] = tmp;
  2774. }
  2775. return value_new_array(r);
  2776. }
  2777. value_t *verb_head(interpreter_t *state, verb_t *self, value_t *x) {
  2778. return verb_take(state, NULL, NUMS[2], x);
  2779. }
  2780. value_t *verb_bin(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2781. if (x->tag != ARRAY)
  2782. x = verb_enlist(state, self, x);
  2783. else if (!x->val.array->data)
  2784. return x;
  2785. if (y->tag != ARRAY)
  2786. y = verb_enlist(state, self, x);
  2787. else if (!y->val.array->data)
  2788. return y;
  2789. size_t xl = x->val.array->length;
  2790. list_t *bins = list_new();
  2791. for (size_t i = 0; i < xl; i++) {
  2792. double s;
  2793. double e;
  2794. value_t *vs = x->val.array->data[i];
  2795. if (vs->tag == NUMBER)
  2796. s = vs->val.number;
  2797. else if (vs->tag == CHAR)
  2798. s = vs->val._char;
  2799. else
  2800. return state->udf;
  2801. value_t *ve =
  2802. i == xl - 1 ? value_new_number(s + 1) : x->val.array->data[i + 1];
  2803. if (ve->tag == NUMBER)
  2804. e = fabs(ve->val.number);
  2805. else if (ve->tag == CHAR)
  2806. e = ve->val._char;
  2807. else
  2808. return state->udf;
  2809. if (bins->data) {
  2810. list_t *pp = list_index(bins, -1);
  2811. double *pe = pp->data[0];
  2812. if (s <= (*pe))
  2813. return state->udf;
  2814. }
  2815. double *sn = malloc_checked(sizeof(double));
  2816. *sn = s;
  2817. double *en = malloc_checked(sizeof(double));
  2818. *en = e;
  2819. list_t *p = list_new();
  2820. list_push(p, sn);
  2821. list_push(p, en);
  2822. list_push(bins, p);
  2823. }
  2824. size_t bl = bins->length;
  2825. list_t *r = list_new();
  2826. size_t yl = y->val.array->length;
  2827. for (size_t i = 0; i < yl; i++) {
  2828. value_t *it = y->val.array->data[i];
  2829. double itv;
  2830. if (it->tag == NUMBER)
  2831. itv = it->val.number;
  2832. else if (it->tag == CHAR)
  2833. itv = it->val._char;
  2834. else
  2835. return state->udf;
  2836. list_t *b = bins->data[0];
  2837. double *s = b->data[0];
  2838. if (itv < (*s)) {
  2839. list_push(r, NNUMS[0]);
  2840. continue;
  2841. }
  2842. b = list_index(bins, -1);
  2843. s = b->data[1];
  2844. if (itv >= (*s)) {
  2845. list_push(r, value_new_number(bl - 1));
  2846. continue;
  2847. }
  2848. double v = NAN;
  2849. for (size_t j = 0; j < bl; j++) {
  2850. b = bins->data[j];
  2851. double *s = b->data[0];
  2852. double *e = b->data[1];
  2853. if (itv >= (*s) && itv < (*e)) {
  2854. v = j;
  2855. break;
  2856. }
  2857. }
  2858. if (!isnan(v))
  2859. list_push(r, value_new_number(v));
  2860. }
  2861. for (size_t j = 0; j < bl; j++) {
  2862. list_t *b = bins->data[j];
  2863. GC_FREE(b->data[0]);
  2864. GC_FREE(b->data[1]);
  2865. GC_FREE(b->data);
  2866. GC_FREE(b);
  2867. }
  2868. GC_FREE(bins->data);
  2869. GC_FREE(bins);
  2870. return value_new_array(r);
  2871. }
  2872. value_t *verb_tail(interpreter_t *state, verb_t *self, value_t *x) {
  2873. if (x->tag != ARRAY)
  2874. return x;
  2875. if (!x->val.array->data)
  2876. return state->udf;
  2877. return list_index(x->val.array, -1);
  2878. }
  2879. value_t *verb_cut(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2880. if (x->tag != ARRAY)
  2881. x = verb_enlist(state, self, x);
  2882. else if (!x->val.array->data)
  2883. return x;
  2884. if (y->tag != ARRAY)
  2885. y = verb_enlist(state, self, x);
  2886. else if (!y->val.array->data)
  2887. return x;
  2888. if (x->val.array->length != 2)
  2889. return state->udf;
  2890. value_t *vs = x->val.array->data[0];
  2891. value_t *ve = x->val.array->data[1];
  2892. if (vs->tag != NUMBER || ve->tag != NUMBER)
  2893. return state->udf;
  2894. size_t s = fabs(vs->val.number);
  2895. size_t e = fabs(ve->val.number);
  2896. list_t *r = list_new();
  2897. size_t l = y->val.array->length;
  2898. list_t *pa = list_new();
  2899. for (size_t i = s; i < e && i < l; i++) {
  2900. value_t *v = list_index(y->val.array, i);
  2901. if (!v)
  2902. break;
  2903. list_push(pa, v);
  2904. }
  2905. list_t *pb = list_new();
  2906. for (size_t i = e; i < l; i++) {
  2907. value_t *v = list_index(y->val.array, i);
  2908. if (!v)
  2909. break;
  2910. list_push(pb, v);
  2911. }
  2912. list_push(r, value_new_array(pa));
  2913. list_push(r, value_new_array(pb));
  2914. return value_new_array(r);
  2915. }
  2916. value_t *verb_prefixes(interpreter_t *state, verb_t *self, value_t *x) {
  2917. if (x->tag != ARRAY)
  2918. x = verb_enlist(state, NULL, x);
  2919. list_t *r = list_new();
  2920. for (size_t i = 0; i < x->val.array->length; i++)
  2921. list_push(r, verb_take(state, NULL, value_new_number(i), x));
  2922. list_push(r, x);
  2923. return value_new_array(r);
  2924. }
  2925. value_t *verb_behead(interpreter_t *state, verb_t *self, value_t *x) {
  2926. return verb_drop(state, NULL, NUMS[1], x);
  2927. }
  2928. value_t *verb_curtail(interpreter_t *state, verb_t *self, value_t *x) {
  2929. return verb_drop(state, NULL, NNUMS[0], x);
  2930. }
  2931. value_t *verb_suffixes(interpreter_t *state, verb_t *self, value_t *x) {
  2932. if (x->tag != ARRAY)
  2933. x = verb_enlist(state, NULL, x);
  2934. list_t *r = list_new();
  2935. for (size_t i = 0; i < x->val.array->length; i++)
  2936. list_push(r, verb_drop(state, NULL, value_new_number(i), x));
  2937. list_push(r, state->unit);
  2938. return value_new_array(r);
  2939. }
  2940. value_t *verb_left(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  2941. return x;
  2942. }
  2943. value_t *verb_same(interpreter_t *state, verb_t *self, value_t *x) { return x; }
  2944. value_t *verb_right(interpreter_t *state, verb_t *self, value_t *x,
  2945. value_t *y) {
  2946. return y;
  2947. }
  2948. value_t *verb_symbol(interpreter_t *state, verb_t *self, value_t *x) {
  2949. char *s = value_str(x);
  2950. return value_new_symbol(s);
  2951. }
  2952. value_t *verb_apply1(interpreter_t *state, verb_t *self, value_t *x,
  2953. value_t *y) {
  2954. return apply_monad(state, x, y);
  2955. }
  2956. value_t *verb_apply2(interpreter_t *state, verb_t *self, value_t *x,
  2957. value_t *y) {
  2958. if (y->tag != ARRAY || y->val.array->length < 2)
  2959. return state->udf;
  2960. return apply_dyad(state, x, y->val.array->data[0], y->val.array->data[1]);
  2961. }
  2962. value_t *verb_shape(interpreter_t *state, verb_t *self, value_t *x) {
  2963. if (x->tag != ARRAY || !x->val.array->data)
  2964. return state->unit;
  2965. if (!is_arrays_array(x->val.array))
  2966. return verb_enlist(state, NULL, verb_count(state, NULL, x));
  2967. if (x->val.array->length < 2)
  2968. return verb_enlist(state, NULL,
  2969. verb_shape(state, NULL, x->val.array->data[0]));
  2970. return verb_enpair(state, NULL, verb_count(state, NULL, x),
  2971. verb_count(state, NULL, x->val.array->data[0]));
  2972. }
  2973. value_t *verb_reshape(interpreter_t *state, verb_t *self, value_t *x,
  2974. value_t *y) {
  2975. if (y->tag != ARRAY)
  2976. y = verb_enlist(state, NULL, y);
  2977. else if (!y->val.array->data)
  2978. return y;
  2979. if (x->tag != ARRAY)
  2980. x = verb_enlist(state, NULL, x);
  2981. else if (!x->val.array->data)
  2982. return state->unit;
  2983. list_t *r;
  2984. if (x->val.array->length < 2) {
  2985. value_t *a = x->val.array->data[0];
  2986. if (a->tag != NUMBER)
  2987. return state->udf;
  2988. size_t k = fabs(a->val.number);
  2989. list_t *t = list_new();
  2990. flatten(y, t);
  2991. r = list_newk(k);
  2992. for (size_t i = 0; i < k; i++)
  2993. r->data[i] = t->data[i % t->length];
  2994. } else if (x->val.array->length > 1) {
  2995. value_t *a = x->val.array->data[0];
  2996. if (a->tag != NUMBER)
  2997. return state->udf;
  2998. value_t *b = x->val.array->data[1];
  2999. if (b->tag != NUMBER)
  3000. return state->udf;
  3001. size_t k = fabs(a->val.number);
  3002. size_t l = fabs(b->val.number);
  3003. y = verb_reshape(state, self,
  3004. verb_enlist(state, NULL, value_new_number(k * l)), y);
  3005. r = list_new();
  3006. size_t yp = 0;
  3007. while (k--) {
  3008. list_t *rw = list_new();
  3009. for (size_t i = 0; i < l; i++)
  3010. list_push(rw, y->val.array->data[yp++]);
  3011. list_push(r, value_new_array(rw));
  3012. }
  3013. } else
  3014. return state->udf;
  3015. return value_new_array(r);
  3016. }
  3017. value_t *verb_repr(interpreter_t *state, verb_t *self, value_t *x) {
  3018. char *s = value_show(x);
  3019. list_t *r = list_new();
  3020. for (size_t i = 0; i < strlen(s); i++)
  3021. list_push(r, value_new_char(s[i]));
  3022. GC_FREE(s);
  3023. return value_new_array(r);
  3024. }
  3025. char *format(char *template, list_t *replaces) {
  3026. buffer_t *text = buffer_new();
  3027. bool skip = false;
  3028. size_t ri = 0;
  3029. size_t tl = strlen(template);
  3030. size_t rl = replaces->length;
  3031. for (size_t i = 0; i < tl; i++) {
  3032. char c = template[i];
  3033. if (skip) {
  3034. buffer_append(text, c);
  3035. skip = false;
  3036. continue;
  3037. }
  3038. if (c == '_') {
  3039. char *s = value_show(list_index(replaces, ri));
  3040. buffer_append_str(text, s);
  3041. GC_FREE(s);
  3042. if (ri < rl - 1)
  3043. ri++;
  3044. continue;
  3045. } else if (c == '{') {
  3046. size_t bi = i;
  3047. buffer_t *n = buffer_new();
  3048. i++;
  3049. while (i < tl && template[i] != '}')
  3050. buffer_append(n, template[i++]);
  3051. if (i >= tl || template[i] != '}') {
  3052. GC_FREE(buffer_read(n));
  3053. buffer_append(text, '{');
  3054. i = bi;
  3055. continue;
  3056. }
  3057. char *s = buffer_read(n);
  3058. ssize_t ind = atoi(s);
  3059. GC_FREE(s);
  3060. value_t *v = list_index(replaces, ind);
  3061. if (!v)
  3062. continue;
  3063. s = value_show(v);
  3064. buffer_append_str(text, s);
  3065. GC_FREE(s);
  3066. continue;
  3067. } else if (c == '~') {
  3068. skip = true;
  3069. continue;
  3070. }
  3071. buffer_append(text, c);
  3072. }
  3073. return buffer_read(text);
  3074. }
  3075. value_t *verb_format(interpreter_t *state, verb_t *self, value_t *x,
  3076. value_t *y) {
  3077. if (y->tag != ARRAY)
  3078. y = verb_enlist(state, NULL, x);
  3079. else if (!y->val.array->data)
  3080. return y;
  3081. char *fmt = value_show(x);
  3082. char *s = format(fmt, y->val.array);
  3083. GC_FREE(fmt);
  3084. size_t z = strlen(s);
  3085. list_t *r = list_newk(z);
  3086. for (size_t i = 0; i < z; i++)
  3087. r->data[i] = CHARS[(int)s[i]];
  3088. return value_new_array(r);
  3089. }
  3090. value_t *verb_insert(interpreter_t *state, verb_t *self, value_t *x,
  3091. value_t *y) {
  3092. if (y->tag != ARRAY)
  3093. y = verb_enlist(state, NULL, y);
  3094. else if (!y->val.array->data)
  3095. return y;
  3096. list_t *r = list_newk(y->val.array->length * 2 - 1);
  3097. size_t rp = 0;
  3098. for (size_t i = 0; i < y->val.array->length; i++) {
  3099. r->data[rp++] = y->val.array->data[i];
  3100. if (i != y->val.array->length - 1)
  3101. r->data[rp++] = x;
  3102. }
  3103. return value_new_array(r);
  3104. }
  3105. uint64_t fibonacci(uint64_t n) {
  3106. uint64_t a = 0;
  3107. uint64_t b = 1;
  3108. while (n-- > 1) {
  3109. uint64_t t = a;
  3110. a = b;
  3111. b += t;
  3112. }
  3113. return b;
  3114. }
  3115. value_t *verb_fibonacci(interpreter_t *state, verb_t *self, value_t *x) {
  3116. if (x->tag == NUMBER)
  3117. return value_new_number(fibonacci((uint64_t)fabs(x->val.number)));
  3118. return _NAN;
  3119. }
  3120. value_t *verb_iota(interpreter_t *state, verb_t *self, value_t *x) {
  3121. if (value_equals(x, NUMS[1]))
  3122. return verb_enlist(state, NULL, NUMS[1]);
  3123. else if (value_equals(x, NUMS[0]))
  3124. return state->unit;
  3125. return verb_range(state, self, NUMS[1], x);
  3126. }
  3127. value_t *verb_range(interpreter_t *state, verb_t *self, value_t *x,
  3128. value_t *y) {
  3129. if ((x->tag == NUMBER || x->tag == CHAR) &&
  3130. (y->tag == NUMBER || y->tag == CHAR)) {
  3131. if (x->tag == NUMBER && is_bad_num(x->val.number))
  3132. return state->udf;
  3133. if (y->tag == NUMBER && is_bad_num(y->val.number))
  3134. return state->udf;
  3135. ssize_t s = get_numeric(x);
  3136. ssize_t e = get_numeric(y);
  3137. if (s == e)
  3138. return verb_enlist(state, NULL, x);
  3139. size_t p = 0;
  3140. list_t *r = list_newk((s > e ? s - e : e - s) + 1);
  3141. if (s > e)
  3142. for (ssize_t i = s; i >= e; i--) {
  3143. if (x->tag == CHAR || y->tag == CHAR)
  3144. r->data[p++] = CHARS[i];
  3145. else
  3146. r->data[p++] = value_new_number(i);
  3147. }
  3148. else
  3149. for (ssize_t i = s; i <= e; i++) {
  3150. if (x->tag == CHAR || y->tag == CHAR)
  3151. r->data[p++] = CHARS[i];
  3152. else
  3153. r->data[p++] = value_new_number(i);
  3154. }
  3155. return value_new_array(r);
  3156. }
  3157. return _NAN;
  3158. }
  3159. value_t *verb_deal(interpreter_t *state, verb_t *self, value_t *x) {
  3160. if (x->tag != ARRAY)
  3161. return x;
  3162. if (!x->val.array->data)
  3163. return state->udf;
  3164. return x->val.array->data[rand() % x->val.array->length];
  3165. }
  3166. value_t *verb_roll(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  3167. if (x->tag == NUMBER && y->tag == NUMBER) {
  3168. size_t k = fabs(x->val.number);
  3169. size_t d = fabs(y->val.number);
  3170. list_t *r = list_newk(k);
  3171. for (size_t i = 0; i < k; i++)
  3172. r->data[i] = value_new_number(rand() % d);
  3173. return value_new_array(r);
  3174. }
  3175. return state->udf;
  3176. }
  3177. value_t *verb_type(interpreter_t *state, verb_t *self, value_t *x) {
  3178. return NUMS[x->tag];
  3179. }
  3180. value_t *verb_cast(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  3181. if (x->tag == NUMBER) {
  3182. int t = fabs(x->val.number);
  3183. if (y->tag == t)
  3184. return y;
  3185. switch (t) {
  3186. case ARRAY:
  3187. if (y->tag == SYMBOL) {
  3188. char *s = y->val.symbol;
  3189. size_t z = strlen(s);
  3190. list_t *r = list_newk(z);
  3191. for (size_t i = 0; i < z; i++)
  3192. r->data[i] = CHARS[(int)s[i]];
  3193. return value_new_array(r);
  3194. }
  3195. break;
  3196. case NUMBER:
  3197. if (y->tag == CHAR)
  3198. return value_new_number(y->val._char);
  3199. else if (y->tag == ARRAY && y->val.array->data &&
  3200. is_char_array(y->val.array)) {
  3201. buffer_t *buf = buffer_new();
  3202. for (size_t i = 0; i < y->val.array->length; i++)
  3203. buffer_append(buf, ((value_t *)y->val.array->data[i])->val._char);
  3204. char *s = buffer_read(buf);
  3205. double r = strtod(s, NULL);
  3206. GC_FREE(s);
  3207. return value_new_number(r);
  3208. }
  3209. break;
  3210. case CHAR:
  3211. if (y->tag == NUMBER)
  3212. return value_new_char(y->val.number);
  3213. break;
  3214. }
  3215. }
  3216. return state->udf;
  3217. }
  3218. value_t *verb_print(interpreter_t *state, verb_t *self, value_t *x) {
  3219. char *s = value_str(x);
  3220. fprintf(stdout, "%s", s);
  3221. GC_FREE(s);
  3222. return state->nil;
  3223. }
  3224. value_t *verb_println(interpreter_t *state, verb_t *self, value_t *x) {
  3225. char *s = value_str(x);
  3226. fprintf(stdout, "%s\n", s);
  3227. GC_FREE(s);
  3228. return state->nil;
  3229. }
  3230. value_t *verb_putch(interpreter_t *state, verb_t *self, value_t *x) {
  3231. if (x->tag != CHAR)
  3232. return state->udf;
  3233. fputc(x->val._char, stdout);
  3234. return state->nil;
  3235. }
  3236. value_t *verb_exit(interpreter_t *state, verb_t *self, value_t *x) {
  3237. if (x->tag != NUMBER)
  3238. return state->udf;
  3239. int code = x->val.number;
  3240. exit(code);
  3241. return state->nil;
  3242. }
  3243. value_t *verb_read(interpreter_t *state, verb_t *self, value_t *x) {
  3244. if (x == NUMS[0]) {
  3245. buffer_t *buf = buffer_new();
  3246. size_t size = 0;
  3247. for (;;) {
  3248. int c = fgetc(stdin);
  3249. if (c < 0)
  3250. break;
  3251. buffer_append(buf, c);
  3252. size++;
  3253. }
  3254. char *s = buffer_read(buf);
  3255. list_t *r = list_newk(size);
  3256. for (size_t i = 0; i < size; i++)
  3257. r->data[i] = CHARS[(int)s[i]];
  3258. GC_FREE(s);
  3259. return value_new_array(r);
  3260. } else if (x == NUMS[1])
  3261. return value_new_char((unsigned char)fgetc(stdin));
  3262. else if (x == NUMS[2]) {
  3263. char line[512];
  3264. if (!fgets(line, sizeof(line), stdin))
  3265. return state->udf;
  3266. size_t z = strlen(line);
  3267. list_t *r = list_newk(z);
  3268. for (size_t i = 0; i < z; i++)
  3269. r->data[i] = CHARS[(int)line[i]];
  3270. return value_new_array(r);
  3271. }
  3272. char *path = value_str(x);
  3273. FILE *fd = fopen(path, "rb");
  3274. if (!fd) {
  3275. GC_FREE(path);
  3276. return state->udf;
  3277. }
  3278. fseek(fd, 0, SEEK_END);
  3279. size_t size = ftell(fd);
  3280. fseek(fd, 0, SEEK_SET);
  3281. unsigned char *buf = malloc_checked(size + 1);
  3282. size = fread(buf, sizeof(unsigned char), size, fd);
  3283. fclose(fd);
  3284. GC_FREE(path);
  3285. list_t *r = list_newk(size);
  3286. for (size_t i = 0; i < size; i++)
  3287. r->data[i] = CHARS[buf[i]];
  3288. GC_FREE(buf);
  3289. return value_new_array(r);
  3290. }
  3291. value_t *verb_write(interpreter_t *state, verb_t *self, value_t *x,
  3292. value_t *y) {
  3293. FILE *fd;
  3294. char *path = NULL;
  3295. if (x->tag != ARRAY)
  3296. x = verb_enlist(state, NULL, x);
  3297. if (y == NUMS[0])
  3298. fd = stderr;
  3299. else {
  3300. path = value_str(y);
  3301. fd = fopen(path, "wb");
  3302. if (!fd) {
  3303. GC_FREE(path);
  3304. return NNUMS[0];
  3305. }
  3306. }
  3307. size_t k = 0;
  3308. for (size_t i = 0; i < x->val.array->length; i++) {
  3309. unsigned char c;
  3310. value_t *v = x->val.array->data[i];
  3311. if (v->tag == NUMBER)
  3312. c = fabs(v->val.number);
  3313. else if (v->tag == CHAR)
  3314. c = v->val._char;
  3315. else
  3316. break;
  3317. fputc(c, fd);
  3318. k++;
  3319. }
  3320. fclose(fd);
  3321. if (path)
  3322. GC_FREE(path);
  3323. return value_new_number(k);
  3324. }
  3325. value_t *verb_system(interpreter_t *state, verb_t *self, value_t *x) {
  3326. char *cmd = value_str(x);
  3327. FILE *pd;
  3328. pd = popen(cmd, "r");
  3329. if (!pd) {
  3330. GC_FREE(cmd);
  3331. return state->udf;
  3332. }
  3333. unsigned char *buffer = NULL;
  3334. size_t buffer_size = 0;
  3335. size_t buffer_allocated = 0;
  3336. size_t bytes_received;
  3337. unsigned char chunk[1024];
  3338. for (;;) {
  3339. bytes_received = fread(chunk, 1, 1024, pd);
  3340. if (bytes_received == 0)
  3341. break;
  3342. size_t head = buffer_size;
  3343. buffer_size += bytes_received;
  3344. if (buffer_size > buffer_allocated) {
  3345. buffer_allocated = buffer_size;
  3346. buffer = realloc_checked(buffer, buffer_allocated);
  3347. }
  3348. for (size_t i = 0; i < bytes_received; i++)
  3349. buffer[head + i] = chunk[i];
  3350. if (feof(pd))
  3351. break;
  3352. }
  3353. pclose(pd);
  3354. GC_FREE(cmd);
  3355. list_t *r = list_newk(buffer_size);
  3356. for (size_t i = 0; i < buffer_size; i++)
  3357. r->data[i] = CHARS[buffer[i]];
  3358. GC_FREE(buffer);
  3359. return value_new_array(r);
  3360. }
  3361. struct files_t {
  3362. FILE *in;
  3363. FILE *out;
  3364. };
  3365. typedef struct files_t files_t;
  3366. struct files_chain_t {
  3367. files_t files;
  3368. pid_t pid;
  3369. struct files_chain_t *next;
  3370. };
  3371. typedef struct files_chain_t files_chain_t;
  3372. static files_chain_t *files_chain;
  3373. void _cleanup_pipe(int *pipe) {
  3374. close(pipe[0]);
  3375. close(pipe[1]);
  3376. }
  3377. static int _do_popen2(files_chain_t *link, const char *command) {
  3378. int child_in[2];
  3379. int child_out[2];
  3380. if (0 != pipe(child_in))
  3381. return -1;
  3382. if (0 != pipe(child_out)) {
  3383. _cleanup_pipe(child_in);
  3384. return -1;
  3385. }
  3386. pid_t cpid = link->pid = fork();
  3387. if (0 > cpid) {
  3388. _cleanup_pipe(child_in);
  3389. _cleanup_pipe(child_out);
  3390. return -1;
  3391. }
  3392. if (0 == cpid) {
  3393. if (0 > dup2(child_in[0], 0) || 0 > dup2(child_out[1], 1))
  3394. _Exit(127);
  3395. _cleanup_pipe(child_in);
  3396. _cleanup_pipe(child_out);
  3397. for (files_chain_t *p = files_chain; p; p = p->next) {
  3398. int fd_in = fileno(p->files.in);
  3399. if (fd_in != 0)
  3400. close(fd_in);
  3401. int fd_out = fileno(p->files.out);
  3402. if (fd_out != 1)
  3403. close(fd_out);
  3404. }
  3405. execl("/bin/sh", "sh", "-c", command, (char *)NULL);
  3406. _Exit(127);
  3407. }
  3408. close(child_in[0]);
  3409. close(child_out[1]);
  3410. link->files.in = fdopen(child_in[1], "w");
  3411. link->files.out = fdopen(child_out[0], "r");
  3412. return 0;
  3413. }
  3414. files_t *popen2(const char *command) {
  3415. files_chain_t *link = (files_chain_t *)malloc(sizeof(files_chain_t));
  3416. if (NULL == link)
  3417. return NULL;
  3418. if (0 > _do_popen2(link, command)) {
  3419. free(link);
  3420. return NULL;
  3421. }
  3422. link->next = files_chain;
  3423. files_chain = link;
  3424. return (files_t *)link;
  3425. }
  3426. int pclose2(files_t *fp) {
  3427. files_chain_t **p = &files_chain;
  3428. int found = 0;
  3429. while (*p) {
  3430. if (*p == (files_chain_t *)fp) {
  3431. *p = (*p)->next;
  3432. found = 1;
  3433. break;
  3434. }
  3435. p = &(*p)->next;
  3436. }
  3437. if (!found)
  3438. return -1;
  3439. if (0 > fclose(fp->out)) {
  3440. free((files_chain_t *)fp);
  3441. return -1;
  3442. }
  3443. int status = -1;
  3444. pid_t wait_pid;
  3445. do {
  3446. wait_pid = waitpid(((files_chain_t *)fp)->pid, &status, 0);
  3447. } while (-1 == wait_pid && EINTR == errno);
  3448. free((files_chain_t *)fp);
  3449. if (wait_pid == -1)
  3450. return -1;
  3451. return status;
  3452. }
  3453. value_t *verb_system2(interpreter_t *state, verb_t *self, value_t *x,
  3454. value_t *y) {
  3455. char *cmd = value_str(y);
  3456. files_t *pd;
  3457. pd = popen2(cmd);
  3458. if (pd == NULL) {
  3459. GC_FREE(cmd);
  3460. return state->udf;
  3461. }
  3462. for (size_t i = 0; i < x->val.array->length; i++) {
  3463. unsigned char c;
  3464. value_t *v = x->val.array->data[i];
  3465. if (v->tag == NUMBER)
  3466. c = fabs(v->val.number);
  3467. else if (v->tag == CHAR)
  3468. c = v->val._char;
  3469. else
  3470. break;
  3471. fputc(c, pd->in);
  3472. }
  3473. fflush(pd->in);
  3474. fclose(pd->in);
  3475. unsigned char *buffer = NULL;
  3476. size_t buffer_size = 0;
  3477. size_t buffer_allocated = 0;
  3478. size_t bytes_received;
  3479. unsigned char chunk[1024];
  3480. for (;;) {
  3481. bytes_received = fread(chunk, 1, 1024, pd->out);
  3482. if (bytes_received == 0)
  3483. break;
  3484. size_t head = buffer_size;
  3485. buffer_size += bytes_received;
  3486. if (buffer_size > buffer_allocated) {
  3487. buffer_allocated = buffer_size;
  3488. buffer = realloc_checked(buffer, buffer_allocated);
  3489. }
  3490. for (size_t i = 0; i < bytes_received; i++)
  3491. buffer[head + i] = chunk[i];
  3492. if (feof(pd->out))
  3493. break;
  3494. }
  3495. pclose2(pd);
  3496. GC_FREE(cmd);
  3497. list_t *r = list_newk(buffer_size);
  3498. for (size_t i = 0; i < buffer_size; i++)
  3499. r->data[i] = CHARS[buffer[i]];
  3500. GC_FREE(buffer);
  3501. return value_new_array(r);
  3502. }
  3503. value_t *verb_shl(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  3504. if (x->tag == NUMBER && y->tag == NUMBER)
  3505. return value_new_number(((int)x->val.number) << ((int)y->val.number));
  3506. return _NAN;
  3507. }
  3508. value_t *verb_shr(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  3509. if (x->tag == NUMBER && y->tag == NUMBER)
  3510. return value_new_number(((int)x->val.number) >> ((int)y->val.number));
  3511. return _NAN;
  3512. }
  3513. value_t *verb_xor(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  3514. if (x->tag == NUMBER && y->tag == NUMBER)
  3515. return value_new_number(((int)x->val.number) ^ ((int)y->val.number));
  3516. return _NAN;
  3517. }
  3518. value_t *verb_band(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  3519. if (x->tag == NUMBER && y->tag == NUMBER)
  3520. return value_new_number(((int)x->val.number) & ((int)y->val.number));
  3521. return _NAN;
  3522. }
  3523. list_t *find_primes(uint64_t limit) {
  3524. bool sieve[limit + 1];
  3525. for (uint64_t i = 0; i <= limit; i++)
  3526. sieve[i] = false;
  3527. if (limit > 2)
  3528. sieve[2] = true;
  3529. if (limit > 3)
  3530. sieve[3] = true;
  3531. for (uint64_t x = 1; x * x <= limit; x++)
  3532. for (uint64_t y = 1; y * y <= limit; y++) {
  3533. uint64_t n = (4 * x * x) + (y * y);
  3534. if (n <= limit && (n % 12 == 1 || n % 12 == 5))
  3535. sieve[n] ^= true;
  3536. n = (3 * x * x) + (y * y);
  3537. if (n <= limit && n % 12 == 7)
  3538. sieve[n] ^= true;
  3539. n = (3 * x * x) - (y * y);
  3540. if (x > y && n <= limit && n % 12 == 11)
  3541. sieve[n] ^= true;
  3542. }
  3543. for (uint64_t r = 5; r * r <= limit; r++)
  3544. if (sieve[r])
  3545. for (int i = r * r; i <= limit; i += r * r)
  3546. sieve[i] = false;
  3547. list_t *r = list_new();
  3548. for (uint64_t a = 1; a <= limit; a++)
  3549. if (sieve[a])
  3550. list_push(r, value_new_number(a));
  3551. return r;
  3552. }
  3553. value_t *verb_primes(interpreter_t *state, verb_t *self, value_t *x) {
  3554. if (x->tag == NUMBER && !is_bad_num(x->val.number))
  3555. return value_new_array(find_primes(fabs(x->val.number) + 1));
  3556. return state->udf;
  3557. }
  3558. value_t *verb_parts(interpreter_t *state, verb_t *self, value_t *x,
  3559. value_t *y) {
  3560. if (x->tag != NUMBER)
  3561. return state->udf;
  3562. if (y->tag != ARRAY)
  3563. y = verb_enlist(state, NULL, y);
  3564. else if (!y->val.array->data)
  3565. return y;
  3566. if (is_bad_num(x->val.number) || x->val.number < 1)
  3567. return y;
  3568. size_t np = fabs(x->val.number);
  3569. list_t *r = list_newk(np);
  3570. size_t rp = 0;
  3571. for (ssize_t i = np; i > 0; i--) {
  3572. size_t k = ceil(((double)y->val.array->length) / (double)i);
  3573. r->data[rp++] = verb_take(state, NULL, value_new_number(k), y);
  3574. y = verb_drop(state, NULL, value_new_number(k), y);
  3575. }
  3576. return value_new_array(r);
  3577. }
  3578. value_t *verb_bor(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  3579. if (x->tag == NUMBER && y->tag == NUMBER)
  3580. return value_new_number(((int)x->val.number) | ((int)y->val.number));
  3581. return _NAN;
  3582. }
  3583. value_t *verb_bnot(interpreter_t *state, verb_t *self, value_t *x) {
  3584. if (x->tag == NUMBER)
  3585. return value_new_number(~(int)x->val.number);
  3586. return _NAN;
  3587. }
  3588. list_t *prime_factors(double n) {
  3589. list_t *factors = list_new();
  3590. double divisor = 2;
  3591. while (n >= 2) {
  3592. if (fmod(n, divisor) == 0) {
  3593. list_push(factors, value_new_number(divisor));
  3594. n /= divisor;
  3595. } else
  3596. divisor++;
  3597. }
  3598. return factors;
  3599. }
  3600. value_t *verb_factors(interpreter_t *state, verb_t *self, value_t *x) {
  3601. if (x->tag == NUMBER && !is_bad_num(x->val.number))
  3602. return value_new_array(prime_factors(x->val.number));
  3603. return state->udf;
  3604. }
  3605. value_t *verb_combine(interpreter_t *state, verb_t *self, value_t *x,
  3606. value_t *y) {
  3607. if (x->tag == NUMBER && y->tag == NUMBER && !is_bad_num(x->val.number) &&
  3608. !is_bad_num(y->val.number)) {
  3609. x = verb_base(state, NULL, NUMS[10], x);
  3610. y = verb_base(state, NULL, NUMS[10], y);
  3611. value_t *n = verb_join(state, NULL, x, y);
  3612. return verb_unbase(state, NULL, NUMS[10], n);
  3613. }
  3614. return _NAN;
  3615. }
  3616. value_t *verb_outof(interpreter_t *state, verb_t *self, value_t *x,
  3617. value_t *y) {
  3618. if (x->tag == NUMBER && y->tag == NUMBER && !is_bad_num(x->val.number) &&
  3619. !is_bad_num(y->val.number)) {
  3620. uint64_t a = (uint64_t)fabs(x->val.number);
  3621. uint64_t b = (uint64_t)fabs(y->val.number);
  3622. if (a == 0)
  3623. return NUMS[1];
  3624. if (b == 0)
  3625. return NUMS[0];
  3626. return value_new_number((double)factorial(b) /
  3627. (factorial(a) * (a >= b ? 1 : factorial(b - a))));
  3628. }
  3629. return _NAN;
  3630. }
  3631. value_t *verb_sort(interpreter_t *state, verb_t *self, value_t *x) {
  3632. value_t *i = verb_gradeup(state, NULL, x);
  3633. return together(state, state->at, x, i, 0, 0, state->at->rank[1],
  3634. state->at->rank[2]);
  3635. }
  3636. value_t *verb_unsort(interpreter_t *state, verb_t *self, value_t *x) {
  3637. value_t *i = verb_gradedown(state, NULL, x);
  3638. return together(state, state->at, x, i, 0, 0, state->at->rank[1],
  3639. state->at->rank[2]);
  3640. }
  3641. value_t *interpreter_run(interpreter_t *state, char *program);
  3642. value_t *verb_eval(interpreter_t *state, verb_t *self, value_t *x) {
  3643. char *s = value_str(x);
  3644. guard_t *g = guard();
  3645. if (setjmp(g->lb)) {
  3646. unguard();
  3647. GC_FREE(s);
  3648. return state->udf;
  3649. }
  3650. value_t *v = interpreter_run(state, s);
  3651. GC_FREE(s);
  3652. unguard();
  3653. return v;
  3654. }
  3655. void jkexec(interpreter_t *state, FILE *fd, bool isrepl, char **s);
  3656. value_t *verb_import(interpreter_t *state, verb_t *self, value_t *x) {
  3657. char *path = value_str(x);
  3658. FILE *fd = fopen(path, "rb");
  3659. if (!fd) {
  3660. GC_FREE(path);
  3661. return state->udf;
  3662. }
  3663. char *s = NULL;
  3664. jkexec(state, fd, false, &s);
  3665. if (s)
  3666. GC_FREE(s);
  3667. fclose(fd);
  3668. GC_FREE(path);
  3669. return state->nil;
  3670. }
  3671. value_t *verb_foreign(interpreter_t *state, verb_t *self, value_t *x,
  3672. value_t *y) {
  3673. if (x->tag != ARRAY)
  3674. x = verb_enlist(state, NULL, x);
  3675. char *pth = value_str(y);
  3676. char *lib;
  3677. char *sig;
  3678. char *fun;
  3679. lib = strtok(pth, ":");
  3680. if (!lib)
  3681. return state->udf;
  3682. sig = strtok(NULL, ":");
  3683. if (!sig)
  3684. return state->udf;
  3685. fun = strtok(NULL, ":");
  3686. if (!fun)
  3687. return state->udf;
  3688. size_t argc = strlen(sig);
  3689. if (argc < 1)
  3690. return state->udf;
  3691. argc--;
  3692. if (argc != x->val.array->length)
  3693. return state->udf;
  3694. ffi_cif cif;
  3695. ffi_type *ret;
  3696. ffi_type *args[argc];
  3697. void *values[argc];
  3698. void *pool[argc];
  3699. size_t fc = 0;
  3700. void *retv = NULL;
  3701. char rett;
  3702. size_t retvsz = 0;
  3703. for (int i = 0; i < strlen(sig); i++) {
  3704. ffi_type *t;
  3705. void *v;
  3706. switch (sig[i]) {
  3707. case '$':
  3708. t = &ffi_type_pointer;
  3709. break;
  3710. case 'p':
  3711. t = &ffi_type_pointer;
  3712. break;
  3713. case 'v':
  3714. if (i != 0)
  3715. goto cleanup;
  3716. t = &ffi_type_void;
  3717. break;
  3718. case 'i':
  3719. t = &ffi_type_sint;
  3720. break;
  3721. case 'l':
  3722. t = &ffi_type_slong;
  3723. break;
  3724. case 'f':
  3725. t = &ffi_type_float;
  3726. break;
  3727. case 'd':
  3728. t = &ffi_type_double;
  3729. break;
  3730. case 'c':
  3731. t = &ffi_type_uchar;
  3732. break;
  3733. default:
  3734. goto cleanup;
  3735. }
  3736. if (i == 0) {
  3737. rett = sig[0];
  3738. ret = t;
  3739. switch (rett) {
  3740. case '$':
  3741. case '@':
  3742. retvsz = sizeof(char *);
  3743. break;
  3744. case 'p':
  3745. retvsz = sizeof(void *);
  3746. break;
  3747. case 'v':
  3748. retvsz = 0;
  3749. break;
  3750. case 'i':
  3751. retvsz = sizeof(int);
  3752. break;
  3753. case 'l':
  3754. retvsz = sizeof(long);
  3755. break;
  3756. case 'f':
  3757. retvsz = sizeof(float);
  3758. break;
  3759. case 'd':
  3760. retvsz = sizeof(double);
  3761. break;
  3762. case 'c':
  3763. retvsz = sizeof(unsigned char);
  3764. break;
  3765. }
  3766. } else {
  3767. switch (sig[i]) {
  3768. case '$':
  3769. case '@': {
  3770. value_t *vt = x->val.array->data[i - 1];
  3771. pool[i - 1] = value_str(vt);
  3772. v = pool[i - 1];
  3773. fc++;
  3774. } break;
  3775. case 'p': {
  3776. void *_pv;
  3777. value_t *vt = x->val.array->data[i - 1];
  3778. if (vt->tag != NUMBER)
  3779. goto cleanup;
  3780. _pv = (void *)(size_t)fabs(vt->val.number);
  3781. pool[i - 1] = malloc_checked(sizeof(void *));
  3782. memcpy(pool[i - 1], &_pv, sizeof(void *));
  3783. v = pool[i - 1];
  3784. fc++;
  3785. } break;
  3786. case 'i': {
  3787. int _iv;
  3788. value_t *vt = x->val.array->data[i - 1];
  3789. if (vt->tag != NUMBER)
  3790. goto cleanup;
  3791. _iv = (int)vt->val.number;
  3792. pool[i - 1] = malloc_checked(sizeof(int));
  3793. memcpy(pool[i - 1], &_iv, sizeof(int));
  3794. v = pool[i - 1];
  3795. fc++;
  3796. } break;
  3797. case 'l': {
  3798. long _lv;
  3799. value_t *_vt = x->val.array->data[i - 1];
  3800. if (_vt->tag != NUMBER)
  3801. goto cleanup;
  3802. _lv = (long)_vt->val.number;
  3803. pool[i - 1] = malloc_checked(sizeof(long));
  3804. memcpy(pool[i - 1], &_lv, sizeof(long));
  3805. v = pool[i - 1];
  3806. fc++;
  3807. } break;
  3808. case 'f': {
  3809. float _fv;
  3810. value_t *_vt = x->val.array->data[i - 1];
  3811. if (_vt->tag != NUMBER)
  3812. goto cleanup;
  3813. _fv = (float)_vt->val.number;
  3814. pool[i - 1] = malloc_checked(sizeof(float));
  3815. memcpy(pool[i - 1], &_fv, sizeof(float));
  3816. v = pool[i - 1];
  3817. fc++;
  3818. } break;
  3819. case 'd': {
  3820. double _dv;
  3821. value_t *_vt = x->val.array->data[i - 1];
  3822. if (_vt->tag != NUMBER)
  3823. goto cleanup;
  3824. _dv = (double)_vt->val.number;
  3825. pool[i - 1] = malloc_checked(sizeof(double));
  3826. memcpy(pool[i - 1], &_dv, sizeof(double));
  3827. v = pool[i - 1];
  3828. fc++;
  3829. } break;
  3830. case 'c': {
  3831. unsigned char _cv;
  3832. value_t *_vt = x->val.array->data[i - 1];
  3833. if (_vt->tag != CHAR)
  3834. goto cleanup;
  3835. _cv = (unsigned char)_vt->val._char;
  3836. pool[i - 1] = malloc_checked(sizeof(unsigned char));
  3837. memcpy(pool[i - 1], &_cv, sizeof(unsigned char));
  3838. v = pool[i - 1];
  3839. fc++;
  3840. } break;
  3841. }
  3842. args[i - 1] = t;
  3843. values[i - 1] = v;
  3844. }
  3845. }
  3846. void *dlh = dlopen(lib, RTLD_LAZY);
  3847. if (!dlh)
  3848. goto cleanup;
  3849. void *exfn = dlsym(dlh, fun);
  3850. char *e = dlerror();
  3851. if (!exfn || e)
  3852. goto cleanup;
  3853. if (ffi_prep_cif(&cif, FFI_DEFAULT_ABI, argc, ret, args) != FFI_OK)
  3854. goto cleanup;
  3855. if (retvsz)
  3856. retv = malloc_checked(retvsz);
  3857. ffi_call(&cif, exfn, retv, values);
  3858. dlclose(dlh);
  3859. value_t *rv = state->nil;
  3860. switch (rett) {
  3861. case 'v':
  3862. break;
  3863. case '$': {
  3864. char *s = *(char **)retv;
  3865. size_t z = strlen(s);
  3866. list_t *l = list_newk(z);
  3867. for (size_t i = 0; i < z; i++)
  3868. l->data[i] = CHARS[(int)s[i]];
  3869. rv = value_new_array(l);
  3870. } break;
  3871. case '@': {
  3872. char *s = *(char **)retv;
  3873. size_t z = strlen(s);
  3874. list_t *l = list_newk(z);
  3875. for (size_t i = 0; i < z; i++)
  3876. l->data[i] = CHARS[(int)s[i]];
  3877. rv = value_new_array(l);
  3878. free(s);
  3879. } break;
  3880. case 'p':
  3881. rv = value_new_number((size_t)*(void **)retv);
  3882. break;
  3883. case 'i':
  3884. rv = value_new_number(*(int *)retv);
  3885. break;
  3886. case 'l':
  3887. rv = value_new_number(*(long *)retv);
  3888. break;
  3889. case 'f':
  3890. rv = value_new_number(*(float *)retv);
  3891. break;
  3892. case 'd':
  3893. rv = value_new_number(*(double *)retv);
  3894. break;
  3895. case 'c':
  3896. rv = value_new_char(*(unsigned char *)retv);
  3897. break;
  3898. }
  3899. GC_FREE(retv);
  3900. for (size_t i = 0; i < fc; i++)
  3901. GC_FREE(pool[i]);
  3902. return rv;
  3903. cleanup:
  3904. for (size_t i = 0; i < fc; i++)
  3905. GC_FREE(pool[i]);
  3906. return state->udf;
  3907. }
  3908. value_t *verb_explode(interpreter_t *state, verb_t *self, value_t *x,
  3909. value_t *y) {
  3910. char *del = value_show(x);
  3911. char *s = value_str(y);
  3912. size_t dell = strlen(del);
  3913. size_t sl = strlen(s);
  3914. list_t *r = list_new();
  3915. list_t *t = list_new();
  3916. for (size_t i = 0; i < sl; i++) {
  3917. if (strncmp(&s[i], del, dell) == 0) {
  3918. list_push(r, value_new_array(t));
  3919. t = list_new();
  3920. i += dell - 1;
  3921. continue;
  3922. }
  3923. list_push(t, CHARS[(int)s[i]]);
  3924. }
  3925. GC_FREE(s);
  3926. GC_FREE(del);
  3927. list_push(r, value_new_array(t));
  3928. return value_new_array(r);
  3929. }
  3930. value_t *verb_implode(interpreter_t *state, verb_t *self, value_t *x,
  3931. value_t *y) {
  3932. if (y->tag != ARRAY || !y->val.array->data)
  3933. return y;
  3934. char *del = value_show(x);
  3935. list_t *r = list_new();
  3936. for (size_t i = 0; i < y->val.array->length; i++) {
  3937. char *s = value_show(y->val.array->data[i]);
  3938. char *_s = s;
  3939. while (*_s)
  3940. list_push(r, CHARS[(int)(*_s++)]);
  3941. GC_FREE(s);
  3942. if (i != y->val.array->length - 1) {
  3943. char *s = del;
  3944. while (*s)
  3945. list_push(r, CHARS[(int)(*s++)]);
  3946. }
  3947. }
  3948. GC_FREE(del);
  3949. return value_new_array(r);
  3950. }
  3951. value_t *verb_tackleft(interpreter_t *state, verb_t *self, value_t *x,
  3952. value_t *y) {
  3953. if (y->tag != ARRAY)
  3954. y = verb_enlist(state, NULL, y);
  3955. list_t *r = list_newk(y->val.array->length + 1);
  3956. r->data[0] = x;
  3957. for (size_t i = 0; i < y->val.array->length; i++)
  3958. r->data[i + 1] = y->val.array->data[i];
  3959. return value_new_array(r);
  3960. }
  3961. value_t *verb_setrecdepth(interpreter_t *state, verb_t *self, value_t *x) {
  3962. if (x->tag != NUMBER)
  3963. return state->udf;
  3964. size_t ov = max_rec_depth;
  3965. size_t v = (size_t)fabs(x->val.number);
  3966. if (v < 1)
  3967. v = 1;
  3968. max_rec_depth = v;
  3969. return value_new_number(ov);
  3970. }
  3971. value_t *verb_tackright(interpreter_t *state, verb_t *self, value_t *x,
  3972. value_t *y) {
  3973. if (y->tag != ARRAY)
  3974. y = verb_enlist(state, NULL, y);
  3975. list_t *r = list_newk(y->val.array->length + 1);
  3976. for (size_t i = 0; i < y->val.array->length; i++)
  3977. r->data[i] = y->val.array->data[i];
  3978. r->data[y->val.array->length] = x;
  3979. return value_new_array(r);
  3980. }
  3981. value_t *verb_eye(interpreter_t *state, verb_t *self, value_t *x) {
  3982. if (x->tag == NUMBER && !is_bad_num(x->val.number)) {
  3983. size_t k = fabs(x->val.number);
  3984. list_t *r = list_newk(k);
  3985. for (size_t i = 0; i < k; i++) {
  3986. list_t *rw = list_newk(k);
  3987. for (size_t j = 0; j < k; j++)
  3988. rw->data[j] = NUMS[i == j];
  3989. r->data[i] = value_new_array(rw);
  3990. }
  3991. return value_new_array(r);
  3992. }
  3993. return state->udf;
  3994. }
  3995. value_t *verb_infix(interpreter_t *state, verb_t *self, value_t *x) {
  3996. return verb_behead(state, NULL, verb_prefixes(state, NULL, x));
  3997. }
  3998. value_t *verb_value(interpreter_t *state, verb_t *self, value_t *x) {
  3999. char *s = value_str(x);
  4000. value_t *r = table_get(state->env, s);
  4001. GC_FREE(s);
  4002. return r ? r : state->udf;
  4003. }
  4004. value_t *verb_lines(interpreter_t *state, verb_t *self, value_t *x) {
  4005. char *s = value_str(x);
  4006. size_t sl = strlen(s);
  4007. list_t *r = list_new();
  4008. list_t *t = list_new();
  4009. for (size_t i = 0; i < sl; i++) {
  4010. if (s[i] == '\n') {
  4011. list_push(r, value_new_array(t));
  4012. t = list_new();
  4013. continue;
  4014. }
  4015. list_push(t, CHARS[(int)s[i]]);
  4016. }
  4017. GC_FREE(s);
  4018. list_push(r, value_new_array(t));
  4019. return value_new_array(r);
  4020. }
  4021. value_t *verb_delete(interpreter_t *state, verb_t *self, value_t *x,
  4022. value_t *y) {
  4023. if (x->tag != NUMBER)
  4024. return state->udf;
  4025. if (y->tag != ARRAY)
  4026. y = verb_enlist(state, NULL, y);
  4027. else if (!y->val.array->data)
  4028. return y;
  4029. size_t z = y->val.array->length;
  4030. ssize_t index = trunc(x->val.number);
  4031. if (index < 0)
  4032. index += ((ssize_t)z);
  4033. if (index < 0 || index >= z)
  4034. return y;
  4035. list_t *r = list_newk(z - 1);
  4036. size_t ri = 0;
  4037. for (size_t i = 0; i < z; i++)
  4038. if (i == index)
  4039. continue;
  4040. else
  4041. r->data[ri++] = y->val.array->data[i];
  4042. return value_new_array(r);
  4043. }
  4044. value_t *verb_deleteInplace(interpreter_t *state, verb_t *self, value_t *x,
  4045. value_t *y) {
  4046. if (x->tag != NUMBER)
  4047. return state->udf;
  4048. if (y->tag != ARRAY)
  4049. return state->udf;
  4050. else if (!y->val.array->data)
  4051. return y;
  4052. size_t z = y->val.array->length;
  4053. ssize_t index = trunc(x->val.number);
  4054. if (index < 0)
  4055. index += ((ssize_t)z);
  4056. if (index < 0 || index >= z)
  4057. return y;
  4058. list_t *r = list_newk(z - 1);
  4059. size_t ri = 0;
  4060. for (size_t i = 0; i < z; i++)
  4061. if (i == index)
  4062. continue;
  4063. else
  4064. r->data[ri++] = y->val.array->data[i];
  4065. y->val.array = r;
  4066. return y;
  4067. }
  4068. value_t *verb_rematch(interpreter_t *state, verb_t *self, value_t *x,
  4069. value_t *y) {
  4070. char *pat = value_str(x);
  4071. char *s = value_str(y);
  4072. pcre *re;
  4073. const char *e;
  4074. int eo;
  4075. if (!(re = pcre_compile(pat, 0, &e, &eo, NULL)))
  4076. goto fail;
  4077. int rc = pcre_exec(re, NULL, s, strlen(s), 0, 0, NULL, 0);
  4078. pcre_free(re);
  4079. GC_FREE(pat);
  4080. GC_FREE(s);
  4081. return NUMS[rc >= 0];
  4082. fail:
  4083. GC_FREE(pat);
  4084. GC_FREE(s);
  4085. return state->udf;
  4086. }
  4087. value_t *verb_show(interpreter_t *state, verb_t *self, value_t *x) {
  4088. char *s = value_str(x);
  4089. list_t *r = list_new();
  4090. for (size_t i = 0; i < strlen(s); i++)
  4091. list_push(r, value_new_char(s[i]));
  4092. GC_FREE(s);
  4093. return value_new_array(r);
  4094. }
  4095. value_t *verb_extract(interpreter_t *state, verb_t *self, value_t *x,
  4096. value_t *y) {
  4097. char *pat = value_str(x);
  4098. char *s = value_str(y);
  4099. size_t len = strlen(s);
  4100. pcre *re;
  4101. const char *e;
  4102. int eo;
  4103. if (!(re = pcre_compile(pat, 0, &e, &eo, NULL)))
  4104. goto fail;
  4105. int ov[128 * 3];
  4106. int rc;
  4107. list_t *r = list_new();
  4108. unsigned int of = 0;
  4109. while (of < len &&
  4110. (rc = pcre_exec(re, 0, s, len, of, 0, ov, sizeof(ov))) >= 0) {
  4111. if (rc == 0)
  4112. rc = sizeof(ov) / 3;
  4113. for (int i = 1; i < rc; i++) {
  4114. char *ss = s + ov[2 * i];
  4115. int sl = ov[2 * i + 1] - ov[2 * i];
  4116. if (sl == 0) {
  4117. list_push(r, _UNIT);
  4118. continue;
  4119. }
  4120. list_t *l = list_newk(sl);
  4121. for (int j = 0; j < sl; j++)
  4122. l->data[j] = CHARS[(int)ss[j]];
  4123. list_push(r, value_new_array(l));
  4124. }
  4125. of = ov[1];
  4126. }
  4127. pcre_free(re);
  4128. GC_FREE(pat);
  4129. GC_FREE(s);
  4130. return value_new_array(r);
  4131. fail:
  4132. GC_FREE(pat);
  4133. GC_FREE(s);
  4134. return state->udf;
  4135. }
  4136. value_t *verb_udf1(interpreter_t *state, verb_t *self, value_t *x) {
  4137. return state->udf;
  4138. }
  4139. value_t *verb_udf2(interpreter_t *state, verb_t *self, value_t *x, value_t *y) {
  4140. return state->udf;
  4141. }
  4142. #define X UINT_MAX
  4143. #define DEFVERB(__symb, __rm, __rl, __rr, __monad, __dyad) \
  4144. {__symb, {__rm, __rl, __rr}, NULL, false, \
  4145. false, verb_##__monad, verb_##__dyad}
  4146. #define DEFVERBD(__symb, __rm, __rl, __rr, __monad, __dyad) \
  4147. {__symb ".", {__rm, __rl, __rr}, NULL, false, \
  4148. false, verb_##__monad, verb_##__dyad}
  4149. #define DEFVERBC(__symb, __rm, __rl, __rr, __monad, __dyad) \
  4150. {__symb ":", {__rm, __rl, __rr}, NULL, false, \
  4151. false, verb_##__monad, verb_##__dyad}
  4152. verb_t VERBS[] = {DEFVERB(":", 0, 0, 0, const, bind),
  4153. DEFVERBC(":", X, 0, 0, unbind, obverse),
  4154. DEFVERB("+", 0, X, X, flip, plus),
  4155. DEFVERBD("+", X, X, X, fibonacci, gcd),
  4156. DEFVERBC("+", X, X, X, sin, combine),
  4157. DEFVERB("-", X, X, X, negate, minus),
  4158. DEFVERBD("-", X, X, X, atan, atan2),
  4159. DEFVERB("*", 0, X, X, first, times),
  4160. DEFVERBD("*", X, X, X, factorial, lcm),
  4161. DEFVERBC("*", X, X, 0, double, replicate),
  4162. DEFVERB("%", X, X, X, reciprocal, divide),
  4163. DEFVERBD("%", X, X, X, sqrt, root),
  4164. DEFVERBC("%", X, X, X, halve, idivide),
  4165. DEFVERB("!", X, X, X, enum, mod),
  4166. DEFVERBD("!", X, X, X, iota, range),
  4167. DEFVERBC("!", 0, X, 0, odometer, chunks),
  4168. DEFVERB("^", X, X, X, exp, power),
  4169. DEFVERBD("^", X, X, X, nlog, log),
  4170. DEFVERB("=", 0, X, X, permute, equals),
  4171. DEFVERBD("=", 0, 0, 0, occurences, mask),
  4172. DEFVERBC("=", 0, 0, 0, classify, equals),
  4173. DEFVERB("~", X, X, X, not, not_equals),
  4174. DEFVERBD("~", X, 0, 0, sign, insert),
  4175. DEFVERBC("~", 0, 0, 0, not, not_equals),
  4176. DEFVERB("<", X, X, X, pred, less),
  4177. DEFVERBD("<", X, X, X, floor, lesseq),
  4178. DEFVERBC("<", 0, X, 0, gradedown, nudge_left),
  4179. DEFVERB(">", X, X, X, succ, greater),
  4180. DEFVERBD(">", X, X, X, ceil, greatereq),
  4181. DEFVERBC(">", 0, X, 0, gradeup, nudge_right),
  4182. DEFVERB(",", 0, 0, 0, enlist, join),
  4183. DEFVERBD(",", X, 0, 0, enlist, enpair),
  4184. DEFVERB("#", 0, X, 0, count, take),
  4185. DEFVERBD("#", 0, 0, 0, where, copy),
  4186. DEFVERBC("#", 0, 0, 0, group, buckets),
  4187. DEFVERB("_", 0, X, 0, nub, drop),
  4188. DEFVERBD("_", 0, X, 0, unbits, unbase),
  4189. DEFVERBC("_", X, X, X, bits, base),
  4190. DEFVERB("?", 0, 0, 0, unique, find),
  4191. DEFVERB("&", 0, X, X, flatten, minand),
  4192. DEFVERB("|", 0, X, X, reverse, maxor),
  4193. DEFVERBD("|", X, X, 0, round, rotate),
  4194. DEFVERBC("|", 0, X, 0, depth, windows),
  4195. DEFVERB("@", X, 0, X, abs, at),
  4196. DEFVERBD("@", 0, 0, 0, shuffle, member),
  4197. DEFVERBC("@", 0, 0, 0, infix, indexof),
  4198. DEFVERB("{", 0, 0, 0, head, bin),
  4199. DEFVERBD("{", 0, 0, 0, tail, cut),
  4200. DEFVERBC("{", 0, X, X, prefixes, shl),
  4201. DEFVERB("}", 0, X, X, behead, xor),
  4202. DEFVERBD("}", 0, 0, 0, curtail, band),
  4203. DEFVERBC("}", 0, X, X, suffixes, shr),
  4204. DEFVERB("[", X, 0, 0, factors, left),
  4205. DEFVERBD("[", X, X, X, bnot, bor),
  4206. DEFVERBC("[", X, X, 0, primes, parts),
  4207. DEFVERB("]", 0, 0, 0, same, right),
  4208. DEFVERBD("]", 0, X, X, sort, outof),
  4209. DEFVERBC("]", 0, 0, 0, unsort, explode),
  4210. DEFVERBD("`", 0, 0, 0, symbol, apply1),
  4211. DEFVERBC("`", 0, 0, 0, square, apply2),
  4212. DEFVERB("$", 0, 0, 0, shape, reshape),
  4213. DEFVERBD("$", 0, 0, 0, repr, format),
  4214. DEFVERBC("$", X, 0, 0, eye, implode),
  4215. DEFVERBD("d", 0, X, 0, udf1, delete),
  4216. DEFVERBD("D", 0, X, 0, udf1, deleteInplace),
  4217. DEFVERBD("p", 0, 0, 0, print, udf2),
  4218. DEFVERBD("P", 0, 0, 0, println, udf2),
  4219. DEFVERBD("c", X, 0, 0, putch, udf2),
  4220. DEFVERBD("s", 0, 0, 0, selfref1, selfref2),
  4221. DEFVERBD("F", 0, 0, 0, read, write),
  4222. DEFVERBD("r", 0, X, X, deal, roll),
  4223. DEFVERBD("t", 0, 0, 0, type, cast),
  4224. DEFVERBD("E", 0, 0, 0, exit, udf2),
  4225. DEFVERBD("y", 0, 0, 0, system, system2),
  4226. DEFVERBD("e", 0, 0, 0, eval, udf2),
  4227. DEFVERBD("i", 0, 0, 0, import, foreign),
  4228. DEFVERBD("L", 0, 0, 0, lines, tackleft),
  4229. DEFVERBD("R", X, 0, 0, setrecdepth, tackright),
  4230. DEFVERBD("v", 0, 0, 0, value, udf2),
  4231. DEFVERBD("x", 0, 0, 0, show, rematch),
  4232. DEFVERBD("X", 0, 0, 0, udf1, extract)};
  4233. value_t *_adverb_fold_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4234. if (x->tag != ARRAY || !x->val.array->data)
  4235. return x;
  4236. value_t *_v = self->bonds->data[0];
  4237. if (_v->tag != VERB)
  4238. return state->udf;
  4239. verb_t *v = _v->val.verb;
  4240. value_t *t = x->val.array->data[0];
  4241. list_t *tx = x->val.array;
  4242. for (size_t i = 1; i < tx->length; i++)
  4243. t = together(state, v, t, tx->data[i], 0, 0, v->rank[1], v->rank[2]);
  4244. return t;
  4245. }
  4246. value_t *_adverb_fold_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4247. value_t *y) {
  4248. if (y->tag != ARRAY)
  4249. y = verb_enlist(state, NULL, y);
  4250. else if (!y->val.array->data)
  4251. return x;
  4252. value_t *_v = self->bonds->data[0];
  4253. if (_v->tag != VERB)
  4254. return state->udf;
  4255. verb_t *v = _v->val.verb;
  4256. value_t *t = x;
  4257. list_t *ty = y->val.array;
  4258. for (size_t i = 0; i < ty->length; i++)
  4259. t = together(state, v, t, ty->data[i], 0, 0, v->rank[1], v->rank[2]);
  4260. return t;
  4261. }
  4262. value_t *_adverb_scan_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4263. if (x->tag != ARRAY || !x->val.array->data)
  4264. return x;
  4265. value_t *_v = self->bonds->data[0];
  4266. if (_v->tag != VERB)
  4267. return state->udf;
  4268. verb_t *v = _v->val.verb;
  4269. list_t *r = list_new();
  4270. value_t *t = x->val.array->data[0];
  4271. list_t *tx = x->val.array;
  4272. list_push(r, t);
  4273. for (size_t i = 1; i < tx->length; i++) {
  4274. t = together(state, v, t, tx->data[i], 0, 0, v->rank[1], v->rank[2]);
  4275. list_push(r, t);
  4276. }
  4277. return value_new_array(r);
  4278. }
  4279. value_t *_adverb_scan_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4280. value_t *y) {
  4281. if (y->tag != ARRAY || !y->val.array->data)
  4282. return y;
  4283. value_t *_v = self->bonds->data[0];
  4284. if (_v->tag != VERB)
  4285. return state->udf;
  4286. verb_t *v = _v->val.verb;
  4287. list_t *r = list_new();
  4288. value_t *t = x;
  4289. list_t *ty = y->val.array;
  4290. list_push(r, t);
  4291. for (size_t i = 1; i < ty->length; i++) {
  4292. t = together(state, v, t, ty->data[i], 0, 0, v->rank[1], v->rank[2]);
  4293. list_push(r, t);
  4294. }
  4295. return value_new_array(r);
  4296. }
  4297. value_t *_adverb_each_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4298. value_t *_v = self->bonds->data[0];
  4299. if (_v->tag != VERB)
  4300. return state->udf;
  4301. verb_t *v = _v->val.verb;
  4302. if (x->tag != ARRAY)
  4303. return each_rank(state, v, x, 0, 1);
  4304. if (!x->val.array->data)
  4305. return x;
  4306. return each_rank(state, v, x, 0, 1);
  4307. }
  4308. value_t *_adverb_each_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4309. value_t *y) {
  4310. value_t *_v = self->bonds->data[0];
  4311. if (_v->tag != VERB)
  4312. return state->udf;
  4313. verb_t *v = _v->val.verb;
  4314. if (x->tag != ARRAY)
  4315. x = verb_enlist(state, NULL, x);
  4316. if (y->tag != ARRAY)
  4317. y = verb_enlist(state, NULL, y);
  4318. list_t *r = list_new();
  4319. list_t *tx = x->val.array;
  4320. list_t *ty = y->val.array;
  4321. for (size_t i = 0; i < tx->length && i < ty->length; i++)
  4322. list_push(r, together(state, v, tx->data[i], ty->data[i], 0, 0, v->rank[1],
  4323. v->rank[2]));
  4324. return value_new_array(r);
  4325. }
  4326. value_t *_adverb_converge_monad(interpreter_t *state, verb_t *self,
  4327. value_t *x) {
  4328. value_t *_v = self->bonds->data[0];
  4329. if (_v->tag != VERB)
  4330. return state->udf;
  4331. verb_t *v = _v->val.verb;
  4332. value_t *t;
  4333. for (;;) {
  4334. t = x;
  4335. x = each_rank(state, v, x, 0, v->rank[0]);
  4336. if (value_equals(x, t))
  4337. break;
  4338. }
  4339. return x;
  4340. }
  4341. verb_t *conjunction_bond(interpreter_t *state, value_t *x, value_t *y);
  4342. value_t *_adverb_converge_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4343. value_t *y) {
  4344. value_t *_v = self->bonds->data[0];
  4345. if (_v->tag != VERB)
  4346. return state->udf;
  4347. verb_t *v = _v->val.verb;
  4348. if (y->tag != ARRAY)
  4349. return together(state, v, y, x, 0, 0, v->rank[1], v->rank[2]);
  4350. if (!y->val.array->data)
  4351. return x;
  4352. v = conjunction_bond(state, value_new_verb(v), x);
  4353. return each_rank(state, v, y, 0, 1);
  4354. }
  4355. value_t *_adverb_converges_monad(interpreter_t *state, verb_t *self,
  4356. value_t *x) {
  4357. value_t *_v = self->bonds->data[0];
  4358. if (_v->tag != VERB)
  4359. return state->udf;
  4360. list_t *r = list_new();
  4361. value_t *t;
  4362. list_push(r, x);
  4363. for (;;) {
  4364. t = x;
  4365. x = apply_monad(state, _v, x);
  4366. if (value_equals(x, t))
  4367. break;
  4368. list_push(r, x);
  4369. }
  4370. return value_new_array(r);
  4371. }
  4372. value_t *_adverb_converges_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4373. value_t *y) {
  4374. value_t *_v = self->bonds->data[0];
  4375. if (_v->tag != VERB)
  4376. return state->udf;
  4377. verb_t *v = _v->val.verb;
  4378. if (y->tag != ARRAY)
  4379. return together(state, v, y, x, 0, 0, v->rank[1], v->rank[2]);
  4380. if (!y->val.array->data)
  4381. return x;
  4382. v = conjunction_bond(state, x, value_new_verb(v));
  4383. return each_rank(state, v, y, 0, 1);
  4384. }
  4385. value_t *_adverb_eachprior_monad(interpreter_t *state, verb_t *self,
  4386. value_t *x) {
  4387. if (x->tag != ARRAY || x->val.array->length < 2)
  4388. return x;
  4389. value_t *_v = self->bonds->data[0];
  4390. if (_v->tag != VERB)
  4391. return state->udf;
  4392. verb_t *v = _v->val.verb;
  4393. list_t *r = list_new();
  4394. for (size_t i = 1; i < x->val.array->length; i++)
  4395. list_push(r, together(state, v, x->val.array->data[i],
  4396. x->val.array->data[i - 1], 0, 0, v->rank[1],
  4397. v->rank[2]));
  4398. return value_new_array(r);
  4399. }
  4400. value_t *_adverb_eachprior_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4401. value_t *y) {
  4402. if (y->tag != ARRAY || !y->val.array->data)
  4403. return y;
  4404. value_t *_v = self->bonds->data[0];
  4405. if (_v->tag != VERB)
  4406. return state->udf;
  4407. verb_t *v = _v->val.verb;
  4408. list_t *r = list_new();
  4409. for (size_t i = 0; i < y->val.array->length; i++)
  4410. list_push(r, together(state, v, y->val.array->data[i],
  4411. i == 0 ? x : y->val.array->data[i - 1], 0, 0,
  4412. v->rank[1], v->rank[2]));
  4413. return value_new_array(r);
  4414. }
  4415. value_t *_adverb_reflex_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4416. value_t *_v = self->bonds->data[0];
  4417. if (_v->tag != VERB)
  4418. return state->udf;
  4419. verb_t *v = _v->val.verb;
  4420. return together(state, v, x, x, 0, 0, v->rank[1], v->rank[2]);
  4421. }
  4422. value_t *_adverb_reflex_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4423. value_t *y) {
  4424. value_t *_v = self->bonds->data[0];
  4425. if (_v->tag != VERB)
  4426. return state->udf;
  4427. verb_t *v = _v->val.verb;
  4428. return together(state, v, y, x, 0, 0, v->rank[1], v->rank[2]);
  4429. }
  4430. value_t *_adverb_amend_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4431. return state->udf;
  4432. }
  4433. value_t *_adverb_amend_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4434. value_t *y) {
  4435. if (x->tag != ARRAY)
  4436. x = verb_enlist(state, NULL, x);
  4437. value_t *v = self->bonds->data[0];
  4438. if (v->tag != ARRAY && v->tag != VERB)
  4439. v = verb_enlist(state, NULL, v);
  4440. if (y->tag != ARRAY)
  4441. y = verb_enlist(state, NULL, y);
  4442. list_t *r = list_copy(y->val.array);
  4443. size_t l = x->val.array->length;
  4444. if (v->tag == VERB) {
  4445. list_t *t = x->val.array;
  4446. for (size_t i = 0; i < t->length; i++) {
  4447. value_t *n = t->data[i];
  4448. if (n->tag != NUMBER)
  4449. break;
  4450. value_t *e = list_index(r, n->val.number);
  4451. if (!e)
  4452. continue;
  4453. list_set(r, n->val.number, apply_monad(state, v, e));
  4454. }
  4455. } else {
  4456. list_t *t = v->val.array;
  4457. for (size_t i = 0; i < t->length; i++) {
  4458. value_t *n = t->data[i];
  4459. if (n->tag != NUMBER)
  4460. break;
  4461. list_set(r, n->val.number, list_index(x->val.array, i < l ? i : l - 1));
  4462. }
  4463. }
  4464. return value_new_array(r);
  4465. }
  4466. value_t *_adverb_filter_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4467. value_t *_v = self->bonds->data[0];
  4468. if (_v->tag != VERB)
  4469. return state->udf;
  4470. if (x->tag != ARRAY)
  4471. x = verb_enlist(state, NULL, x);
  4472. else if (!x->val.array->data)
  4473. return x;
  4474. verb_t *v = _v->val.verb;
  4475. list_t *r = list_new();
  4476. for (size_t i = 0; i < x->val.array->length; i++) {
  4477. value_t *b = each_rank(state, v, x->val.array->data[i], 0, v->rank[0]);
  4478. if (value_is_truthy(b))
  4479. list_push(r, x->val.array->data[i]);
  4480. }
  4481. return value_new_array(r);
  4482. }
  4483. value_t *_adverb_filter_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4484. value_t *y) {
  4485. return state->udf;
  4486. }
  4487. value_t *_adverb_span_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4488. value_t *v = self->bonds->data[0];
  4489. if (v->tag != VERB)
  4490. return state->udf;
  4491. if (x->tag != ARRAY)
  4492. x = verb_enlist(state, NULL, x);
  4493. else if (!x->val.array->data)
  4494. return x;
  4495. list_t *r = list_new();
  4496. list_t *p = list_new();
  4497. for (size_t i = 0; i < x->val.array->length; i++) {
  4498. value_t *b = apply_monad(state, v, x->val.array->data[i]);
  4499. if (value_is_truthy(b)) {
  4500. list_push(r, value_new_array(p));
  4501. p = list_new();
  4502. } else
  4503. list_push(p, x->val.array->data[i]);
  4504. }
  4505. list_push(r, value_new_array(p));
  4506. return value_new_array(r);
  4507. }
  4508. value_t *_adverb_span_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4509. value_t *y) {
  4510. value_t *_v = self->bonds->data[0];
  4511. if (_v->tag != VERB)
  4512. return state->udf;
  4513. verb_t *v = _v->val.verb;
  4514. value_t *r = verb_windows(state, NULL, x, y);
  4515. return each_rank(state, v, r, 0, 1);
  4516. }
  4517. value_t *_adverb_inverse_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4518. value_t *_v = self->bonds->data[0];
  4519. if (_v->tag != VERB)
  4520. return state->udf;
  4521. verb_t *v = _v->val.verb;
  4522. verb_t *iv = table_get(Inverses, v->name);
  4523. if (!iv)
  4524. return state->udf;
  4525. return each_rank(state, iv, x, 0, iv->rank[0]);
  4526. }
  4527. value_t *_adverb_inverse_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4528. value_t *y) {
  4529. value_t *_v = self->bonds->data[0];
  4530. if (_v->tag != VERB)
  4531. return state->udf;
  4532. verb_t *v = _v->val.verb;
  4533. verb_t *iv = table_get(Inverses, v->name);
  4534. if (!iv)
  4535. return state->udf;
  4536. value_t *a = each_rank(state, iv, x, 0, iv->rank[0]);
  4537. value_t *b = each_rank(state, iv, y, 0, iv->rank[0]);
  4538. return apply_dyad(state, _v, a, b);
  4539. }
  4540. #define ADVERB(__name, __symb) \
  4541. verb_t *adverb_##__name(interpreter_t *state, value_t *v) { \
  4542. verb_t *nv = verb_new(); \
  4543. nv->bonds = list_newk(1); \
  4544. nv->bonds->data[0] = v; \
  4545. char *r = value_show(v); \
  4546. size_t l = strlen(r) + strlen(__symb) + 1; \
  4547. nv->name = malloc_checked_atomic(l); \
  4548. snprintf(nv->name, l, "%s" __symb, r); \
  4549. GC_FREE(r); \
  4550. nv->rank[0] = 0; \
  4551. nv->monad = _adverb_##__name##_monad; \
  4552. nv->dyad = _adverb_##__name##_dyad; \
  4553. return nv; \
  4554. }
  4555. ADVERB(fold, "/");
  4556. ADVERB(converge, "/.");
  4557. ADVERB(scan, "\\");
  4558. ADVERB(converges, "\\.");
  4559. ADVERB(each, "\"");
  4560. ADVERB(eachprior, "\".");
  4561. ADVERB(reflex, ";.");
  4562. ADVERB(amend, "`");
  4563. ADVERB(filter, "&.");
  4564. ADVERB(span, "/:");
  4565. ADVERB(inverse, "-:");
  4566. adverb_t ADVERBS[] = {
  4567. {"/", adverb_fold, NULL}, {"/.", adverb_converge, NULL},
  4568. {"\\", adverb_scan, NULL}, {"\\.", adverb_converges, NULL},
  4569. {"\"", adverb_each, NULL}, {"\".", adverb_eachprior, NULL},
  4570. {";.", adverb_reflex, NULL}, {"`", adverb_amend, NULL},
  4571. {"&.", adverb_filter, NULL}, {"/:", adverb_span, NULL},
  4572. {"-:", adverb_inverse, NULL}};
  4573. value_t *_conjunction_bond_monad(interpreter_t *state, verb_t *self,
  4574. value_t *x) {
  4575. value_t *v1 = self->bonds->data[0];
  4576. value_t *v2 = self->bonds->data[1];
  4577. if (v1->tag == VERB && v2->tag == VERB)
  4578. return apply_monad(state, v1, apply_monad(state, v2, x));
  4579. else if (v1->tag == VERB)
  4580. return apply_dyad(state, v1, x, v2);
  4581. else if (v2->tag == VERB)
  4582. return apply_dyad(state, v2, v1, x);
  4583. else
  4584. return state->udf;
  4585. }
  4586. value_t *_conjunction_bond_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4587. value_t *y) {
  4588. value_t *v1 = self->bonds->data[0];
  4589. value_t *v2 = self->bonds->data[1];
  4590. if (v1->tag == VERB && v2->tag == VERB)
  4591. return apply_monad(state, v1, apply_dyad(state, v2, x, y));
  4592. else if (v1->tag == VERB)
  4593. return apply_dyad(state, v1, apply_dyad(state, v1, x, y), v2);
  4594. else if (v2->tag == VERB)
  4595. return apply_dyad(state, v2, v1, apply_dyad(state, v2, x, y));
  4596. else
  4597. return state->udf;
  4598. }
  4599. value_t *_conjunction_pick_monad(interpreter_t *state, verb_t *self,
  4600. value_t *x) {
  4601. value_t *v1 = self->bonds->data[0];
  4602. value_t *v2 = self->bonds->data[1];
  4603. if (v1->tag != VERB || v2->tag != ARRAY)
  4604. return state->udf;
  4605. value_t *n = apply_monad(state, v1, x);
  4606. value_t *f = verb_at(state, NULL, v2, n);
  4607. return apply_monad(state, f, x);
  4608. }
  4609. value_t *_conjunction_pick_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4610. value_t *y) {
  4611. value_t *v1 = self->bonds->data[0];
  4612. value_t *v2 = self->bonds->data[1];
  4613. if (v1->tag != VERB || v2->tag != ARRAY)
  4614. return state->udf;
  4615. value_t *n = apply_dyad(state, v1, x, y);
  4616. value_t *f = verb_at(state, NULL, v2, n);
  4617. return apply_dyad(state, f, x, y);
  4618. }
  4619. value_t *_conjunction_while_monad(interpreter_t *state, verb_t *self,
  4620. value_t *x) {
  4621. value_t *v1 = self->bonds->data[0];
  4622. value_t *v2 = self->bonds->data[1];
  4623. if (v1->tag == VERB) {
  4624. for (;;) {
  4625. if (!value_is_truthy(apply_monad(state, v1, x)))
  4626. break;
  4627. x = apply_monad(state, v2, x);
  4628. }
  4629. } else if (v1->tag == NUMBER) {
  4630. size_t k = (size_t)fabs(v1->val.number);
  4631. for (size_t i = 0; i < k; i++)
  4632. x = apply_monad(state, v2, x);
  4633. }
  4634. return x;
  4635. }
  4636. value_t *_conjunction_while_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4637. value_t *y) {
  4638. value_t *v1 = self->bonds->data[0];
  4639. value_t *v2 = self->bonds->data[1];
  4640. if (v1->tag == VERB) {
  4641. for (;;) {
  4642. if (!value_is_truthy(apply_dyad(state, v1, x, y)))
  4643. break;
  4644. x = apply_dyad(state, v2, x, y);
  4645. }
  4646. } else if (v1->tag == NUMBER) {
  4647. size_t k = (size_t)fabs(v1->val.number);
  4648. for (size_t i = 0; i < k; i++)
  4649. x = apply_dyad(state, v2, x, y);
  4650. }
  4651. return x;
  4652. }
  4653. value_t *_conjunction_rank_monad(interpreter_t *state, verb_t *self,
  4654. value_t *x) {
  4655. value_t *v1 = self->bonds->data[0];
  4656. value_t *v2 = self->bonds->data[1];
  4657. if (v1->tag != VERB || v2->tag != NUMBER)
  4658. return state->udf;
  4659. unsigned int rank =
  4660. v2->val.number == INFINITY ? UINT_MAX : fabs(v2->val.number);
  4661. return each_rank(state, v1->val.verb, x, 0, rank);
  4662. }
  4663. value_t *_conjunction_rank_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4664. value_t *y) {
  4665. value_t *v1 = self->bonds->data[0];
  4666. value_t *v2 = self->bonds->data[1];
  4667. if (v1->tag != VERB)
  4668. return state->udf;
  4669. unsigned int rl;
  4670. unsigned int rr;
  4671. if (v2->tag == NUMBER)
  4672. rl = rr = v2->val.number == INFINITY ? UINT_MAX : fabs(v2->val.number);
  4673. else if (v2->tag == ARRAY && v2->val.array->length == 2) {
  4674. value_t *a = v2->val.array->data[0];
  4675. value_t *b = v2->val.array->data[1];
  4676. if (a->tag != NUMBER)
  4677. return state->udf;
  4678. rl = a->val.number == INFINITY ? UINT_MAX : fabs(a->val.number);
  4679. if (b->tag != NUMBER)
  4680. return state->udf;
  4681. rr = b->val.number == INFINITY ? UINT_MAX : fabs(b->val.number);
  4682. } else
  4683. return state->udf;
  4684. return together(state, v1->val.verb, x, y, 0, 0, rl, rr);
  4685. }
  4686. value_t *_conjunction_monaddyad_monad(interpreter_t *state, verb_t *self,
  4687. value_t *x) {
  4688. value_t *v = self->bonds->data[0];
  4689. if (v->tag != VERB)
  4690. return state->udf;
  4691. return each_rank(state, v->val.verb, x, 0, v->val.verb->rank[0]);
  4692. }
  4693. value_t *_conjunction_monaddyad_dyad(interpreter_t *state, verb_t *self,
  4694. value_t *x, value_t *y) {
  4695. value_t *v = self->bonds->data[1];
  4696. if (v->tag != VERB)
  4697. return state->udf;
  4698. return together(state, v->val.verb, x, y, 0, 0, v->val.verb->rank[1],
  4699. v->val.verb->rank[2]);
  4700. }
  4701. value_t *_conjunction_if_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4702. value_t *v1 = self->bonds->data[0];
  4703. value_t *v2 = self->bonds->data[1];
  4704. if (v1->tag != VERB || v2->tag != VERB)
  4705. return state->udf;
  4706. value_t *b = apply_monad(state, v2, x);
  4707. if (value_is_truthy(b))
  4708. return x;
  4709. return apply_monad(state, v1, x);
  4710. }
  4711. value_t *_conjunction_if_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4712. value_t *y) {
  4713. value_t *v1 = self->bonds->data[0];
  4714. value_t *v2 = self->bonds->data[1];
  4715. if (v1->tag != VERB || v2->tag != VERB)
  4716. return state->udf;
  4717. value_t *b = apply_dyad(state, v2, x, y);
  4718. if (value_is_truthy(b))
  4719. return y;
  4720. return apply_dyad(state, v1, x, y);
  4721. }
  4722. value_t *_conjunction_under_monad(interpreter_t *state, verb_t *self,
  4723. value_t *x) {
  4724. value_t *v1 = self->bonds->data[0];
  4725. value_t *v2 = self->bonds->data[1];
  4726. if (v1->tag != VERB || v2->tag != VERB)
  4727. return state->udf;
  4728. verb_t *iv = table_get(Inverses, v2->val.verb->name);
  4729. if (!iv)
  4730. return state->udf;
  4731. value_t *v = apply_monad(state, v2, x);
  4732. v = apply_monad(state, v1, v);
  4733. return each_rank(state, iv, v, 0, iv->rank[0]);
  4734. }
  4735. value_t *_conjunction_under_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4736. value_t *y) {
  4737. value_t *v1 = self->bonds->data[0];
  4738. value_t *v2 = self->bonds->data[1];
  4739. if (v1->tag != VERB || v2->tag != VERB)
  4740. return state->udf;
  4741. verb_t *iv = table_get(Inverses, v2->val.verb->name);
  4742. if (!iv)
  4743. return state->udf;
  4744. value_t *a = apply_monad(state, v2, x);
  4745. value_t *b = apply_monad(state, v2, y);
  4746. value_t *v = apply_dyad(state, v1, a, b);
  4747. return each_rank(state, iv, v, 0, iv->rank[0]);
  4748. }
  4749. value_t *_conjunction_collect_monad(interpreter_t *state, verb_t *self,
  4750. value_t *x) {
  4751. value_t *v1 = self->bonds->data[0];
  4752. value_t *v2 = self->bonds->data[1];
  4753. list_t *r = list_new();
  4754. if (v1->tag == VERB) {
  4755. for (;;) {
  4756. if (!value_is_truthy(apply_monad(state, v1, x)))
  4757. break;
  4758. list_push(r, x);
  4759. x = apply_monad(state, v2, x);
  4760. }
  4761. } else if (v1->tag == NUMBER) {
  4762. size_t k = (size_t)fabs(v1->val.number);
  4763. for (size_t i = 0; i < k; i++) {
  4764. list_push(r, x);
  4765. x = apply_monad(state, v2, x);
  4766. }
  4767. }
  4768. return value_new_array(r);
  4769. }
  4770. value_t *_conjunction_collect_dyad(interpreter_t *state, verb_t *self,
  4771. value_t *x, value_t *y) {
  4772. value_t *v1 = self->bonds->data[0];
  4773. value_t *v2 = self->bonds->data[1];
  4774. list_t *r = list_new();
  4775. if (v1->tag == VERB) {
  4776. for (;;) {
  4777. if (!value_is_truthy(apply_dyad(state, v1, x, y)))
  4778. break;
  4779. list_push(r, x);
  4780. x = apply_dyad(state, v2, x, y);
  4781. }
  4782. } else if (v1->tag == NUMBER) {
  4783. size_t k = (size_t)fabs(v1->val.number);
  4784. for (size_t i = 0; i < k; i++) {
  4785. list_push(r, x);
  4786. x = apply_dyad(state, v2, x, y);
  4787. }
  4788. }
  4789. return x;
  4790. }
  4791. #define CONJUNCTION(__name, __symb) \
  4792. verb_t *conjunction_##__name(interpreter_t *state, value_t *x, value_t *y) { \
  4793. verb_t *nv = verb_new(); \
  4794. nv->bonds = list_newk(2); \
  4795. nv->bonds->data[0] = x; \
  4796. nv->bonds->data[1] = y; \
  4797. char *rx = value_show(x); \
  4798. char *ry = value_show(y); \
  4799. size_t l = strlen(rx) + strlen(ry) + strlen(__symb) + 1; \
  4800. nv->name = malloc_checked_atomic(l); \
  4801. snprintf(nv->name, l, "%s" __symb "%s", rx, ry); \
  4802. GC_FREE(rx); \
  4803. GC_FREE(ry); \
  4804. nv->rank[0] = 0; \
  4805. nv->rank[1] = 0; \
  4806. nv->rank[1] = 0; \
  4807. nv->monad = _conjunction_##__name##_monad; \
  4808. nv->dyad = _conjunction_##__name##_dyad; \
  4809. return nv; \
  4810. }
  4811. CONJUNCTION(bond, ";");
  4812. CONJUNCTION(pick, "?.");
  4813. CONJUNCTION(while, "?:");
  4814. CONJUNCTION(rank, "\":");
  4815. CONJUNCTION(monaddyad, ";:");
  4816. CONJUNCTION(if, "&:");
  4817. CONJUNCTION(under, "^:");
  4818. CONJUNCTION(collect, "\\:");
  4819. adverb_t CONJUNCTIONS[] = {
  4820. {";", NULL, conjunction_bond}, {"?.", NULL, conjunction_pick},
  4821. {"?:", NULL, conjunction_while}, {"\":", NULL, conjunction_rank},
  4822. {";:", NULL, conjunction_monaddyad}, {"&:", NULL, conjunction_if},
  4823. {"^:", NULL, conjunction_under}, {"\\:", NULL, conjunction_collect}};
  4824. #define countof(x) (sizeof(x) / sizeof((x)[0]))
  4825. #define FINDER(kind, rname, table) \
  4826. kind *find_##rname(char *s) { \
  4827. for (size_t i = 0; i < countof(table); i++) { \
  4828. if (strcmp(table[i].name, s) == 0) \
  4829. return &table[i]; \
  4830. } \
  4831. return NULL; \
  4832. }
  4833. FINDER(verb_t, verb, VERBS);
  4834. FINDER(adverb_t, adverb, ADVERBS);
  4835. FINDER(adverb_t, conjunction, CONJUNCTIONS);
  4836. node_t *node_new(enum node_tag_t tag) {
  4837. node_t *node = malloc_checked(sizeof(node_t));
  4838. node->tag = tag;
  4839. return node;
  4840. }
  4841. node_t *node_new_strand(list_t *l) {
  4842. node_t *node = malloc_checked(sizeof(node_t));
  4843. node->tag = N_STRAND;
  4844. node->l = l;
  4845. return node;
  4846. }
  4847. node_t *node_new_literal(value_t *v) {
  4848. node_t *node = malloc_checked(sizeof(node_t));
  4849. node->tag = N_LITERAL;
  4850. node->v = v;
  4851. return node;
  4852. }
  4853. node_t *node_new1(enum node_tag_t tag, node_t *a) {
  4854. node_t *node = malloc_checked(sizeof(node_t));
  4855. node->tag = tag;
  4856. node->a = a;
  4857. return node;
  4858. }
  4859. node_t *node_new2(enum node_tag_t tag, node_t *a, node_t *b) {
  4860. node_t *node = malloc_checked(sizeof(node_t));
  4861. node->tag = tag;
  4862. node->a = a;
  4863. node->b = b;
  4864. return node;
  4865. }
  4866. node_t *node_new3(enum node_tag_t tag, node_t *a, node_t *b, node_t *c) {
  4867. node_t *node = malloc_checked(sizeof(node_t));
  4868. node->tag = tag;
  4869. node->a = a;
  4870. node->b = b;
  4871. node->c = c;
  4872. return node;
  4873. }
  4874. typedef struct {
  4875. lexer_t *lexer;
  4876. interpreter_t *state;
  4877. size_t pos;
  4878. size_t end;
  4879. size_t dp;
  4880. bool bn;
  4881. } parser_t;
  4882. parser_t *parser_new(interpreter_t *state) {
  4883. parser_t *parser = malloc_checked(sizeof(parser_t));
  4884. parser->state = state;
  4885. return parser;
  4886. }
  4887. void parser_error(parser_t *parser, char *s) { fatal(s); }
  4888. bool parser_done(parser_t *parser) { return parser->pos >= parser->end; }
  4889. token_t *parser_lookahead(parser_t *parser, size_t offset) {
  4890. size_t pos = parser->pos + offset;
  4891. if (pos >= parser->end)
  4892. return NULL;
  4893. return list_index(parser->lexer->tokens, pos);
  4894. }
  4895. bool parser_stop(parser_t *parser) {
  4896. token_t *tok = parser_lookahead(parser, 0);
  4897. if (!tok)
  4898. return true;
  4899. return tok->tag == T_RPAR;
  4900. }
  4901. void parser_eat(parser_t *parser) {
  4902. if (!parser_done(parser))
  4903. parser->pos++;
  4904. }
  4905. node_t *parser_parse_expr(parser_t *parser);
  4906. node_t *parser_parse_verb(parser_t *parser) {
  4907. token_t *tok = parser_lookahead(parser, 0);
  4908. if (!tok || tok->tag != T_PUNCT)
  4909. return NULL;
  4910. verb_t *verb = find_verb(tok->text);
  4911. if (!verb)
  4912. return NULL;
  4913. return node_new_literal(value_new_verb(verb));
  4914. }
  4915. value_t *_adverb_wrapper_monad(interpreter_t *state, verb_t *self, value_t *x) {
  4916. adverb_t *av = self->bonds->data[0];
  4917. if (x->tag != VERB)
  4918. return state->udf;
  4919. return value_new_verb(av->adverb(state, x));
  4920. }
  4921. value_t *_adverb_wrapper_dyad(interpreter_t *state, verb_t *self, value_t *x,
  4922. value_t *y) {
  4923. adverb_t *av = self->bonds->data[0];
  4924. if (x->tag != VERB)
  4925. return state->udf;
  4926. verb_t *v = av->adverb(state, x);
  4927. return each_rank(state, v, y, 0, v->rank[0]);
  4928. }
  4929. node_t *parser_parse_adverb_atom(parser_t *parser) {
  4930. token_t *tok = parser_lookahead(parser, 0);
  4931. if (!tok || tok->tag != T_PUNCT)
  4932. return NULL;
  4933. adverb_t *adverb = find_adverb(tok->text);
  4934. if (!adverb)
  4935. return NULL;
  4936. verb_t *nv = verb_new();
  4937. nv->name = strdup_checked(tok->text);
  4938. nv->bonds = list_newk(1);
  4939. nv->bonds->data[0] = adverb;
  4940. nv->rank[0] = 0;
  4941. nv->rank[1] = 0;
  4942. nv->rank[2] = 0;
  4943. nv->monad = _adverb_wrapper_monad;
  4944. nv->dyad = _adverb_wrapper_dyad;
  4945. return node_new_literal(value_new_verb(nv));
  4946. }
  4947. value_t *_conjunction_wrapper_dyad(interpreter_t *state, verb_t *self,
  4948. value_t *x, value_t *y) {
  4949. adverb_t *av = self->bonds->data[0];
  4950. return value_new_verb(av->conjunction(state, x, y));
  4951. }
  4952. node_t *parser_parse_conjunction_atom(parser_t *parser) {
  4953. token_t *tok = parser_lookahead(parser, 0);
  4954. if (!tok || tok->tag != T_PUNCT)
  4955. return NULL;
  4956. adverb_t *adverb = find_conjunction(tok->text);
  4957. if (!adverb)
  4958. return NULL;
  4959. verb_t *nv = verb_new();
  4960. nv->name = strdup_checked(tok->text);
  4961. nv->bonds = list_newk(1);
  4962. nv->bonds->data[0] = adverb;
  4963. nv->rank[0] = 0;
  4964. nv->rank[1] = 0;
  4965. nv->rank[2] = 0;
  4966. nv->monad = NULL;
  4967. nv->dyad = _conjunction_wrapper_dyad;
  4968. return node_new_literal(value_new_verb(nv));
  4969. }
  4970. node_t *parser_parse_atom(parser_t *parser) {
  4971. token_t *tok = parser_lookahead(parser, 0);
  4972. node_t *node = NULL;
  4973. switch (tok->tag) {
  4974. case T_RPAR:
  4975. parser_error(parser, "unmatched");
  4976. case T_LPAR:
  4977. parser_eat(parser);
  4978. tok = parser_lookahead(parser, 0);
  4979. if (tok && tok->tag == T_RPAR) {
  4980. node = node_new_literal(parser->state->unit);
  4981. break;
  4982. }
  4983. parser->dp++;
  4984. node = parser_parse_expr(parser);
  4985. if (!node)
  4986. parser_error(parser, "unmatched");
  4987. if (parser->bn)
  4988. node->dp = 2;
  4989. else
  4990. node->dp = parser->dp;
  4991. parser->dp--;
  4992. tok = parser_lookahead(parser, 0);
  4993. if (!tok || tok->tag != T_RPAR)
  4994. parser_error(parser, "unmatched");
  4995. break;
  4996. case T_PUNCT:
  4997. node = parser_parse_verb(parser);
  4998. if (!node)
  4999. node = parser_parse_adverb_atom(parser);
  5000. if (!node)
  5001. node = parser_parse_conjunction_atom(parser);
  5002. if (!node)
  5003. parser_error(parser, "parse");
  5004. break;
  5005. case T_NUMBER:
  5006. node = node_new_literal(value_new_number(strtod(tok->text, NULL)));
  5007. break;
  5008. case T_BNUMBER: {
  5009. if (!tok->text[1])
  5010. parser_error(parser, "trailing-base");
  5011. int base = tok->text[0] == 'x' ? 16 : tok->text[0] == 'b' ? 2 : 8;
  5012. node =
  5013. node_new_literal(value_new_number(strtol(tok->text + 1, NULL, base)));
  5014. } break;
  5015. case T_NAME:
  5016. node = node_new_literal(value_new_symbol(strdup_checked(tok->text)));
  5017. break;
  5018. case T_QUOTE:
  5019. if (!*tok->text)
  5020. node = node_new_literal(parser->state->unit);
  5021. else if (!*(tok->text + 1))
  5022. node = node_new_literal(value_new_char(tok->text[0]));
  5023. else {
  5024. size_t z = strlen(tok->text);
  5025. list_t *r = list_newk(z);
  5026. for (size_t i = 0; i < z; i++)
  5027. r->data[i] = CHARS[(int)tok->text[i]];
  5028. node = node_new_literal(value_new_array(r));
  5029. }
  5030. break;
  5031. }
  5032. if (!node)
  5033. parser_error(parser, "parse");
  5034. parser_eat(parser);
  5035. return node;
  5036. }
  5037. bool is_unbound(interpreter_t *state, char *s) {
  5038. if (state->args->data) {
  5039. list_t *args = list_index(state->args, -1);
  5040. size_t argc = args->length - 1;
  5041. if (argc == 2 && strcmp(s, "y") == 0)
  5042. return false;
  5043. else if (strcmp(s, "x") == 0)
  5044. return false;
  5045. } else if (table_has(state->env, s))
  5046. return false;
  5047. return true;
  5048. }
  5049. node_t *parser_parse_sequence(parser_t *parser, node_t *a,
  5050. enum token_tag_t tag) {
  5051. token_t *tok;
  5052. if ((tok = parser_lookahead(parser, 0)) &&
  5053. (tok->tag == tag || (tag == T_NUMBER && tok->tag == T_BNUMBER))) {
  5054. if (tag == T_NAME && !is_unbound(parser->state, tok->text))
  5055. return NULL;
  5056. list_t *as = list_new();
  5057. list_push(as, a->v);
  5058. do {
  5059. if (tag == T_NAME && tok->tag == T_NAME &&
  5060. !is_unbound(parser->state, tok->text))
  5061. break;
  5062. a = parser_parse_atom(parser);
  5063. list_push(as, a->v);
  5064. } while ((tok = parser_lookahead(parser, 0)) &&
  5065. (tok->tag == tag || (tag == T_NUMBER && tok->tag == T_BNUMBER)));
  5066. return node_new_literal(value_new_array(as));
  5067. }
  5068. return NULL;
  5069. }
  5070. node_t *_parser_parse_noun(parser_t *parser) {
  5071. node_t *n;
  5072. node_t *a = parser_parse_atom(parser);
  5073. if (a->tag == N_LITERAL && a->v->tag == NUMBER &&
  5074. (n = parser_parse_sequence(parser, a, T_NUMBER)))
  5075. return n;
  5076. else if (a->tag == N_LITERAL && a->v->tag == SYMBOL &&
  5077. is_unbound(parser->state, a->v->val.symbol) &&
  5078. (n = parser_parse_sequence(parser, a, T_NAME)))
  5079. return n;
  5080. else if (a->tag == N_LITERAL &&
  5081. ((a->v->tag == ARRAY && is_char_array(a->v->val.array)) ||
  5082. a->v->tag == CHAR) &&
  5083. (n = parser_parse_sequence(parser, a, T_QUOTE)))
  5084. return n;
  5085. return a;
  5086. }
  5087. node_t *parser_parse_noun(parser_t *parser, bool flat) {
  5088. node_t *a = flat ? parser_parse_atom(parser) : _parser_parse_noun(parser);
  5089. token_t *tok;
  5090. if ((tok = parser_lookahead(parser, 0)) && tok->tag == T_PUNCT &&
  5091. strcmp(tok->text, ",:") == 0) {
  5092. parser_eat(parser);
  5093. list_t *l = list_new();
  5094. list_push(l, a);
  5095. for (;;) {
  5096. if (parser_stop(parser))
  5097. parser_error(parser, "trailing-strand");
  5098. a = flat ? parser_parse_atom(parser) : _parser_parse_noun(parser);
  5099. list_push(l, a);
  5100. if (!((tok = parser_lookahead(parser, 0)) && tok->tag == T_PUNCT &&
  5101. strcmp(tok->text, ",:") == 0))
  5102. break;
  5103. parser_eat(parser);
  5104. }
  5105. return node_new_strand(l);
  5106. }
  5107. return a;
  5108. }
  5109. bool parser_node_is_verbal(parser_t *parser, node_t *n) {
  5110. value_t *v;
  5111. if (n->tag == N_FUN)
  5112. return true;
  5113. else if (n->tag == N_ADV || n->tag == N_CONJ || n->tag == N_PARTIAL_CONJ)
  5114. return true;
  5115. else if (n->tag == N_FORK || n->tag == N_HOOK || n->tag == N_BOND ||
  5116. n->tag == N_OVER)
  5117. return true;
  5118. else if (n->tag == N_LITERAL && n->v->tag == VERB)
  5119. return true;
  5120. else if (n->tag == N_LITERAL && n->v->tag == SYMBOL &&
  5121. (v = table_get(parser->state->env, n->v->val.symbol)) &&
  5122. v->tag == VERB)
  5123. return true;
  5124. return false;
  5125. }
  5126. node_t *parser_parse_adverb(parser_t *parser, node_t *v, bool *flag) {
  5127. token_t *tok;
  5128. adverb_t *adv;
  5129. node_t *t;
  5130. for (;;) {
  5131. tok = parser_lookahead(parser, 0);
  5132. if (!tok || tok->tag != T_PUNCT)
  5133. break;
  5134. if ((adv = find_adverb(tok->text))) {
  5135. if (flag)
  5136. *flag = true;
  5137. parser_eat(parser);
  5138. t = node_new(N_ADV);
  5139. t->av = adv;
  5140. t->a = v;
  5141. v = t;
  5142. } else
  5143. break;
  5144. }
  5145. return v;
  5146. }
  5147. node_t *parser_parse_conjunction(parser_t *parser, node_t *v, bool *flag) {
  5148. token_t *tok;
  5149. adverb_t *adv;
  5150. node_t *t;
  5151. for (;;) {
  5152. tok = parser_lookahead(parser, 0);
  5153. if (!tok || tok->tag != T_PUNCT)
  5154. break;
  5155. if ((adv = find_conjunction(tok->text))) {
  5156. if (flag)
  5157. *flag = true;
  5158. parser_eat(parser);
  5159. if (parser_stop(parser)) {
  5160. t = node_new(N_PARTIAL_CONJ);
  5161. t->av = adv;
  5162. t->a = v;
  5163. } else {
  5164. t = node_new(N_CONJ);
  5165. t->av = adv;
  5166. t->a = v;
  5167. t->b = parser_parse_noun(parser, true);
  5168. }
  5169. v = t;
  5170. } else
  5171. break;
  5172. }
  5173. return v;
  5174. }
  5175. bool is_apply(node_t *n) {
  5176. return n->tag == N_LITERAL && n->v->tag == VERB &&
  5177. (strcmp(n->v->val.verb->name, "`.") == 0 ||
  5178. strcmp(n->v->val.verb->name, "`:") == 0);
  5179. }
  5180. bool is_obverse(node_t *n) {
  5181. return n->tag == N_LITERAL && n->v->tag == VERB &&
  5182. strcmp(n->v->val.verb->name, "::") == 0;
  5183. }
  5184. node_t *parser_parse_expr(parser_t *parser) {
  5185. token_t *tmp;
  5186. list_t *ns = list_new();
  5187. while (!parser_stop(parser)) {
  5188. if (!ns->data && (tmp = parser_lookahead(parser, 0)) &&
  5189. tmp->tag == T_PUNCT && strcmp(tmp->text, ":") == 0) {
  5190. parser_eat(parser);
  5191. node_t *r = parser_parse_expr(parser);
  5192. if (!r)
  5193. r = node_new_literal(parser->state->nil);
  5194. return node_new1(N_FUN, r);
  5195. }
  5196. node_t *n = parser_parse_noun(parser, false);
  5197. if (!ns->data && n->tag == N_LITERAL && n->v->tag == SYMBOL &&
  5198. (tmp = parser_lookahead(parser, 0)) && tmp->tag == T_PUNCT &&
  5199. strcmp(tmp->text, ":") == 0) {
  5200. parser_eat(parser);
  5201. bool t = parser->bn;
  5202. parser->bn = true;
  5203. node_t *r = parser_parse_expr(parser);
  5204. parser->bn = t;
  5205. return node_new2(N_BIND, n, r);
  5206. }
  5207. for (;;) {
  5208. bool flag = false;
  5209. n = parser_parse_adverb(parser, n, &flag);
  5210. n = parser_parse_conjunction(parser, n, &flag);
  5211. if (!flag)
  5212. break;
  5213. }
  5214. list_push(ns, n);
  5215. }
  5216. size_t len;
  5217. node_t *l, *m, *r;
  5218. for (;;) {
  5219. len = ns->length;
  5220. if (len < 2)
  5221. break;
  5222. if (len >= 3 &&
  5223. (is_apply(list_index(ns, -2)) || is_obverse(list_index(ns, -2))) &&
  5224. parser_node_is_verbal(parser, list_index(ns, -3))) {
  5225. r = list_pop(ns);
  5226. m = list_pop(ns);
  5227. l = list_pop(ns);
  5228. list_push(ns, node_new3(N_DYAD, m, l, r));
  5229. } else if (len >= 3 && !parser_node_is_verbal(parser, list_index(ns, -1)) &&
  5230. parser_node_is_verbal(parser, list_index(ns, -2)) &&
  5231. !parser_node_is_verbal(parser, list_index(ns, -3))) {
  5232. r = list_pop(ns);
  5233. m = list_pop(ns);
  5234. l = list_pop(ns);
  5235. list_push(ns, node_new3(N_DYAD, m, l, r));
  5236. } else if (len >= 3 && parser_node_is_verbal(parser, list_index(ns, -1)) &&
  5237. parser_node_is_verbal(parser, list_index(ns, -2)) &&
  5238. parser_node_is_verbal(parser, list_index(ns, -3))) {
  5239. r = list_pop(ns);
  5240. m = list_pop(ns);
  5241. l = list_pop(ns);
  5242. list_push(ns, node_new3(N_FORK, l, m, r));
  5243. } else if (len >= 3 && parser_node_is_verbal(parser, list_index(ns, -1)) &&
  5244. parser_node_is_verbal(parser, list_index(ns, -2)) &&
  5245. !parser_node_is_verbal(parser, list_index(ns, -3))) {
  5246. r = list_pop(ns);
  5247. m = list_pop(ns);
  5248. l = list_pop(ns);
  5249. list_push(ns, node_new3(N_OVER, l, m, r));
  5250. } else if (len >= 2 && is_apply(list_index(ns, -1))) {
  5251. r = list_pop(ns);
  5252. l = list_pop(ns);
  5253. list_push(ns, node_new2(N_BOND, r, l));
  5254. } else if (len >= 2 && !parser_node_is_verbal(parser, list_index(ns, -1)) &&
  5255. parser_node_is_verbal(parser, list_index(ns, -2))) {
  5256. r = list_pop(ns);
  5257. l = list_pop(ns);
  5258. list_push(ns, node_new2(N_MONAD, l, r));
  5259. } else if (len >= 2 && parser_node_is_verbal(parser, list_index(ns, -1)) &&
  5260. parser_node_is_verbal(parser, list_index(ns, -2))) {
  5261. r = list_pop(ns);
  5262. l = list_pop(ns);
  5263. list_push(ns, node_new2(N_HOOK, l, r));
  5264. } else if (len >= 2 && parser_node_is_verbal(parser, list_index(ns, -1)) &&
  5265. !parser_node_is_verbal(parser, list_index(ns, -2))) {
  5266. r = list_pop(ns);
  5267. l = list_pop(ns);
  5268. list_push(ns, node_new2(N_BOND, r, l));
  5269. } else if (len >= 3) {
  5270. r = list_pop(ns);
  5271. m = list_pop(ns);
  5272. l = list_pop(ns);
  5273. list_push(ns, node_new3(N_INDEX2, m, l, r));
  5274. } else if (len >= 2) {
  5275. r = list_pop(ns);
  5276. l = list_pop(ns);
  5277. list_push(ns, node_new2(N_INDEX1, l, r));
  5278. }
  5279. }
  5280. return ns->data ? ns->data[0] : NULL;
  5281. }
  5282. node_t *parser_parse(parser_t *parser, lexer_t *lexer) {
  5283. parser->lexer = lexer;
  5284. parser->pos = 0;
  5285. parser->end = parser->lexer->tokens->length;
  5286. node_t *node = parser_parse_expr(parser);
  5287. if (!parser_done(parser)) {
  5288. token_t *tok = parser_lookahead(parser, 0);
  5289. if (tok && tok->tag == T_RPAR)
  5290. parser_error(parser, "unmatched");
  5291. parser_error(parser, "parse");
  5292. }
  5293. return node;
  5294. }
  5295. value_t *interpreter_run(interpreter_t *state, char *program) {
  5296. lexer_t *lexer = lexer_new();
  5297. lexer_lex(lexer, program);
  5298. parser_t *parser = parser_new(state);
  5299. node_t *node = parser_parse(parser, lexer);
  5300. list_t *t = lexer->tokens;
  5301. for (size_t i = 0; i < t->length; i++) {
  5302. token_t *tok = t->data[i];
  5303. if (tok->text)
  5304. GC_FREE(tok->text);
  5305. GC_FREE(tok);
  5306. }
  5307. GC_FREE(t->data);
  5308. GC_FREE(t);
  5309. value_t *r = interpreter_walk(state, node);
  5310. GC_FREE(parser);
  5311. return r;
  5312. }
  5313. #include "help.h"
  5314. const char *VSTR = VER " " __DATE__;
  5315. void jkexec(interpreter_t *state, FILE *fd, bool isrepl, char **s) {
  5316. value_t *v = NULL;
  5317. list_t *r;
  5318. if (!isrepl)
  5319. r = list_new();
  5320. for (;;) {
  5321. buffer_t *buffer;
  5322. char line[256];
  5323. buffer = buffer_new();
  5324. if (isrepl)
  5325. putc('\t', stdout);
  5326. if (!fgets(line, sizeof(line), fd))
  5327. break;
  5328. if (isrepl) {
  5329. if (strcmp(line, "\\\\\n") == 0)
  5330. break;
  5331. else if (strcmp(line, "\\\n") == 0) {
  5332. printf("%s", HELP);
  5333. continue;
  5334. } else if (strcmp(line, "\\0\n") == 0) {
  5335. printf("%s", SHELP);
  5336. continue;
  5337. } else if (strcmp(line, "\\+\n") == 0) {
  5338. printf("%s", VHELP);
  5339. continue;
  5340. } else if (strcmp(line, "\\a\n") == 0) {
  5341. printf("%s", V2HELP);
  5342. continue;
  5343. } else if (strcmp(line, "\\\"\n") == 0) {
  5344. printf("%s", AHELP);
  5345. continue;
  5346. } else if (strcmp(line, "\\;\n") == 0) {
  5347. printf("%s", CHELP);
  5348. continue;
  5349. } else if (strcmp(line, "\\-:\n") == 0) {
  5350. printf("%s", IHELP);
  5351. continue;
  5352. }
  5353. }
  5354. while (strlen(line) > 2 && strcmp(line + strlen(line) - 3, "..\n") == 0) {
  5355. line[strlen(line) - 3] = 0;
  5356. buffer_append_str(buffer, line);
  5357. if (isrepl)
  5358. putc('\t', stdout);
  5359. if (!fgets(line, sizeof(line), fd))
  5360. return;
  5361. }
  5362. buffer_append_str(buffer, line);
  5363. *s = buffer_read(buffer);
  5364. v = interpreter_run(state, *s);
  5365. GC_FREE(*s);
  5366. *s = NULL;
  5367. if (isrepl && v->tag != NIL) {
  5368. table_set(state->env, "it", v);
  5369. char *s = value_show(v);
  5370. fputs(s, stdout);
  5371. GC_FREE(s);
  5372. if (isrepl)
  5373. putc('\n', stdout);
  5374. } else if (!isrepl && v && v->tag != NIL)
  5375. list_push(r, v);
  5376. }
  5377. if (!isrepl && r->data) {
  5378. char *s = value_show(list_index(r, -1));
  5379. fputs(s, stdout);
  5380. GC_FREE(s);
  5381. }
  5382. }
  5383. int main(int argc, char **argv) {
  5384. GC_INIT();
  5385. guards = list_new();
  5386. is_interactive = isatty(0);
  5387. HASH_SEED = time(NULL);
  5388. srand(HASH_SEED);
  5389. VCACHE = table_new();
  5390. SCACHE = table_new();
  5391. for (size_t i = 0; i < countof(VERBS); i++) {
  5392. value_t *v = value_new_const(VERB);
  5393. v->val.verb = &VERBS[i];
  5394. table_set(VCACHE, VERBS[i].name, v);
  5395. }
  5396. _UNIT = value_new(ARRAY);
  5397. _UNIT->val.array = list_new();
  5398. interpreter_t *state = interpreter_new();
  5399. for (int i = 1; i <= 8; i++) {
  5400. NNUMS[i - 1] = value_new_const(NUMBER);
  5401. NNUMS[i - 1]->val.number = -i;
  5402. }
  5403. for (int i = 0; i < 256; i++) {
  5404. NUMS[i] = value_new_const(NUMBER);
  5405. NUMS[i]->val.number = i;
  5406. }
  5407. list_t *cs = list_newk(256);
  5408. for (int i = 0; i < 256; i++) {
  5409. CHARS[i] = value_new_const(CHAR);
  5410. CHARS[i]->val._char = i;
  5411. cs->data[i] = CHARS[i];
  5412. }
  5413. _NAN = value_new_const(NUMBER);
  5414. _NAN->val.number = NAN;
  5415. INF = value_new_const(NUMBER);
  5416. INF->val.number = INFINITY;
  5417. NINF = value_new_const(NUMBER);
  5418. NINF->val.number = -INFINITY;
  5419. list_t *vs = list_new();
  5420. for (size_t i = 0; i < strlen(VSTR); i++)
  5421. list_push(vs, CHARS[(int)VSTR[i]]);
  5422. table_set(state->env, "A", value_new_array(cs));
  5423. table_set(state->env, "JKV", value_new_array(vs));
  5424. table_set(state->env, "E", value_new_number(exp(1)));
  5425. table_set(state->env, "pi", value_new_number(M_PI));
  5426. table_set(state->env, "tau", value_new_number(M_PI * 2));
  5427. table_set(state->env, "nan", _NAN);
  5428. table_set(state->env, "inf", INF);
  5429. table_set(state->env, "nil", state->nil);
  5430. table_set(state->env, "udf", state->udf);
  5431. Inverses = table_new();
  5432. table_set(Inverses, "+", find_verb("+"));
  5433. table_set(Inverses, "-", find_verb("-"));
  5434. table_set(Inverses, "|", find_verb("|"));
  5435. table_set(Inverses, "~", find_verb("~"));
  5436. table_set(Inverses, "%", find_verb("%"));
  5437. table_set(Inverses, "]", find_verb("]"));
  5438. table_set(Inverses, "*:", find_verb("%:"));
  5439. table_set(Inverses, "%:", find_verb("*:"));
  5440. table_set(Inverses, ">", find_verb("<"));
  5441. table_set(Inverses, "<", find_verb(">"));
  5442. table_set(Inverses, "_.", find_verb("_:"));
  5443. table_set(Inverses, "_:", find_verb("_."));
  5444. table_set(Inverses, "^.", find_verb("^"));
  5445. table_set(Inverses, "^", find_verb("^."));
  5446. table_set(Inverses, "+;.", find_verb("%:"));
  5447. table_set(Inverses, "*/", find_verb("["));
  5448. table_set(Inverses, "[", interpreter_run(state, "*/")->val.verb);
  5449. table_set(Inverses, "!", interpreter_run(state, ">|/")->val.verb);
  5450. table_set(Inverses, "!.", interpreter_run(state, "|/")->val.verb);
  5451. table_set(Inverses, "]@>:", interpreter_run(state, "]@<:")->val.verb);
  5452. table_set(Inverses, "]@<:", interpreter_run(state, "]@>:")->val.verb);
  5453. list_t *args = list_new();
  5454. for (int i = 1; i < argc; i++) {
  5455. list_t *arg = list_new();
  5456. char *s = argv[i];
  5457. while (*s)
  5458. list_push(arg, CHARS[(int)(*s++)]);
  5459. list_push(args, value_new_array(arg));
  5460. }
  5461. table_set(state->env, "args", value_new_array(args));
  5462. if (is_interactive)
  5463. printf("jk\t\\\\ to exit \\ for help\n");
  5464. char *s = NULL;
  5465. if (is_interactive)
  5466. setjmp(interactive_checkpoint);
  5467. if (s) {
  5468. GC_FREE(s);
  5469. s = NULL;
  5470. }
  5471. jkexec(state, stdin, is_interactive, &s);
  5472. }