Its main trick is intense (and super expensive) control flow analysis to aggressively monomorphise and unbox values.
i did not realize i needed to have something like that in my life
Of course, in Futhark it should be called the thorn ( ᚦ) operator.
2 replies
That's why I've started using
& in Haskell! It's probably not idiomatic, but it is nicer than $ imo.
I did some preliminary testing for the lexer VS flex, flex reaches about 150-180 MB/s depending on the compiler, on my Ryzen 3700X. Running futhark on a RTX 2080 TI yielded about 5 GB/s, though i cant actually double check since my uni broke that machine. On my RX 580 it does about 650 MB/s, and futhark c does about 80 MB/s.
The latter can probably be explained by the large lookup table i use, which isn't very cache friendly of course.
I don't think its physically broken, but my faculty has been messing with these machines ever since the start of this academic year.
1 reply
I haven't tried other flex alternatives, although the name (fast lexical analyzer generator) would lead me to believe that it should be decently competitive. Its gnu software though, so its methods may be dated by now.
This isn't even a cluster system. There are just a few machines you can use as student, and another few you can use as staff. All fair-use, although that wont stop students from running R programms that run for several months and require 700GB of memory <grumble>
Up until recently students could log into staff machines so i had just been using those.
Every once in a while a machine breaks magically and usually its fixed within a few hours but not this machine i guess
Is it intentional that
(a).b is an invalid module expression? I guess its kind of weird, but i'm trying to import a parametric module and apply it simultaneously
so i have something like
a.fut:
module a (T: integer) = {
let add (a: T.t) (b: T.t): T.t = a + b
}and then in b.fut:
module a_u8 = (import "a").a u8(x).y is a special expression in the term language (because x.y is handled as a single lexeme), and I guess we just never added it similarly to the module language.
Just curious; what does the 'bikeshed' label mean? It's only on the AD issue as far as I can tell?
Based on this: http://bikeshed.org/
And regarding mailing lists, I sure wish they had that functionality sometimes
Second link does not work for me
I shall check the archives
Hello everyone :). I'm new to GPU programming and futhark but I want to try to speed up a rendering process of points. And would like some input on how it could be handled. I have a large set of 2d points where each coordinate is a float32. I want to "rasterize" these points and then count how many of them are visible. I understand I essentially have to implement a rasterization algorithm. However I'm wondering if this can be done efficiently at all on a GPU. Essentially I seem to have 2 choices. Either loop through every point and set the correspondin cell to 1. Or loop through the grid and check whether each cell contains at least a single point. But option 1 has concurrent write access to the grid which afaik is bad. And option 2 is slow because it's hard to determine whether a point falls into a cell.
Are there any other ways to attack this?
In fact I don't really need a final image. I just want to see have the number of cells which contain at least a single point.
@rowanG077 GPUs can indeed be used for graphics. Some of them even allow you to connect a monitor! Most modern 3D is based rasterisation, too. Futhark doesn't expose most of the graphics facilities provided by GPUs, but people have written various kinds of rasterisers. One of my colleagues wrote a paper with some techniques: https://futhark-lang.org/publications/array19.pdf
I think your option 1 is best. Concurrent write accesses are not great, but in Futhark it would be a "generalised histogram" (
reduce_by_index), for which Futhark has a high-quality implementation.
I am using reduce_by_index to update cells in the neighborhoods around many cells. I am trying to generalize to allow the user to pass in the radius of the neighborhood, but I think I am stuck because I need more expressive shape types.
Right now, the radius is hardcoded to 1, so we have 9 cells in the neighborhood, and my type is
neigborhood (i: i64) (j: i64): [9](i64, i64)
I would like to write
neigborhood (i: i64) (j: i64) (radius: i64) : [(2*radius + 1)**2](i64, i64)
which is currently impossible. Is there a workaround for this? Radius would be ultimately passed in at my entry point.
Can I make it so that I can use Futhark code as a library that will first attempt to use OpenCL and then fall back to sequential C code if OpenCL cannot be used?
I'm mostly worried about link failures