k tutorial¶
Alexander Belopolsky and Dennis Shasha
date: December 04, 2019 (with frequent updates)
Since 1992, Arthur Whitney's k and its derivatives have served a small number of highly skilled programmers, enabling them to create high performance (and high value) applications for finance and other data-intensive domains. While implementation efforts in other languages such as C++ and Java (and to a lesser extent Python) often involve thousands of lines of code, much of it built on top of libraries, a typical k application is on the order of scores of lines of code without the need for libraries. The expressive power is in the language itself.
This tutorial aims to familiarize users of conventional programming languages with the latest iteration of k, Shakti k. The language is capable of managing streaming, in-memory, historical, relational, and time-series data. The distributed model extends out to multiple machines whether on-premise or in the cloud.
Shakti k provides connectivity via Python, HTTP, SSL/TLS, and json. Shakti k supports compression and encryption for data, whether in-memory, in-flight or on disk. Shakti k also has primitives for blockchain operations.
The tutorial introduces language concepts, then presents examples. A good way to learn the language is to try to program the examples on your own.
Readers are invited to suggest corrections or new examples. You can contact the authors at tutorial@shakti.com.
First encounter¶
Installing k¶
This can change, but as of March 2019, Shakti Software distributes a free evaluation version of k through Anaconda.org. The use of free version is subject to Evaluation Agreement.
Please follow the installation instructions at https://anaconda.org/shaktidb/shakti.
Shakti k does not have any dependencies and once you install it, you are ready to go. Simply type k at the command prompt and you will see the k banner and the prompt will change to a single space:
$ k
2019-09-25 14:44:44 20core 64gb avx2 © shakti m test
█
The banner starts with the timestamp corresponding to the modification time of the k program. The timestamp is followed by a letter that will be M if you are using macOS and L if you are using Linux. The Ncore and Mgb parts show how many CPU cores (N) your k session will use and how much memory (M gigabytes) is earmarked for k use.
Power tip: for better interactive experience we recommend installing
the rlwrap
utility and define an alias
alias k="rlwrap k"
this will allow editing the expression that you enter at the k prompt and to recall any previous input from history.
Using k as a calculator¶
You can start using k as a powerful calculator: enter an expression at the prompt, press Enter and k will evaluate the expression and print the result. Many avaliable operations will look familiar, but you will soon discover some features that are unique to k.
Arithmetics¶
In k, +
, -
, and *
work as the usual addition, subtraction and
multiplication operations, e.g.,
3*4
12
but the division operator is %
while /
has several uses including
serving as a prefix for comments that k will ignore:
10%3 / 10 divided by 3
3.333333
The next feature that may come as a surprise is that k does not use the traditional order of operations
3*2+4 / addition is performed first
18
Instead of the (P)EMDAS order, k consistently evaluates its expressions from right to left with only the parentheses having higher precedence
(3*2)+4 / multiplication is performed first
10
Elementary functions¶
As any good scientific calculator, k comes with a number of built-in functions. You can apply these functions by simpy typing their names before the argument separated by a space
sqrt 2
1.414214
Trigonometric functions operate on arguments in radians and you will often need the π constant to convert from degrees. The π constant is built in in k and if you are using k on a Mac, you can type it using the alt-p key combination
sin π%2
1f
Special values¶
Unlike some other languages that are quick to give up and report an error
when given invalid input, k tries hard to provide useful answers. Thus
if the result of a function is infinite, k will return a special ∞
value
and indicate the sign of the infinity and such invalid result may dissapear
in the subsequent computations
log 0
-∞
exp log 0
0f
When the result is completely undefined, k will return ø
, which stands for
missing data
(log 0) + 1 % 0
ø
Unary operators¶
In the previous sections, we have seen operators that take two numbers as
operands and functions that take one number as an argument. From elementary
math, you are familiar with a unary -
operator that takes a single operand
and returns its negation. Not surprisingly, -
does the same in k:
- 42
-42
However, k takes this idea of the same operator having both binary and unary
forms to the whole new level. Each operator in k has both unary and binary
forms. For example, similar to -
, the k division operator, %
, has a unary
form that computes the inverse
% 3
0.3333333
Whetting your appetite¶
As you experiment with k as a calculator, you will soon notice that it does not
complain about seemingly meaningless keystrokes. Thus a +
sign on its own or
following a number are simply echoed back by k
+
+
2+
2+
Entering the entire top row of shifted symbols is somehow understood by k as a valid expression
~!@#$%^&*()_+
~:!:@:#:$:%:^:&:*:()_+
We will explain what is happening here in future sections.
Reference card¶
Within a k session, you can type \h
to get
a summary of the basic operations.
\h
$k [-p 1234] [-v] [f.k] .z.i(pid) .z.x(arg) .z.e(env)
Verb Adverb Noun Type
: assign ' each n bar char " ab" `c
+ add flip / over n div name ``a`b `n
- subtract negate \ scan n mod int Ø 0 2 `i
* multiply first ': eachprior peach float ø 2.3 π ∞ `f
% divideby inverse /: eachright sv|join date 2019-06-28 `D .z.D
& min|and where \: eachleft vs|split time 12:34:56.789 `t .z.t
| max|or reverse
< less up I/O (0:i close) list (2;3.4;`c) `.
> more down 0: read/write line dict ``b!(2;`c) `a
= equal group 1: read/write byte tabl +{b:2 3 4} `A
~ match not 2: read/write data expr :32+9*f%5 `0
! key enum 3: set/conn (.Z.s) func {(+/x)%#x} `1
, catenate enlist 4: get/http (.Z.gGP)
^ except asc
# take(f#) count #[t;c;b[;a]] select \l f.k load \h help
_ drop(f_) floor _[t;c;b[;a]] update \t[:n]x milli/time \l log
? find unique ?[x;i;f[;y]] splice \u[:n]x micro/trace
@ index type @[x;i;f[;y]] amend \v [d] vars \f [d] fns
. apply value .[x;i;f[;y]] dmend \lf file \lc char \ll line
$ cast|pad string $[c;t;f] cond \cd [d] get[set] \gr x file
util: in within bin like find freq prm cmb
math: sqrt sin cos abs [n]log [n]exp [n]rand
aggr: count first last min max sum avg var dev med
tabl: [f]asc [f]dsc [f]key [select|update]A by B from T where C
datetime:YMDHRSTUV+duration:ymdhrstuv; .z.D+.z.t day:7\ mth:12\ bar:5'
`year`month`date`hour`minute`second`millisecond`microsecond`nanosecond
Ktab: (+`a!2 3)!+`b!3 4
Xtab: `a`b!+`c`d!(2 3;4 5)
1way: `p`m`crc`sha`rip`bad e.g. hash: `sha@"crypto"
2way:``j`k`csv`b64`hex`aes e.g. json: `j?`j@{a:2;b:`c}
if[c;..];while[c;..]
2+3 /comment
2+ \3 /display
\\ exit
trap: [e]f/:x
time: [n]f\:x
attr: ^x ?x
syns: count first last sum min max asc
Lists¶
Much of the expressiveness of k derives from the fact that most operations that operate on single values (atoms) generalize nicely to lists.
The simplest list is simply a sequence of numbers separated by spaces. When you apply one of the arithmetic operations between an atom and a list, k computes the result of the operation between each element of the list and the atom. For example,
1 2 3 4 * 10
10 20 30 40
You can also perform operations on the lists of the same length and k will pair up the corresponding elements and apply the operation to each pair.
1 2 3 + 3 2 1
4 4 4
However, if the lengths don't match, k will signal an error
1 2 3 + 3 2
1 2 3 + 3 2
^
length error
>
We will explain the meaning of error displays later, but for now
you just need to know that entering \
at the
>
prompt will clear the error and allow you to continue.
Generating lists¶
Entering long lists into k can soon become tedious and k provides
nice ways to generate lists either deterministically or randomly.
You can generate a uniform list of any length by placing &
, !
or rand
in
front of a number:
&N
- N zeros!N
- 0 through N-1rand N
- N random numbers drawn uniformly from [0, 1]rand -N
- N random numbers drawn from a normal distribution with μ=0 and σ=1
Zeros¶
Instead of typing 0
twelve times for a list of twelve zeros,
we can tell k to generate a list for us
&12
0 0 0 0 0 0 0 0 0 0 0 0
Index¶
x: !10
x / this is the enumeration of 0 1 2 3 4 5 6 7 8 9 (but summarized):
!10
5 + x
5+!10
Random lists¶
10 rand 50 / generate randomly with replacement (so there can be duplicates)
45 13 28 28 8 43 9 17 30 7
(because of randomness, your result might not exactly match the above)
10 rand 10
8 6 8 2 5 2 0 1 0 2
-10 rand 10 / generate randomly without replacement (no duplicates)
9 7 0 1 5 3 6 4 2 8
There are more advanced ways to generate random arrays:
20 rand 100 / recall uniform random with replacement
11 0 55 17 36 50 73 85 62 93 16 58 75 81 81 72 36 90 21 78
rand 6 / uniform random with replacement between 0 and 1
0.3937393 0.4717788 0.2755357 0.1641728 0.9861826 0.8855961
/ Others are on the way. (To be added)
Saving your work¶
Once we generated some data, we would want to save it and give it a name by which it can be recalled later. This is done by using an assignment expression that in k is a colon:
a: &12
From now on, a
will refer to a list of 12 zeros until we reuse this name
by assigning it to something else.
a
0 0 0 0 0 0 0 0 0 0 0 0
Cut and reshape¶
From simple lists, k can create lists of lists by cutting the lists into
chunks using the _
operator. For example, we can take a list of 12 zeros
that we saved in the a
varible before and cut it at positions 0, 2 and 6:
0 2 6 _ a
0 0
0 0 0 0
0 0 0 0 0 0
If we want to cut a list into chunks of equal size, we can use the #
(reshape) operator:
3 4 # a
0 0 0 0
0 0 0 0
0 0 0 0
The reshape operator generalizes even further allowing cutting lists
into lists of lists of lists and so on ( a
is still 12 zeros):
3 2 2 # a
(0 0;0 0)
(0 0;0 0)
(0 0;0 0)
Note that if k had a 3d display, it could show this as a stack of 2x2 matrices, but because we only have two dimensions, k shows stacked matrices in a linear notation. The same notation can be used for input
(1 2;3 4)
1 2
3 4
If the reshape operator is given a list on the right that contains fewer items than is necessary to fill the shape, items from the front of the list will be reused
10#!4
0 1 2 3 0 1 2 3 0 1
This property can be used to build a unit matrix by filling an n x n shape with a length n+1 unit vector:
5 5 # 1,&5
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
0 0 0 0 1
While we have not introduced the ,
operation, you have probably guessed
that in the above example 1,
prepends 1
to the following list. We will
discuss ,
in greater detail next.
Enlist and join¶
K's simple syntax for entering lists has one drawback. It is not obvious how
to enter lists with zero or one element. In fact, k offers no literal syntax for
that and such lists have to be generated. You already know one method: simply
reshape an atom to a list using 0#
or 1#
. When you create a list of one element,
you will see that it is displayed as follows:
1#42
,42
That leading ,
in front of the number is the enlist function that turns the
number 42
into a 1-element list. When applied to a list, ,
turns it into
a 1-element list containing a list. Note the difference between
2#1 2 3 / take first 2 elements
1 2
and
2#,1 2 3 / repeat twice
1 2 3
1 2 3
In the first case, 2#
reshape gets a 3-element list and cuts it to length 2, but
in the second case it gets a 1-element list, so it recycles the first element (which
itself is a list 1 2 3
) and makes a 2 x 3 matrix.
When ,
is placed between two lists or between a list and an atom it joins the elements
together
(1,2 3 4;1 2,3 4;1 2 3,4)
1 2 3 4
1 2 3 4
1 2 3 4
Joining lists of lists or matrices joins the rows:
x,x: e:(!3)=/:!3
1 0 0
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1
(recall that e:(!3)=/:!3
above is a unit 3 x 3 matrix)
To join columns, we can use the ,'
operator:
x,'x: e
1 0 0 1 0 0
0 1 0 0 1 0
0 0 1 0 0 1
To flatten a list of lists, use the ,/
function
,/(1;2 3;4 5 6)
1 2 3 4 5 6
To recursively flatten a nested list, use the ,//
function. Compare
,/(1;(0 1;1 0))
1
0 1
1 0
and
,//(1;(0 1;1 0))
1 0 1 1 0
Lists from atoms and back¶
We can split integers into digits
10\:123456789
1 2 3 4 5 6 7 8 9
and put them back together
10/:1 2 3 4 5 6 7 8 9
123456789
using the pair of vector from scalar (\:
) and scalar from vector (/:
) operators.
For binary expansion - just use 2 instead of 10
2\:42
1 0 1 0 1 0
and we can use the same operator to split time in seconds into days, hours and minutes:
24 60 60\:12345
3 25 45
24 60 60/:3 25 45
12345
Generating 2d data¶
Several primitives exist in k to generate lists of lists in a single operation:
!v
- odometerprm N
- all permutations of numbers from 0 through N-1M cmb N
- all combinations of M numbers drawn from 0 through N-1
Odometer¶
Read this column by column from left to right. The lowest row is counting in base 3. The second row is base 2 and the top row is base 1. The second row changes its values when the lowest row reaches 10 base 3.
!1 2 3
^0 0 0 0 0 0
^0 0 0 1 1 1
0 1 2 0 1 2
Combinatorics¶
prm 3
?0 1 2
?1 0 2
?1 2 0
?0 2 1
?2 0 1
?2 1 0
2 cmb 3
0 1
0 2
1 2
List operations¶
Atomic operations¶
The four arithmetic operator +
, -
, *
and %
can operate
on both scalars and arrays of arbitrary shape. We call these operator
"atomic" because for any scalar function that operates on atoms, this operator
works on each element of an array.
For example, let's
create a triangular shape array that we've seen before and give it
a name b
:
b:0 2 6 _ &12
b
0 0
0 0 0 0
0 0 0 0 0 0
We can add a scalar to b
and it will be added to each element
b + 2
2 2
2 2 2 2
2 2 2 2 2 2
or we can add a vector and its elements will be added to the rows
of b
:
b + 1 2 3
1 1
2 2 2 2
3 3 3 3 3 3
It may take some practice to understand how these rules generalize to deeply nested lists. For example,
(3 2 # !6) + 3 2 2 # &12
(0 0;1 1)
(2 2;3 3)
(4 4;5 5)
In addition to the four arithmetic operations, k applies the same rules to
&
- and/min|
- or/max<
,>
, and=
- comparison
Reduction operations¶
There are ways to reduce lists to single values. Thus first
and count
applied
to a list return the first element and the count of elements respectively
(because of the use of randomness, your results may not be the same as those you see here)
r: 10 rand 300 / generate randomly with replacement
r
243 125 8 155 17 9 4 36 207 230
count r / count of elements
10
first r / first element
243
last r / the last element
230
max r / the largest element (maximum)
243
min r / the smallest element (minimum)
4
avg r / average (arithmetic mean)
103.4
Full list operations¶
In k, we rarely need to process lists one element at a time because we have powerful operations that can transorm the entire list in one go. Thus you can sort a list in either ascending or descending order
asc r
^4 8 9 17 36 125 155 207 230 243
dsc r
243 230 207 155 125 36 17 9 8 4
or you can reverse the list using |
and compare two lists using ~
:
(dsc r) ~ |asc r
1
note that when we compare two lists using ~
, we get a single 1 when they
match and a single 0 when they don't.
Selection operations¶
There are ways to index arrays.
z: 22 + !10
z
22+!10
z[0]
22
z[3]
25
z[3 5]
25 27
z[2+!6]
24 25 26 27 28 29
z[(#z)-1]
31
z[_ (#z) % 2]
27
Now consider multi-dimensional arrays.
mymulti: (1 2 3; 4 5 6; 7 8 9; 10 11 12)
mymulti
1 2 3
4 5 6
7 8 9
10 11 12
/ we interpret mymulti as a four row three column matrix
mymulti[0;0]
1
mymulti[2;0]
7
mymulti[1;2]
6
/ Now we can get full rows
mymulti[1]
4 5 6
/ and full columns
mymulti[;1]
2 5 8 11
Verbs, Adverbs, and User-defined Functions¶
The operations in k are called 'verbs' and often have two meanings depending on whether they are 'unary' (applied to a single argument) or 'binary' (applied to a pair of arguments). Normally, the binary verb will be the more familiar one. Much of the power comes from applying verbs to arrays.
Verbs¶
unary + (flip or transpose)¶
x: (1 2 3 4; 5 6 7 8) / assign to x a two row array whose first row is 1 2 3 4
+x / a four row array (transpose of x) whose first row is 1 5
1 5
2 6
3 7
4 8
binary + (plus)¶
2 + 3 / scalar (single element) addition
5
x + x / array addition
2 4 6 8
10 12 14 16
2 + x / element to array addition
3 4 5 6
7 8 9 10
Just as most human languages have verb modifiers called adverbs,
k does too. They apply to most unary and binary operators.
Thus, the /
adverb (called 'over'), instead of indicating a comment,
can cause the binary version of the verb to apply to the elements
in the array in sequence and yields a single result.
The \
adverb (called 'scan') does the same but keeps all the intermediate
results.
over and scan¶
+/ 1 2 3 4 / Apply the + operator between every pair of elements; produce sum
10
+\ 1 2 3 4 / Same as above but produce all partial sums
1 3 6 10
Adverbs can modify verbs directly but can also modify verb-adverb combinations (which are lifted to verb status). The ' (each) adverb takes both roles.
x
1 2 3 4
5 6 7 8
+/'x / Apply +/ to each row of x
10 26
+\'x / Apply +\ to each row of y
1 3 6 10
5 11 18 26
Adverbs can modify user-defined functions as well.
f:{[a] (a*a)+3 }
f[4]
19
Now we can apply f to each element of an array using the each adverb.
f'1 2 3 4
4 7 12 19
While verbs combined with \ and / have the syntactic form of unary verbs, verbs combined with \: and /: have the syntactic form of binary verbs. Examples:
There is \: (each left):
1 2 3 4 +\: 10
11 12 13 14
There is each right:
20 +/: 1 2 3 4
21 22 23 24
There is each left each right (which should be interpreted as performing an each left on successive elements of the right array):
1 2 3 4 +\:/: 10 20 30 40 50
11 12 13 14
21 22 23 24
31 32 33 34
41 42 43 44
51 52 53 54
Eachright eachleft considers each element of the left array one at a time and applies +/: to the that element and the right array
1 2 3 4 +/:\: 10 20 30 40 50
11 21 31 41 51
12 22 32 42 52
13 23 33 43 53
14 24 34 44 54
This also applies to binary (and even user-defined) verbs.
lambdas¶
g:{[a;b] a + (7 * b)}
g[2;3]
23
Eachleft considers each element of the left array one at a time and applies g to that element and to the entire right array.
1 2 3 4 g\: 10 20
71 141
72 142
73 143
74 144
Eachright considers each element of the right array one at a time and applies g to the left array and that element.
1 2 3 4 g/: 10 20
71 72 73 74
141 142 143 144
Eachleft eachright considers each element of the right array one at a time and applies g\: to the left array and that element.
1 2 3 4 g\:/: 10 20
71 72 73 74
141 142 143 144
(Try for example 1 2 3 4 g\: 20)
Eachright eachleft considers each element of the left array one at a time and applies g/: to that element and the right array
1 2 3 4 g/:\: 10 20
71 141
72 142
73 143
74 144
(Try for example 3 g/: 10 20)
Finally, each can apply to just one argument
x1: 1 2 3 4
x2: 50
g[;x2]'x1
351 352 353 354
g[x1]'x2
351 352 353 354
Extended Example: matrix multiplication¶
Recall that matrix multiplication involves the dot products between rows of the left matrix and the columns of the right matrix.
leftmat: (1 2 3; 4 5 6; 7 8 9; 10 11 12)
leftmat
1 2 3
4 5 6
7 8 9
10 11 12
rightmat: (100 200 300 400 500; 1000 2000 3000 4000 5000; 10000 20000 30000 40000 50000)
rightmat
100 200 300 400 500
1000 2000 3000 4000 5000
10000 20000 30000 40000 50000
dot:{[v1;v2] +/ v1 * v2} / dot product function
dot[4 5 6; 300 3000 30000]
196200
matmult:{[m1;m2] m1 dot/:\: +m2}
matmult[leftmat; rightmat]
32100 64200 96300 128400 160500
65400 130800 196200 261600 327000
98700 197400 296100 394800 493500
132000 264000 396000 528000 660000
Adverbs Replace Loops¶
K programmers tend not to need loops. In fact, some of them disdain loops. The reason is simply that the language uses adverbs instead of loops.
For example, the loop
result = 0
for i = 1 to len(myarray)
result += f(myarray[i])
becomes
result: +/f'array
In principle each invocation of f could (and will eventually) be done in parallel.
By contrast, the loop
result = 0
for i = 1 to len(myarray)
if result = f(result, myarray[i])
becomes / or \ at least for some f's
myarray: 2 2 2 2
f:{[x;y] x + 2*y}
f\ myarray
2 6 10 14
f/ myarray
14
Second example where k's initializations can be useful:
g: {[x;y] x*y}
g\myarray
2 4 8 16
Beyond numbers¶
Names, characters and strings¶
Character strings are simply an array of characters
x: "fast, cool, and really concise"
#x
30
x[2 4]
"s,"
x[<x]
" ,,aaacccdeefilllnnooorssty"
Whereas character strings occupy one byte per character, symbols are hashed and therefore take less space, a useful feature in a large data application in which a symbol is repeated many times.
x: (`abc; `defg)
x
`abc`defg
#x
2
Here is some guidance in choosing between symbols and characters. If there are a few distinct character sequences and they are repeated many times (e.g. a history of all trades where there are only a few thousand stock symbols but millions of trades), then symbols are best for operations like sorting and matching. Otherwise char vectors are probably better especially if you need to do substring matching.
Dates, times and durations¶
Date format is year-month-day and you can get the day by .z.D
x: .z.D
x
2019-12-04
x + 44
2020-01-17
.z.D+/:10 rand 24:00:00 / generate random datetimes
2019-12-04T22:55:55 2019-12-04T03:09:39 2019-12-04T02:03:52 2019-12-04T02:48:02 2019-12-04T18:44:05 2019-12-04T10:35:13 2019-12-04T18:19:04 2019-12-04T07:30:19 2019-12-04T05:40:20 2019-12-04T15:48:08
Time format is hour:minutes:minutes.milliseconds
x: .z.t / greenwich mean time
x
03:51:28.435
x+ 609 / add to milliseconds
03:51:29.044
09:30+10 rand 06:30 / generate random times
10:04 12:55 14:38 10:42 09:53 14:48 11:05 11:49 14:11 12:34
Dictionaries and tables¶
/ Dictionaries are key to value structures. There are many ways to create a dictionary.
/ From partitions
mypart: ="many sentences have the letter e very often"
mypart
|4 14 19 23 30 32 37
a|1 16
c|,11
e|6 9 12 18 22 25 28 31 34 41
f|,39
h|15 21
l|,24
m|,0
n|2 7 10 42
o|,38
r|29 35
s|5 13
t|8 20 26 27 40
v|17 33
y|3 36
mypart["e"]
6 9 12 18 22 25 28 31 34 41
/ Creating them directly
mydict: `bob`carol!(2;3)
mydict
bob :2
carol:3
mydict[`bob]
2
/ dictionaries can be heteogeneous in their values
mydict2: `alice`bill`tom`judy`carol!(1 2 3 4; 8 7; 5; "we the people"; `abc)
mydict2[`carol]
`abc
mydict2[`judy]
"we the people"
/ dictionaries can be heterogeneous in their keys too
mydict3: (`marie;`jeremie;2)!(34; "hello, world"; 7)
mydict3
`marie |34
`jeremie|"hello, world"
2 |7
/ find the keys of dictionaries is easy:
! mydict3
`marie
`jeremie
2
/ and the values
. mydict3
34
"hello, world"
7
Control Flow¶
Before we get to control flow in the classical sense, it's important to remember how to read a k expression. As mentioned earlier, there is no precedence as there is in some languages where for example * binds more than + or -. Instead the precedence is right to left which by the way conforms with mathematical usage (e.g., for sum of x * y we first multiply x and y and then take the sum)
Order of operations¶
20*4-3
20
(20*4) - 3 / to make * bind closer than -, you need parentheses
77
20*(4-3) / otherwise, precedence is right to left.
20
let's say we want the sum of the elements having values greater than 35
x: 90 30 60 40 20 19
x > 35 / put a 1 where values are greater than 35
1 0 1 1 0 0
& x > 35 / indexes that are greater than 35
^?0 2 3
x[&x > 35] / elements in x that are greater than 35
90 60 40
+/ x[&x > 35] / sum of elements of x whose values are greater than 35
190
$[c;t;f]
(Conditional)¶
if c is true then execute the t branch else the f branch)
x: 3 4 5
y: 10 20 30
$[5 > 3; +/x; +/y]
12
$[5 < 3; +/x; +/y]
60
?[x;I;[f;]y] (replace the index positions by what comes afterwards)¶
x: 3 4 5 6 7 8 9 10
y: 100 200 300 400
?[x;3;y] / replace what's in position 3 by y
3 4 5 100 200 300 400 6 7 8 9 10
?[x;3 4; y] / starting in position 3 and counting 4, replace by y
3 4 5 100 200 300 400 10
Input/Output and Interprocess Communication¶
Create a file having these three lines and call it tmp:
We the people
of the
United States
Read in the file
x: 0: "tmp"
x
We the people
of the
United States
x[1] / x is just an array so x[1] is the second element, viz.
"of the"
y: (x[0]; x[1]; x[2]; x[1])
y
We the people
of the
United States
of the
"tmp2" 0: y
Now look at tmp2 and see that you have:
We the people
of the
United States
of the
1:
(write binary image)¶
x: 1 2 3 4
"tmp3" 2: x
y: 2: "tmp3"
y
1 2 3 4
Text input/output is 0:
¶
"foo"2:("This is line 1\n This is line 2")
2:"foo"
"This is line 1\n This is line 2"
Interactive prompt is 1:
¶
name:1:""1:"What is your name? "
What is your name? Carol
name
"Carol"
Data Structures¶
k has atom, list (2;`c), dict [a:2;b:`c] and func {[x;y]x+y} dict [a:2;b:`c] view f::32+1.8*c TODO
Debugging¶
on error(inspect variables and assign to them)
f:{[x;y] x + y}
f[5;6]
11
f[5;`abc]
{[x;y] x + y}
^
type error
> x
5
> y
`abc
> y:7
> x+y
12
'
up a level in the call stack (e.g. if in function f, go to caller of f)
\
out of debugging mode
2+ \3 trace
Scripts and Meta-functions¶
\l a.k
load¶
create a file foo.k with the two lines
x: 1 2 3 4
f: {[x] x*x}
Then start a k session and then load foo.k in that session using the \l command:
\l foo.k
x
1 2 3 4
f
{[x] x*x}
\v
variables \f
functions¶
\l foo.k
\v
,`x
\f
``f
\w
workspace¶
how much memory are you using
\w
1072128
Tables¶
In the row-wise vs. column-wise table debate, k comes out as columnwise. We'll work up to this slowly.
Consider a list:
x: 1 2 3 4 5 10 15
x
1 2 3 4 5 10 15
Create a one column table from this:
xtab: +`numcol!x
xtab
numcol
------
1
2
3
4
5
10
15
Here the table has one column and its header is numcol.
select numcol from xtab
numcol
------
1
2
3
4
5
10
15
select sum numcol from xtab
numcol
------
40
select sum numcol from xtab where numcol > 4
numcol
------
30
Ok, now let's create a multiple column table. Your results may differ because of the use of randomness.
n: 7
newtab: +(`stock`date`price`vol)!(n rand `ibm`goog`hp;.z.D+/:n rand 16:00:00;100 + n rand 200; n rand 5000)
newtab
stock date price vol
----- ------------------- ----- ----
`hp 2019-12-04T10:47:25 127 464
`goog 2019-12-04T00:30:22 116 676
`goog 2019-12-04T02:36:39 163 1711
`ibm 2019-12-04T00:33:30 192 2867
`goog 2019-12-04T09:03:26 145 3420
`goog 2019-12-04T06:08:15 264 1370
`ibm 2019-12-04T03:50:29 195 2498
select sum price*vol by stock from newtab
stock|vol
-----|-------
`goog|1214889
`hp |58928
`ibm |1037574
select sum price*vol by stock from newtab where date > 2018-08-25T10:00:00
stock|vol
-----|-------
`goog|1214889
`hp |58928
`ibm |1037574
User-defined functions:
f:{[x] 1.5*x}
select sum f[price*vol] by stock from newtab where date > 2018-08-25T10:00:00
stock|vol
-----|--------
`goog|1822334f
`hp |88392f
`ibm |1556361f
Extracting Data from Tables into Other Structures
/ select always gives a table
select date from newtab
date
-------------------
2019-12-04T10:47:25
2019-12-04T00:30:22
2019-12-04T02:36:39
2019-12-04T00:33:30
2019-12-04T09:03:26
2019-12-04T06:08:15
2019-12-04T03:50:29
newtab: +`stock`date`price`vol!(n rand `ibm`goog`hp;.z.D+/:n rand 16:00:00;100 + n rand 200; n rand 5000)
newtab
stock date price vol
----- ------------------- ----- ----
`hp 2019-12-04T10:02:30 250 3409
`ibm 2019-12-04T00:32:28 209 4379
`goog 2019-12-04T04:52:07 115 555
`hp 2019-12-04T15:33:23 266 362
`ibm 2019-12-04T00:40:24 267 479
`hp 2019-12-04T11:51:26 158 3163
`ibm 2019-12-04T14:50:35 179 895
Importing from a csv file
Create a small csv file mytrade.csv whose schema is: tradeid,stock,timeindicator,price,vol
mytrade.csv:
1,goog,50,1237,100
2,msft,51,109,100
3,goog,52,1240,200
4,msft,53,112,200
("iniii";",")0:"mytrade.csv"
1 2 3 4
`goog`msft`goog`msft
50 51 52 53
1237 109 1240 112
100 100 200 200
mytrade1: +(`tradeid`stock`timeindicator`price`vol)!("iniii";",")0:"mytrade.csv"
select tradeid, price, vol from mytrade1
tradeid price vol
------- ----- ---
1 1237 100
2 109 100
3 1240 200
4 112 200
select tradeid, price, vol from mytrade1 where price > 500
tradeid price vol
------- ----- ---
1 1237 100
3 1240 200
Then one with proper datetimestamps called mytradebac2.csv
1,goog,15:16:50,1237,100
2,msft,15:16:51,109,100
3,goog,15:18:50,1240,200
4,msft,15:18:52,112,200
mytrade2: +(`tradeid`stock`time`price`vol)!("intfi";",")0:"mytradebac2.csv"
select tradeid, time, price, vol from mytrade2 where price > 500
tradeid time price vol
------- ------------ ----- ---
1 15:16:50.000 1237f 100
3 15:18:50.000 1240f 200
/ Note that it is also possible to import from just a set of lists, / because that is what 0:"somefile.csv" gives.
0:"mytradebac2.csv"
1,goog,15:16:50,1237,100
2,msft,15:16:51,109,100
3,goog,15:18:50,1240,200
4,msft,15:18:52,112,200
select sum price * vol by `minute$time from mytrade2
time |vol
-----|-------
15:16|134600f
15:18|270400f
select `minute$time from mytrade2
time
-----
15:16
15:16
15:18
15:18
Modifying Tables
select vol from mytrade2
vol
---
100
100
200
200
update vol:1+vol from mytrade2
tradeid stock time price vol
------- ----- ------------ ----- ---
1 goog 15:16:50.000 1237 101
2 msft 15:16:51.000 109 101
3 goog 15:18:50.000 1240 201
4 msft 15:18:52.000 112 201
/ but the table itself is not changed because there no assignment:
mytrade2
tradeid stock time price vol
------- ----- ------------ ----- ---
1 `goog 15:16:50.000 1237f 100
2 `msft 15:16:51.000 109f 100
3 `goog 15:18:50.000 1240f 200
4 `msft 15:18:52.000 112f 200
/ On the other hand
mytrade2updated: update vol:1+vol from mytrade2
mytrade2updated
tradeid stock time price vol
------- ----- ------------ ----- ---
1 goog 15:16:50.000 1237 101
2 msft 15:16:51.000 109 101
3 goog 15:18:50.000 1240 201
4 msft 15:18:52.000 112 201
delete from mytrade2 where vol > 100
tradeid stock time price vol
------- ----- ------------ ----- ---
1 goog 15:16:50.000 1237 100
2 msft 15:16:51.000 109 100
/ Would need to assign to mytrade2 to see this effect. / e.g. mytrade2: delete from mytrade2 where vol > 100
/ Here is a row to insert.
x: {tradeid:15;stock:`goog;time:15:26:50.123;price:2337f;vol:200}
/ An insert:
mytrade2: mytrade2,x
/ There is a notion of a keyed table where each key value is supposed / to occur only once. tradeid is an example.
u:`tradeid key mytrade2
u
tradeid|stock time price vol
-------|----- ------------ ----- ---
1 |goog 15:16:50.000 1237 100
2 |msft 15:16:51.000 109 100
3 |goog 15:18:50.000 1240 200
4 |msft 15:18:52.000 112 200
15 |goog 15:26:50.123 2337 200
/ Note that the , operator will insert if the new row has a new key / (tradeid of 45)
y: {tradeid:45;stock:`goog;time:16:26:50.123;price:3337f;vol:75}
u,y
tradeid|stock time price vol
-------|----- ------------ ----- ---
1 |goog 15:16:50.000 1237 100
2 |msft 15:16:51.000 109 100
3 |goog 15:18:50.000 1240 200
4 |msft 15:18:52.000 112 200
15 |goog 15:26:50.123 2337 200
45 |goog 16:26:50.123 3337 75
/ but update if the new row has an existing key (tradeid of 4). / So the , operator is called an upsert.
y: {tradeid:4;stock:`goog;time:16:26:50.123;price:3337f;vol:75}
u,y
tradeid|stock time price vol
-------|----- ------------ ----- ---
1 |goog 15:16:50.000 1237 100
2 |msft 15:16:51.000 109 100
3 |goog 15:18:50.000 1240 200
4 |goog 16:26:50.123 3337 75
15 |goog 15:26:50.123 2337 200
/ This is an upsert because this is an operation that specifies / the key and all fields.
u
tradeid|stock time price vol
-------|----- ------------ ----- ---
1 |goog 15:16:50.000 1237 100
2 |msft 15:16:51.000 109 100
3 |goog 15:18:50.000 1240 200
4 |msft 15:18:52.000 112 200
15 |goog 15:26:50.123 2337 200
Shortcuts¶
1) We would be remiss to fail to mention some shortcuts that k afficionados love to use, even though some of us feel that they reduce clarity. For example, unary functions implicitly perform "each" when applied to arrays. For example,
f:{[a] (a*a)+3 }
f'1 2 3 4
4 7 12 19
f 1 2 3 4
4 7 12 19
x:(1 2 3 4;5 6 7 8)
f'x
4 7 12 19
28 39 52 67
f x
4 7 12 19
28 39 52 67
default function parameters are x y z, e.g. {z+x*y}[3;2;1] is 7
f:{z + x*y}
f[10; 20; 30]
230
Eval¶
It is possible to make strings be evaluated as k expressions using the dot operator:
. "2+3"
5
System Calls (in progress)¶
e.g. instead of \ls \du \wc -l try
Unix ls is just \ff ?.c
Unix du is \fk ?.c
Unix wc -l is \fl ?.c
A gallery of exercises/examples¶
These examples are going to start from a database of trades.
n: 100
secid: n rand (`goog;`facebook;`ibm;`msft)
price: 100 + n rand 200
vol: 10 + n rand 1000
time: !n
Find all trades such that the price is over 175.
ii: & price > 175
z: secid[ii] ,' price[ii] ,' vol[ii] ,' time[ii]
z[5]
`facebook
289
75
10
Find the high and low of each security.
mydict: =secid
names: rand secid
maxmin: {[name] name , (|/price[mydict[name]]), (&/price[mydict[name]])}
maxmin'names
`goog
282
104
/ or to be able to take an arbitrary dictionary
maxmin2:{[name;somedict] name , (|/price[somedict[name]]), (&/price[somedict[name]])}
maxmin2[;mydict]'names
`goog
282
104
Get the moving average of each security
mavg:{[name] name, (+\price[mydict[name]])%'(1+!#mydict[name])}
mavg'names
`goog
164f
179.5
154.3333
171.25
166.6
178.1667
187.4286
185.5
196.2222
195.8
202.6364
198.3333
200.6154
197.9286
202.6
200f
198f
200.6111
200.9474
203.05
206f
Determining the finishing time of each task if you do them in earliest deadline first order.
n: 10
taskid: !n
tasktime: 2 + n rand 20
deadlines: 40 + n rand 50
tasktime
21 4 6 9 14 10 16 18 17 2
deadlines
65 46 75 71 79 70 74 87 62 81
/ Put them all together
taskid,'tasktime,'deadlines
0 21 65
1 4 46
2 6 75
3 9 71
4 14 79
5 10 70
6 16 74
7 18 87
8 17 62
9 2 81
/ Now determine the order of deadline indexes / for the deadlines to be in order.
inddead: < deadlines
deadlines[inddead]
46 62 65 70 71 74 75 79 81 87
/ Put the tasks in the same order
taskid[inddead]
1 8 0 5 3 6 2 4 9 7
/ Put the task times in the same order
tasktime[inddead]
4 17 21 10 9 16 6 14 2 18
/ Put all of them in the same order
taskid[inddead],'tasktime[inddead],'deadlines[inddead]
1 4 46
8 17 62
0 21 65
5 10 70
3 9 71
6 16 74
2 6 75
4 14 79
9 2 81
7 18 87
/ Find end point if tasks are executed in this order
+\tasktime[inddead]
4 21 42 52 61 77 83 97 99 117
/ Determine which deadlines are met (1) and which aren't (0)
deadlines[inddead] > +\tasktime[inddead]
1 1 1 1 1 0 0 0 0 0
/ Determine which taskids have their deadlines met
taskid[inddead][& deadlines[inddead] > +\tasktime[inddead]]
1 8 0 5 3
Debugging¶
fib:{[n] fib[n-1] + fib[n-2]}
fib[5]
{[n] fib[n-1] + fib[n-2]}
^
stack error
To debug this, we can do several thing. First, just query the variables, e.g.
> n
-189
We might realize that n should never be negative. Another thing we can do is store all the values of n in this recursive function.
out: ()
fib:{[n] out,: n; fib[n-1] + fib[n-2]}
fib[5]
{[n] out,: n; fib[n-1] + fib[n-2]}
^
stack error
>
But now we can query out:
> out
5 3 1 -1 -3 -5 -7 -9 -11 -13 -15 -17 -19 -21 -23 -25 -27 -29 -31 -33 -35 -37 -39 -41 -43 -45 -47 -49 -51 -53 -55 -57 -59 -61 -63 -65 -67 -69 -71 -73 -75 -77 -79 -81 -83 -85 -87 -89 -91 -93 -95 -97 -99 -101 -103 -105 -107 -109 -111 -113 -115 -117 -119 -121 -123 -125 -127 -129 -131 -133 -135 -137 -139 -141 -143 -145 -147 -149 -151 -153 -155 -157 -159 -161 -163 -165 -167 -169 -171 -173 -175 -177 -179 -181 -183 -185
Order-based Relational Algebra¶
This version is build just on arrables (tables consisting of lists of ordered lists (arrays)). First we review selects, projects, and moving aggregates. Then we show equi-joins then general joins. file: arrable.k
/ functions
/ given a set of indexes give me those values of a vector x
i:1 5 7
x:0 10 20 30 40 50 60 70 80
x @ i
10 50 70
/ given a bunch of lists alllist / a selection string on alllist selstr / (a typical selstr might be "(alllist[1] > 3 ) & (alllist[0] < 40)" ) / and a set of output columns outlist / Output: after selecting based on selstr, output columns cols of alllist
mysel:{[alllist;selstr;outlist] ii: & . selstr;alllist[outlist]@\:ii}
/ given a bunch of lists alllist / a subset to sort by sortby / other lists to follow that sort outlist / Output: based on the sort order of alllist[sortby] / create rows of outlist in order
myasc:{[alllist;sortby;outlist] myind: $[1 < #sortby; < +alllist[sortby]; < alllist[*sortby]]; alllist[outlist]@\:myind}
/ This does moving sum on numpoints (e.g. three point moving sum) / of the array myarray / If there are fewer than numpoint in myarray, it does the moving sum up / to the number of points in myarray.
movsum:{[numpoints;myarray] x: +\myarray; xsub: $[numpoints < #myarray;(numpoints # 0),x[!(#x)-numpoints];0]; x - xsub}
/ This does moving average on numpoints (e.g. three point moving average) / of the array myarray / If there are fewer than numpoint in myarray, it does the moving average up / to the number of points in myarray.
movavg:{[numpoints;myarray] x: movsum[numpoints; myarray]; mydivs: $[numpoints < #myarray; (1+!numpoints), numpoints _ (#myarray) # numpoints;(1+!#myarray)]; x%mydivs}
/ relational equijoin on one attribute
/ given two lists of lists LL1 and LL2
/ index from LL1 indLL1
/ index from LL2 indLL2
/ indexes from LL1 outLL1
/ indexes from LL2 outLL2
/ Find all indexes of LL1 and indexes of LL2 that match
/ based on the values in indLL1 and indLL2 and then take the cross product
/ for the columns outLL1 of LL1 and outLL2 of LL2
eqjoin:{[LL1;LL2;indLL1;indLL2;outLL1;outLL2] mymatch: &:' LL1[indLL1] =\: LL2[indLL2]; outindLL1: ,/ ((#:')mymatch) #' !#LL1[indLL1]; outindLL2: ,/ mymatch; (LL1[outLL1]@\:outindLL1),LL2[outLL2]@\:outindLL2}
/ Below is unused (and needs to be debugged) until we can get longer functions
/ mycross:{[pair; mydict1; mydict2] x: pair[0]; y: pair[1]; ,/mydict1[x] ,/:\: mydict2[y]}
/ fin:{[allmatches; Ll1; LL2; outLL1; outLL2] (indtoval[;allmatches[;0]]'LL1[outLL1]), (indtoval[;allmatches[;1]]'LL2[outLL2])}
/ eqjoindict:{[LL1;LL2;indLL1;indLL2;outLL1;outLL2] d1: = LL1[indLL1]; keys1: ! d1; d2: = LL2[indLL2]; keys2: ! d2; mym: keys1 rand keys2; pairs: (mym,'keys2); pairs@: & pairs[;0] < #keys1; mm: ,/mycross[;d1;d2]'(keys1[pairs[;0]],' pairs[;1]); fin[mm; LL1; LL2; outLL1; outLL2]}
/ data
indata: (3 4 3 4 3 9 9 9 9 9 9;30 40 30 40 30 90 90 90 90 90 90;7 9 1 2 1 2 3 4 8 7 3)
/ execution
mysel[indata; "indata[1] > 60 "; 0 1 2]
9 9 9 9 9 9
90 90 90 90 90 90
2 3 4 8 7 3
x: myasc[indata; 2 1; ,0]
x
,3 3 4 9 9 9 9 3 9 9 4
movavg[3;x]
,3 3 4 9 9 9 9 3 9 9 4f
/ can combine these
x1: mysel[indata; "indata[1] > 35 "; 0 1 2]
x2: myasc[x1; 2 1; ,0]
movavg[3;x2]
,4 9 9 9 9 9 9 4f
/ Now look at the equijoin
LL1: indata
LL2: (300 400 300 400 300 800 800 800 301 401 402;30 40 30 40 30 80 80 80 30 40 40;7 9 1 2 1 2 3 4 8 7 3)
eqjoin[LL1; LL2; 1; 1; 0 1; 0 1]
3 3 3 3 4 4 4 4 3 3 3 3 4 4 4 4 3 3 3 3
30 30 30 30 40 40 40 40 30 30 30 30 40 40 40 40 30 30 30 30
300 300 300 301 400 400 401 402 300 300 300 301 400 400 401 402 300 300 300 301
30 30 30 30 40 40 40 40 30 30 30 30 40 40 40 40 30 30 30 30
/ Sometimes we want to perform a join, but get results only for one argument / e.g. LL1 in this case:
y1: eqjoin[LL1; LL2; 1; 1; 0 1 2; !0]
y2: mysel[y1; "y1[1] > 35 "; 0 1 2]
y3: myasc[y2; 1; ,2]
movavg[3;y3]
,9 9 9 9 2 2 2 2f
String matching (dynamic programming example)
/ data
s1: "rent"
s2: "let"
initval: 1 + ((#s1) * (#s2))
mat: (1+(#s1);1+(#s2)) # initval
mat[0]: 0, 1+!#s2
mat: mat@[;0;:;]'0, 1+!#s1
k reference¶
Until now, we have touched on only a few of the verbs and types. Here is Arthur Whitney's full list. From what you understand already, these won't be hard to learn.
+ plus flip
- minus negate
* times first
% divide inverse
& min|and where
| max|or reverse
< less up
> more down
= equal group
~ match not
! dict|mod key|enum
, concat enlist
^ except null
$ pad|cast string
# take|select count
_ drop|delete floor
? find|rand uniq|rand
@ index type
. apply value
First, let's look at how to read this table. In each row, the binary meaning precedes the unary meaning.
Let's go through each in turn.
/ this is a comment \\ if alone on a line exits the k session or a debugging environment
:
gets¶
x: 1 2 3 4 / gets indicates assignment
+
plus flip¶
2 + 3
5
/ Unary + (transpose)
+ (1 2 3 4; 5 6 7 8)
1 5
2 6
3 7
4 8
-
minus negate¶
2 - 3
-1
/ Unary - (negation)
- 3 4 5
-3 -4 -5
*
times first¶
4*5
20
/ Unary * (first in list)
* 15 24 19 10
15
%
divide inverse¶
5 % 3
1.666667
/ Unary % (inverse)
% 81
0.01234568
!
mod|div enum¶
/ The following is 14 mod 3 (it turns out to be convenient to put the modulus first)
3\ 14
2
/ Unary ! enumerate either by integer (in this case numbers 0, 1, 2, ... 19) or float
!20
!20
&
min|and (and, if 1 is interpreted as true and 0 as false) where¶
5 & 3
3
1 & 1
1
1 & 0
0
/ Unary & (indexes where a there is a non-zero)
x: 4 8 9 2 9 8 4
x = 9 / 1 will indicate a match
0 0 1 0 1 0 0
& x = 9 / locations in the list above that are 1
^?2 4
x > 4
0 1 1 0 1 1 0
& x > 4
^?1 2 4 5
|
max|or reverse¶
5 | 3
5
1 | 0
1
/ Unary | reverses lists
| 2 3 4 5 6
6 5 4 3 2
<
less asc¶
5 < 3 / returns 0 because false
0
3 < 5 / returns 1 because true
1
/ Unary < says which order of indices gives data in ascending order
x: 6 2 4 1 10 4
xind: < x / index locations from smallest value to highest
xind / Notice that x[3] is 1, the lowest value in the list
?3 1 2 5 0 4
x[xind] / sort the values in ascending order
1 2 4 4 6 10
>
more dsc¶
5 >gt; 3
1
3 > 5
0
x: 6 2 4 1 10 4
/ Unary > says which order of indices gives data in descending order
xind: > x / index locations from highest value to smallest
xind
?4 0 2 5 1 3
x[xind] / descending sorted order
10 6 4 4 2 1
=
equal group¶
5 = 5
1
1 2 13 10 = 1 2 13 10 / element by element
1 1 1 1
1 2 14 10 = 1 2 13 10 / 0 at position 2 indicates inequality
1 1 0 1
/ Unary = gives a dictionary (more on dictionaries below) mapping values to indexes where / values are present
= 30 20 50 60 30 50 20 20
20|1 6 7
30|0 4
50|2 5
60|,3
/ In the above example 20 is at positions 1 6 and 7
x: = 30 20 50 60 30 50 20 20
~
match not¶
1 2 13 10 ~ 1 2 13 10 / should match
1
1 2 10 13 ~ 1 2 13 10 / does not match (order matters)
0
/ Unary ~ (like the Boolean not operator except that any non-zero becomes zero)
~ 1
0
~ 0
1
~ 5
0
~ -5
0
,
concat enlist¶
x: 10 11 12 13
y: 4 3
x,y / concatenate one way
10 11 12 13 4 3
y,x / concatenate the other
4 3 10 11 12 13
z: y,x
z[4] / can index these as if these were an array
12
z[2+!3] / can fetch many indexes
10 11 12
/ Unary , converts a scalar (atom) into a list or a list into a deeper list
x: 5
x / x is an atom
5
*x / first on an atom is just the atom itself
5
x: ,5 / x is now a list
x / Note the comma in front indicating that x is list
,5
y: x,15
y / the comma goes away since two or more elements already form a list
5 15
*y / The * operator takes the first element of a list
5
*x / The * operator takes the first element of a list even a singleton
5
^
except null¶
x: 8 7 6 5 2
x ^ y / elements of x (preseving order in x) that are not in y
8 7 6 2
y ^ x
,15
#
take|shape count¶
x: 10 20 30 40 50 60
3 # x
10 20 30
10 # x / Notice x is of length 6; result wraps
10 20 30 40 50 60 10 20 30 40
3 2 # x / creates a three row, two column matrix
10 20
30 40
50 60
3 10 # x / creates a three row, 10 column matrix (with wrapping)
10 20 30 40 50 60 10 20 30 40
50 60 10 20 30 40 50 60 10 20
30 40 50 60 10 20 30 40 50 60
/ Unary # counts the lengths of lists
x
10 20 30 40 50 60
# x / number of elements in x
6
y: 9 8 5
z: (x;y)
z
10 20 30 40 50 60
9 8 5
#:'z / counts each list
6 3
_
drop|cut floor¶
x: 10 20 30 40 50 60
3 _ x / cut away 3 elements from the beginning
40 50 60
10 _ x / Notice x is of length 6; so this eliminates more than necessary
!0
/ !0 means an empty list
/ Now unary _ is the floor operator
15 % 4
3.75
_ 15 % 4
3
$
cast|+/* string¶
` $ "abc" / cast string to symbol (name)
`abc
. "18" / cast string to int
18
. "18.2" / cast string to float
18.2
/ Unary form
$ `abc / cast symbol to string
"abc"
rand
rand|find unique¶
/ Random as discussed in section 1
10 rand 12 / with replacement (can be duplicates) from 0 to 11
10 1 9 2 4 4 7 10 10 10
15 rand 12 / there can be more elements (15) than the domain (0 to 11)
9 6 3 2 10 1 0 6 10 5 8 1 10 1 3
-10 rand 12 / random and uniform without replacemnt (no duplicsates)
1 7 2 5 0 3 4 8 10 6
/ With list as left argument we can find the index of the first match
40 20 30 10 20 30 ? 30
2
/ Unary rand
for atoms (scalars)
rand 7 / random between 0 and 1 (uniform)
0.9327079 0.2525184 0.9101379 0.6989249 0.6978369 0.238826 0.8294105
/ Unary ? for arrays removes duplicates but preserves order
? 40 20 30 10 20 30
?40 20 30 10
@
at type¶
x: 40 20 30 10 20 30
@[x;2 4 5]
30 20 30
x / unchanged
40 20 30 10 20 30
@[x; 2 4 5;: ; -17 -12 -8]
40 20 -17 10 -12 -8
x / still unchanged
40 20 30 10 20 30
f: {[x] x * x}
@[x; 2 4 5; f] / squares locations 2 4 and 5
40 20 900 10 400 900
x / still unmodified
40 20 30 10 20 30
@[x; 2 4 5; f] / squares locations 2 4 and 5
40 20 900 10 400 900
/ Unary @ finds the type of an object
@ 18
`i
@ "18"
`c
@ `abc
`n
.
dot value¶
/ Unary . Can evaluate a string
. "18 + 5"
23
. "f: {[x] x * x * x}"
. "f[5]"
125
21) abs (absolute value)
abs -3.2
3.2
abs -3.2 4 5.3
3.2 4 5.3
log
(natural log, also known as ln or log base e)¶
log 8
2.079442
exp
(exponential on e)¶
exp 1
2.718282
exp 3
20.08554
log 20.08554
3f
sin
(takes its argument in radians)¶
mypi: 3.14 / a crude approximation of pi
sin[mypi] / should be approximately 0
0.001592653
sin[mypi % 2] / sin pi/2 is 1, so this is close
0.9999997
cos
also takes its argument in radians¶
cos[mypi] / close to -1
-0.9999987
cos[mypi %2]
0.0007963267
in
(membership test)¶
5 in 10 20 30 5 6 9 10
1
5 in 10 20 30 5 6 5 10 / even if there are two instances, still return 1
1
5 in 10 20 30 50 6 15 10
0
bin
(binary search assumes ascending sorted order))¶
x: 5 * (3 + !10)
x
15 20 25 30 35 40 45 50 55 60
x bin 33 / index location that is less than or equal to 33
3
x bin 35 / index location that is less than or equal to 35 (here, equal)
4
within
¶
(upper bound for singleton right hand side lists and closed lower and open upper bound for binary right hand side lists)
x: 40 10 20 23 15 16 18
x within 0 20 / less than 20
0 1 0 0 1 1 1
x
40 10 20 23 15 16 18
x within 16 23 / between 16 (inclusive) and 23 (exclusive)
0 0 1 0 0 1 1
/ No unary version
find
/ substring looking for an exact match¶
x: "abcdef"
x find "cde" / look for beginning and length of match
,2 3
x find "bd" / If there is no match, then is this the return value I want???
0#,Ø Ø
y: x, x
y
"abcdefabcdef"
y find "cde" / get all matches as a list
2 3
8 3
like
(string match with wildcards)¶
x: "abcdef"
x like "ab"
0
x like "ab*" / allows wildcards
1
x like "*def"
1
x like "*bcd*" / arbitrary length wildcards with *
1
x like "?bcd*" / ? is a single character substitution
1