# Circular programming

### From HaskellWiki

EndreyMark (Talk | contribs) ((1) Link to HaWiki's CircularProgramming (2) Initiating Attribute grammar page) |
EndreyMark (Talk | contribs) (→Porting from HaWiki's CircularProgramming: but only my writings are ported, because I do not know the copyright status of the other contributions) |
||

Line 1: | Line 1: | ||

− | HaWiki's | + | == Porting from HaWiki's CircularProgramming == |

− | + | I port only my material from there (because I know only the copyright status of my writings on HaWiki). | |

+ | So, the previous materials (to which mine is connected) are on | ||

+ | HaWiki's [http://www.haskell.org/hawiki/CircularProgramming CircularProgramming] | ||

+ | |||

+ | <haskell> | ||

+ | -- | ||

+ | -- | ||

+ | -- THE LOOP OF THE ARROW --- BEING ITSELF ALSO AN ARROW | ||

+ | -- | ||

+ | -- ,..............................................., | ||

+ | -- : : | ||

+ | -- : : | ||

+ | -- : +---------+ : | ||

+ | -- : INPUT |:::::::::| OUTPUT : | ||

+ | -- ==========================>+:::::::::+========================> | ||

+ | -- : |::ARROW::| : | ||

+ | -- : +====>+:::::::::+=====+ : | ||

+ | -- : | |:::::::::| | : | ||

+ | -- : | +---------+ | : | ||

+ | -- : | | : | ||

+ | -- : +==========<==========+ : | ||

+ | -- : FEEDBACK : | ||

+ | -- : : | ||

+ | -- ;...............................................; | ||

+ | -- | ||

+ | </haskell> | ||

+ | |||

+ | Ross Paterson: [http://www.soi.city.ac.uk/~ross/papers/fop.html Arrows and Computation] introduces the <hask>loop</hask> construct of arrows (on page 11) | ||

+ | <haskell> | ||

+ | class Arrow arrow => ArrowLoop arrow where | ||

+ | loop :: arrow (input, feedback) (output, feedback) -> arrow input output | ||

+ | </haskell> | ||

+ | with an ordinary function called <hask>trace</hask>: | ||

+ | <haskell> | ||

+ | trace :: ((input, feedback) -> (output, feedback)) -> input -> output | ||

+ | trace f input = let (output, feedback) = f (input, feedback) | ||

+ | in output | ||

+ | </haskell> | ||

+ | (I renamed the type variables and their corresponding argument names, too.) | ||

+ | |||

+ | Of course <hask>trace</hask> can be seen as <hask>loop</hask> for ordinary functions as arrows: | ||

+ | <haskell> | ||

+ | instance ArrowLoop (->) where | ||

+ | loop = trace | ||

+ | </haskell> | ||

+ | |||

+ | Let us use a modified version of <hask>trace</hask> here: | ||

+ | our version of <hask>trace</hask> should require not a function from an ordered pair (to an ordered pair), | ||

+ | but a curried function (to an ordered pair): | ||

+ | |||

+ | <haskell> | ||

+ | trace :: (input -> feedback -> (output, feedback)) -> input -> output | ||

+ | trace f input = let (output, feedback) = f input feedback | ||

+ | in output | ||

+ | </haskell> | ||

+ | |||

+ | What can <hask>trace</hask> be good for? | ||

+ | For example, for modularizing that magic feedback step in circular programs! | ||

+ | By the concept of modularity I mean as described in [http://www.cs.chalmers.se/~rjmh/ John Hughes]' article [http://www.math.chalmers.se/~rjmh/Papers/whyfp.html Why Functional Programming Matters]. | ||

+ | |||

+ | Let us see the tree minimum problem (described above), and let us modularize the magic circular step by this feedback view of <hask>trace</hask>: | ||

+ | |||

+ | <haskell> | ||

+ | repmin :: Ord a => Tree a -> Tree a | ||

+ | repmin = trace repIImin | ||

+ | |||

+ | repIImin :: Ord a => Tree a -> a -> (Tree a, a) | ||

+ | repIImin (Leaf minval) rep = (Leaf rep, minval) | ||

+ | repIImin (Branch left right) rep = let (left', min_left) = repIImin left rep | ||

+ | (right', min_right) = repIImin right rep | ||

+ | in (Branch left' right', min min_left min_right) | ||

+ | </haskell> | ||

+ | |||

+ | Notations: in | ||

+ | <math>\mathrm{rep}\Vert\mathrm{min}</math>, | ||

+ | my -<code>II</code>- infix notation means parallelity, | ||

+ | referring to the fact that the wonderful achievement of these functions is being able to work in one traversal. | ||

+ | Sometimes I use -<code>X</code>- infix notation instead of -<code>II</code>-, referring to cross product: | ||

+ | <math>\mathrm{rep}\times\!\mathrm{min}</math> | ||

+ | |||

+ | As we can see, we have things enough to confuse: two arguments, a returned tuple, how they depend from each other... | ||

+ | A hint: let us think of the scheme of <hask>trace</hask>: | ||

+ | * input -- what we work with | ||

+ | * output -- what we get | ||

+ | * feedback -- a state-like thing, yes, a feedback, that can be thought of through both aspects | ||

+ | ** feedback-as-produced: thinking of it as it leaves the box | ||

+ | ** feedback-as-fed-back: thinking of it as it is fed back into the box again | ||

+ | If we decide, what the roles of input, output, feedback plays, then these circular programming funtions will look easier. E.g. here, at the repmin problem | ||

+ | * input is a tree filled with numbers, a heterogenous tree | ||

+ | * output is the repmin'ed tree, a homogenized tree | ||

+ | * feedback (something state-like): here that strange numeric value plays the role of feedback! Its two aspects are | ||

+ | ** feedback-as-produced -- in this aspect, it is a minimum value, a ''statistic'' of the tree | ||

+ | ** feedback-as-fed-back -- in this aspect, it is a setter value, a rather ''seed''-like thing: it is spred homogenously in the shape of a given tree | ||

+ | Now it is easier to look at | ||

+ | <math> | ||

+ | \mathrm{rep}\Vert\mathrm{min} | ||

+ | </math>, we have a road-map for our eyes in the jungle of parameters and returned tuple. | ||

+ | |||

+ | <haskell> | ||

+ | -- | ||

+ | -- | ||

+ | -- How repmin is born: wrapping rep||min with trace | ||

+ | -- | ||

+ | -- | ||

+ | -- ,.........................repmin............................, | ||

+ | -- : : | ||

+ | -- : : | ||

+ | -- : +------------+ : | ||

+ | -- : Heterogenous tree |::::::::::::| Homogenized tree : | ||

+ | -- =============================>+::::::::::::+================================> | ||

+ | -- : |::rep||min::| : | ||

+ | -- : +====>+::::::::::::+=====+ : | ||

+ | -- : | |::::::::::::| | : | ||

+ | -- : setter-val | +------------+ | min-val : | ||

+ | -- : a seed | | a statistic : | ||

+ | -- : +===========<============+ : | ||

+ | -- : FEEDBACK : | ||

+ | -- : : | ||

+ | -- ;...........................................................; | ||

+ | -- | ||

+ | </haskell> | ||

+ | |||

+ | |||

+ | Normalizing vectors by <hask>trace</hask> is not so straightforward, but a small trick can help us: | ||

+ | |||

+ | <haskell> | ||

+ | normalize :: Floating a => [a] -> [a] | ||

+ | normalize = trace divideIInorm | ||

+ | |||

+ | divideIInorm :: Floating a => [a] -> a -> ([a], a) | ||

+ | divideIInorm vector a = let (divided, norm2) = scaleIInorm2 vector (recip a) | ||

+ | in (divided, sqrt norm2) | ||

+ | |||

+ | scaleIInorm2 :: Floating a => [a] -> a -> ([a], a) | ||

+ | scaleIInorm2 [] _ = ([], 0) | ||

+ | scaleIInorm2 (x : xs) a = let (scaled, norm2) = scaleIInorm2 xs a | ||

+ | in (a * x : scaled, x * x + norm2) | ||

+ | </haskell> | ||

+ | |||

+ | (I mean my notation <hask>scaleIInorm2</hask> be read as | ||

+ | <math>\mathrm{scale}\Vert\mathrm{norm}^2</math>). | ||

+ | We can see what the problem was and how it was solved here: | ||

+ | * <hask>trace</hask> requires a rather strict scheme: the function it takes as an argument must be exactly of the same scheme as seen in the above feedback picture | ||

+ | * the recursion by which we can solve the problem eventually (here: math>\mathrm{scale}\Vert\mathrm{norm}^2</math>) usually provides another scheme. E.g. here we cannot compute norm directly by recursion. | ||

+ | So we have to join together the needs of <hask>trace</hask> and the possibilities of recursion together, | ||

+ | we have to convert between these two plugs if they are incompatible, | ||

+ | that is why we created | ||

+ | <math>\mathrm{divide}\Vert\mathrm{norm}</math>. | ||

+ | |||

+ | This is not the only way to force the problem to the shape of <hask>trace</hask>. | ||

+ | There are other ways, but they do not seem so natural, | ||

+ | they look somewhat artificial, | ||

+ | albeit in efficiency aspects, I think, they are the same: | ||

+ | |||

+ | <haskell> | ||

+ | normalize :: Floating a => [a] -> [a] | ||

+ | normalize = trace divideBySqrtIInorm2 | ||

+ | |||

+ | divideBySqrtIInorm2 :: Floating a => [a] -> a -> ([a], a) | ||

+ | divideBySqrtIInorm2 vector scalar = scaleIInorm2 vector (1 / sqrt scalar) | ||

+ | </haskell> | ||

+ | |||

+ | works too, and there is another artificial solution yet: | ||

+ | |||

+ | <haskell> | ||

+ | normalize :: Floating a => [a] -> a -> ([a], a) | ||

+ | normalize = trace scaleIIrecipNorm | ||

+ | |||

+ | scaleIIrecipNorm :: [a] -> a -> ([a], a) | ||

+ | scaleIIrecipNorm vector a = let (scaled, norm2) = scaleIInorm2 vector a | ||

+ | in (scaled, recip $ sqrt norm2) | ||

+ | </haskell> | ||

+ | |||

+ | this is good, too. All these three solutions converge, they terminate. | ||

+ | |||

+ | But there are circular programming examples | ||

+ | where such possible solutions will not be so equivalent: | ||

+ | one of them may diverge, while the other converge. | ||

+ | It seems to me, the diff example (described in Wouter Swiestra's [http://www.haskell.org/tmrwiki/WhyAttributeGrammarsMatter Why Attribute Grammars Matter]) is exactly such circular programming problem. | ||

+ | |||

+ | This diff example is: | ||

+ | * how to compute the average of a list | ||

+ | * and (in the same traversal of the list) decrease all elements of the list by this average value. | ||

+ | |||

+ | Here, the gap between | ||

+ | * the needs of `trace` | ||

+ | * and possibilities of (single traversal) recursion | ||

+ | seems even wider than in the previous problem. Let us try: | ||

+ | |||

+ | <haskell> | ||

+ | diff :: Fractional a => [a] -> [a] | ||

+ | diff = trace decrementIIaverage | ||

+ | |||

+ | decrementIIaverage :: Fractional a => [a] -> a -> ([a], a) | ||

+ | decrementIIaverage list delta = let (decremented, (sum', length')) = decrementIIsumLength list delta | ||

+ | in (decremented, sum' `divideByInt` length') | ||

+ | |||

+ | decrementIIsumLength :: Fractional a => [a] -> a -> ([a], (a, Integer)) | ||

+ | decrementIIsumLength [] _ = ([], (0, 0)) | ||

+ | decrementIIsumLength (a : as) delta = let (decremented, (sum', length')) = decrementIIsumLength as delta | ||

+ | in (a - delta : decremented, (a + sum', succ length')) | ||

+ | |||

+ | divideByInt :: Fractional a => a -> Integer -> a | ||

+ | a `divideByInt` n = a / fromRational (toRational n) | ||

+ | </haskell> | ||

+ | |||

+ | As we can see, recursion cannot deal with average (indirectly), that is why | ||

+ | <math>\mathrm{decrement}\Vert\mathrm{lengthSum}</math> | ||

+ | provides (besides the list with all its element decremented by the delta value) not an average value, but a length---sum pair. | ||

+ | |||

+ | There is no real reason, why the return value of | ||

+ | <math>\mathrm{decrement}\Vert\mathrm{lengthSum}</math> | ||

+ | should be <hask>([a], (a, Integer))</hask> instead of the more simple <hask>([a], a, Integer)</hask>: | ||

+ | let us untangle the nested tuples, using the | ||

+ | <math>\left(A \times B\right) \times C \equiv A \times B \times C</math> | ||

+ | law! | ||

+ | Let us use a | ||

+ | <math>\mathrm{decrement}\Vert\mathrm{sum}\Vert\mathrm{length}</math> | ||

+ | instead of | ||

+ | <math>\mathrm{decrement}\Vert\mathrm{sumLength}</math>: | ||

+ | |||

+ | <haskell> | ||

+ | diff :: Fractional a => [a] -> [a] | ||

+ | diff = trace decrementIIaverage | ||

+ | |||

+ | decrementIIaverage :: Fractional a => [a] -> a -> ([a], a) | ||

+ | decrementIIaverage list delta = let (decremented, sum', length') = decrementIIsumIIlength list delta | ||

+ | in (decremented, sum' `divideByInt` length') | ||

+ | |||

+ | decrementIIsumIIlength :: Fractional a => [a] -> a -> ([a], a, Integer) | ||

+ | decrementIIsumIIlength [] _ = ([], 0, 0) | ||

+ | decrementIIsumIIlength (a : as) delta = let (decremented, sum', length') = decrementIIsumIIlength as delta | ||

+ | in (a - delta : decremented, a + sum', succ length') | ||

+ | |||

+ | divideByInt :: Fractional a => a -> Integer -> a | ||

+ | a `divideByInt` n = a / fromRational (toRational n) | ||

+ | </haskell> | ||

+ | |||

+ | As in the previous problem (normalizing a vector), also at this probem we could imagine other possible solutions: | ||

+ | <haskell> | ||

+ | divdiff :: Fractional a => [a] -> [a] | ||

+ | divdiff = trace decquotIIsumLength | ||

+ | |||

+ | decquotIIsumLength :: Fractional a => [a] -> (a, Integer) -> ([a], (a, Integer)) | ||

+ | decquotIIsumLength list (a, n) = let (dec, sum', length') = decrementIIsumIIlength list (a `divideByInt` n) | ||

+ | in (dec, (sum', length')) | ||

+ | </haskell> | ||

+ | but this is divergent (even for the empty list): it does not terminate at all, albeit this solution seemed for me standing in a duality relation with the previous one. | ||

+ | |||

+ | Excuse me for my toys. | ||

+ | Of course in Wouter Swiestra's [http://www.haskell.org/tmrwiki/WhyAttributeGrammarsMatter Why Attribute Grammars Matter] there are much deeper solutions to this and other examples. | ||

+ | And also very deep connections with | ||

+ | * catamorhisms, category theory | ||

+ | * attribute grammars | ||

+ | are described there. | ||

+ | And there are links from that article to other materials providing deep insights in these powerful theories: attribute grammar tools and articles of the Essential Haskell Compiler project. | ||

+ | |||

+ | It is no wonder that it is exactly a compiler project homepage that provides good circular programming and attribute grammar tutorials | ||

+ | * both on its [http://www.cs.uu.nl/wiki/Ehc/WebHome project page] | ||

+ | * and on a separate [http://www.cs.uu.nl/wiki/HUT/AttributeGrammarSystem attribute grammar] page: | ||

+ | the reason may be that a compiler project is complex enough to require good separation of concerns, modularity, reuse (goals of Aspect Oriented Programming, too), and attribute grammars are powerful tools to achieve these goals. | ||

+ | |||

+ | == Related concepts == | ||

+ | |||

+ | See also [[Attribute grammar]] page. | ||

[[Category:Idioms]] | [[Category:Idioms]] |

## Revision as of 15:04, 5 June 2006

## 1 Porting from HaWiki's CircularProgramming

I port only my material from there (because I know only the copyright status of my writings on HaWiki).

So, the previous materials (to which mine is connected) are on HaWiki's CircularProgramming

-- -- -- THE LOOP OF THE ARROW --- BEING ITSELF ALSO AN ARROW -- -- ,..............................................., -- : : -- : : -- : +---------+ : -- : INPUT |:::::::::| OUTPUT : -- ==========================>+:::::::::+========================> -- : |::ARROW::| : -- : +====>+:::::::::+=====+ : -- : | |:::::::::| | : -- : | +---------+ | : -- : | | : -- : +==========<==========+ : -- : FEEDBACK : -- : : -- ;...............................................; --

class Arrow arrow => ArrowLoop arrow where loop :: arrow (input, feedback) (output, feedback) -> arrow input output

trace :: ((input, feedback) -> (output, feedback)) -> input -> output trace f input = let (output, feedback) = f (input, feedback) in output

(I renamed the type variables and their corresponding argument names, too.)

Of courseinstance ArrowLoop (->) where loop = trace

but a curried function (to an ordered pair):

trace :: (input -> feedback -> (output, feedback)) -> input -> output trace f input = let (output, feedback) = f input feedback in output

For example, for modularizing that magic feedback step in circular programs! By the concept of modularity I mean as described in John Hughes' article Why Functional Programming Matters.

Let us see the tree minimum problem (described above), and let us modularize the magic circular step by this feedback view ofrepmin :: Ord a => Tree a -> Tree a repmin = trace repIImin repIImin :: Ord a => Tree a -> a -> (Tree a, a) repIImin (Leaf minval) rep = (Leaf rep, minval) repIImin (Branch left right) rep = let (left', min_left) = repIImin left rep (right', min_right) = repIImin right rep in (Branch left' right', min min_left min_right)

Notations: in
,
my -`II`

- infix notation means parallelity,
referring to the fact that the wonderful achievement of these functions is being able to work in one traversal.
Sometimes I use -`X`

- infix notation instead of -`II`

-, referring to cross product:

As we can see, we have things enough to confuse: two arguments, a returned tuple, how they depend from each other...

A hint: let us think of the scheme of- input -- what we work with
- output -- what we get
- feedback -- a state-like thing, yes, a feedback, that can be thought of through both aspects
- feedback-as-produced: thinking of it as it leaves the box
- feedback-as-fed-back: thinking of it as it is fed back into the box again

If we decide, what the roles of input, output, feedback plays, then these circular programming funtions will look easier. E.g. here, at the repmin problem

- input is a tree filled with numbers, a heterogenous tree
- output is the repmin'ed tree, a homogenized tree
- feedback (something state-like): here that strange numeric value plays the role of feedback! Its two aspects are
- feedback-as-produced -- in this aspect, it is a minimum value, a
*statistic*of the tree - feedback-as-fed-back -- in this aspect, it is a setter value, a rather
*seed*-like thing: it is spred homogenously in the shape of a given tree

- feedback-as-produced -- in this aspect, it is a minimum value, a

Now it is easier to look at , we have a road-map for our eyes in the jungle of parameters and returned tuple.

-- -- -- How repmin is born: wrapping rep||min with trace -- -- -- ,.........................repmin............................, -- : : -- : : -- : +------------+ : -- : Heterogenous tree |::::::::::::| Homogenized tree : -- =============================>+::::::::::::+================================> -- : |::rep||min::| : -- : +====>+::::::::::::+=====+ : -- : | |::::::::::::| | : -- : setter-val | +------------+ | min-val : -- : a seed | | a statistic : -- : +===========<============+ : -- : FEEDBACK : -- : : -- ;...........................................................; --

normalize :: Floating a => [a] -> [a] normalize = trace divideIInorm divideIInorm :: Floating a => [a] -> a -> ([a], a) divideIInorm vector a = let (divided, norm2) = scaleIInorm2 vector (recip a) in (divided, sqrt norm2) scaleIInorm2 :: Floating a => [a] -> a -> ([a], a) scaleIInorm2 [] _ = ([], 0) scaleIInorm2 (x : xs) a = let (scaled, norm2) = scaleIInorm2 xs a in (a * x : scaled, x * x + norm2)

). We can see what the problem was and how it was solved here:

- requires a rather strict scheme: the function it takes as an argument must be exactly of the same scheme as seen in the above feedback picturetrace
- the recursion by which we can solve the problem eventually (here: math>\mathrm{scale}\Vert\mathrm{norm}^2</math>) usually provides another scheme. E.g. here we cannot compute norm directly by recursion.

we have to convert between these two plugs if they are incompatible, that is why we created .

This is not the only way to force the problem to the shape ofThere are other ways, but they do not seem so natural, they look somewhat artificial, albeit in efficiency aspects, I think, they are the same:

normalize :: Floating a => [a] -> [a] normalize = trace divideBySqrtIInorm2 divideBySqrtIInorm2 :: Floating a => [a] -> a -> ([a], a) divideBySqrtIInorm2 vector scalar = scaleIInorm2 vector (1 / sqrt scalar)

works too, and there is another artificial solution yet:

normalize :: Floating a => [a] -> a -> ([a], a) normalize = trace scaleIIrecipNorm scaleIIrecipNorm :: [a] -> a -> ([a], a) scaleIIrecipNorm vector a = let (scaled, norm2) = scaleIInorm2 vector a in (scaled, recip $ sqrt norm2)

this is good, too. All these three solutions converge, they terminate.

But there are circular programming examples where such possible solutions will not be so equivalent: one of them may diverge, while the other converge. It seems to me, the diff example (described in Wouter Swiestra's Why Attribute Grammars Matter) is exactly such circular programming problem.

This diff example is:

- how to compute the average of a list
- and (in the same traversal of the list) decrease all elements of the list by this average value.

Here, the gap between

- the needs of `trace`
- and possibilities of (single traversal) recursion

seems even wider than in the previous problem. Let us try:

diff :: Fractional a => [a] -> [a] diff = trace decrementIIaverage decrementIIaverage :: Fractional a => [a] -> a -> ([a], a) decrementIIaverage list delta = let (decremented, (sum', length')) = decrementIIsumLength list delta in (decremented, sum' `divideByInt` length') decrementIIsumLength :: Fractional a => [a] -> a -> ([a], (a, Integer)) decrementIIsumLength [] _ = ([], (0, 0)) decrementIIsumLength (a : as) delta = let (decremented, (sum', length')) = decrementIIsumLength as delta in (a - delta : decremented, (a + sum', succ length')) divideByInt :: Fractional a => a -> Integer -> a a `divideByInt` n = a / fromRational (toRational n)

As we can see, recursion cannot deal with average (indirectly), that is why provides (besides the list with all its element decremented by the delta value) not an average value, but a length---sum pair.

There is no real reason, why the return value of

should belet us untangle the nested tuples, using the law! Let us use a instead of :

diff :: Fractional a => [a] -> [a] diff = trace decrementIIaverage decrementIIaverage :: Fractional a => [a] -> a -> ([a], a) decrementIIaverage list delta = let (decremented, sum', length') = decrementIIsumIIlength list delta in (decremented, sum' `divideByInt` length') decrementIIsumIIlength :: Fractional a => [a] -> a -> ([a], a, Integer) decrementIIsumIIlength [] _ = ([], 0, 0) decrementIIsumIIlength (a : as) delta = let (decremented, sum', length') = decrementIIsumIIlength as delta in (a - delta : decremented, a + sum', succ length') divideByInt :: Fractional a => a -> Integer -> a a `divideByInt` n = a / fromRational (toRational n)

As in the previous problem (normalizing a vector), also at this probem we could imagine other possible solutions:

divdiff :: Fractional a => [a] -> [a] divdiff = trace decquotIIsumLength decquotIIsumLength :: Fractional a => [a] -> (a, Integer) -> ([a], (a, Integer)) decquotIIsumLength list (a, n) = let (dec, sum', length') = decrementIIsumIIlength list (a `divideByInt` n) in (dec, (sum', length'))

but this is divergent (even for the empty list): it does not terminate at all, albeit this solution seemed for me standing in a duality relation with the previous one.

Excuse me for my toys. Of course in Wouter Swiestra's Why Attribute Grammars Matter there are much deeper solutions to this and other examples. And also very deep connections with

- catamorhisms, category theory
- attribute grammars

are described there. And there are links from that article to other materials providing deep insights in these powerful theories: attribute grammar tools and articles of the Essential Haskell Compiler project.

It is no wonder that it is exactly a compiler project homepage that provides good circular programming and attribute grammar tutorials

- both on its project page
- and on a separate attribute grammar page:

the reason may be that a compiler project is complex enough to require good separation of concerns, modularity, reuse (goals of Aspect Oriented Programming, too), and attribute grammars are powerful tools to achieve these goals.

## 2 Related concepts

See also Attribute grammar page.