javascript - Use Emscripten with Fortran: LAPACK binding - Stack Overflow

admin管理员组

文章数量:1333386

My goal is to use LAPACK with Emscripten.

My question is: how to port LAPACK to JS? The are two ways I can think of: CLAPACK to JS where my question is: does anybody know an unofficial version that is later than 3.2.1? And the other way to think of is: how to port FORTRAN to JS?

Emscripten is capable of transforming C code to JavaScript. But unfortunately, LAPACK 3.5.0 (/) is only available in FORTRAN95.

The CLAPACK project (/) is basically what I want: a C version of LAPACK. But this one is outdated; the latest is 3.2.1.

F2C only works up to FORTRAN 77. LAPACK 3.5.0 was written in FORTRAN 95.

So my question now is: why is there no newer port of LAPACK to C?

The optimal way would be to directly transform the FORTRAN95 code of LAPACK to javascript with clang and emscripten. But I just don't know where to start.

Emscripten currently does not support FORTRAN. But it handles LLVM bitcode, so it should not be a problem to use clang to generate LLVM bc from a FORTRAN file.

For testing purpose, I have this file:

      program hello
      print *, "Hello World!"
      end program hello

It piles just fine with "clang hello.f -o hello -lgfortran". I am not capable of transforming this into valid bitcode.

clang -c -emit-llvm hello.f      
clang -S -emit-llvm hello.f -o hello.bc -lgfortran

None of these approaches works, because emscripten keeps telling me

emcc -c hello.o -o hello.js
hello.o is not valid LLVM bitcode

I am not sure anyways if this would be even possible, because LAPACK obviously needs libgfortran to work. And I can't merge a library into javascript code...

Thanks in advance!

Edit:

I almost managed it to convert BLAS from LAPACK 3.5.0 to JS. I used dragonegg to acplish this.

gfortran caxpy.f -flto -S -fplugin=/usr/lib/gcc/x86_64-linux-gnu/4.6/plugin/dragonegg.so 
gfortran cgerc.f ...
...

After gaining LLVM bitcode from that:

emcc caxpy.s.ll cgerc.s.ll cher.s.ll ... -o blas.js -s EXPORTED_FUNCTIONS="['_caxpy_', ... , '_ztpsv_']"

But emscripten still leaves me with the following errors:

warning: unresolved symbol: _gfortran_st_write
warning: unresolved symbol: _gfortran_string_len_trim
warning: unresolved symbol: _gfortran_transfer_character_write
warning: unresolved symbol: _gfortran_transfer_integer_write
warning: unresolved symbol: _gfortran_st_write_done
warning: unresolved symbol: _gfortran_stop_string
warning: unresolved symbol: cabs
warning: unresolved symbol: cabsf

AssertionError: Did not receive forwarded data in an output - process failed?

The problem is that lgfortran is prepiled I think.

My goal is to use LAPACK with Emscripten.

My question is: how to port LAPACK to JS? The are two ways I can think of: CLAPACK to JS where my question is: does anybody know an unofficial version that is later than 3.2.1? And the other way to think of is: how to port FORTRAN to JS?

Emscripten is capable of transforming C code to JavaScript. But unfortunately, LAPACK 3.5.0 (http://wwwlib/lapack/) is only available in FORTRAN95.

The CLAPACK project (http://wwwlib/clapack/) is basically what I want: a C version of LAPACK. But this one is outdated; the latest is 3.2.1.

F2C only works up to FORTRAN 77. LAPACK 3.5.0 was written in FORTRAN 95.

So my question now is: why is there no newer port of LAPACK to C?

The optimal way would be to directly transform the FORTRAN95 code of LAPACK to javascript with clang and emscripten. But I just don't know where to start.

Emscripten currently does not support FORTRAN. But it handles LLVM bitcode, so it should not be a problem to use clang to generate LLVM bc from a FORTRAN file.

For testing purpose, I have this file:

      program hello
      print *, "Hello World!"
      end program hello

It piles just fine with "clang hello.f -o hello -lgfortran". I am not capable of transforming this into valid bitcode.

clang -c -emit-llvm hello.f      
clang -S -emit-llvm hello.f -o hello.bc -lgfortran

None of these approaches works, because emscripten keeps telling me

emcc -c hello.o -o hello.js
hello.o is not valid LLVM bitcode

I am not sure anyways if this would be even possible, because LAPACK obviously needs libgfortran to work. And I can't merge a library into javascript code...

Thanks in advance!

Edit:

I almost managed it to convert BLAS from LAPACK 3.5.0 to JS. I used dragonegg to acplish this.

gfortran caxpy.f -flto -S -fplugin=/usr/lib/gcc/x86_64-linux-gnu/4.6/plugin/dragonegg.so 
gfortran cgerc.f ...
...

After gaining LLVM bitcode from that:

emcc caxpy.s.ll cgerc.s.ll cher.s.ll ... -o blas.js -s EXPORTED_FUNCTIONS="['_caxpy_', ... , '_ztpsv_']"

But emscripten still leaves me with the following errors:

warning: unresolved symbol: _gfortran_st_write
warning: unresolved symbol: _gfortran_string_len_trim
warning: unresolved symbol: _gfortran_transfer_character_write
warning: unresolved symbol: _gfortran_transfer_integer_write
warning: unresolved symbol: _gfortran_st_write_done
warning: unresolved symbol: _gfortran_stop_string
warning: unresolved symbol: cabs
warning: unresolved symbol: cabsf

AssertionError: Did not receive forwarded data in an output - process failed?

The problem is that lgfortran is prepiled I think.

• javascript
• c
• fortran
• lapack
• emscripten

Share Improve this question Follow edited Feb 24, 2014 at 15:20 japedo asked Feb 24, 2014 at 14:13 japedojapedo 39033 silver badges1616 bronze badges 10

• CLAPACK says "f2c'ed version of LAPACK" in the title. I haven't tried this, but why don't you run f2c on LPACK 3.5.0? – Carsten Commented Feb 24, 2014 at 14:16
• 1 You write So my question now is: why is there no newer port of LAPACK to C? but that is not, any longer, your question. That line is a hang over from your earlier version of your question and has already been answered. I can't see, anywhere in this version of your question, the statement that f2c only works up to FORTRAN77; your response to @Carsten's ment was a bit snippy. – High Performance Mark Commented Feb 24, 2014 at 14:40
• 1 You are right. I did not mention that again, that is a copy and paste error from the previous question. Sorry Carsten ;) My question is: how to port LAPACK to JS? The are two ways I can think of: CLAPACK to JS where my question is: does anybody know an unofficial version that is later than 3.2.1? And the other way to think of is: how to port FORTRAN to JS? – japedo Commented Feb 24, 2014 at 14:56
• I'm wondering, what you final goal is. Do you want to write a numbercrunching app in JS? Sounds like a bad idea to me... – Stefan Commented Feb 24, 2014 at 17:53
• 1 Did you make progress on this? Did you see this github project for getting the sources to gfortran and piling it. github./alexleach/gfortran-mlion. There is a SO thread discussing this here: stackoverflow./questions/12316780/… – Trausti Kristjansson Commented Apr 23, 2014 at 1:29

 |  Show 5 more ments

2 Answers 2

Sorted by: Reset to default 6

Thank you for your reply! Indeed I did make progress on this. Finally it is working. I was very close, just follow these steps:

gfortran caxpy.f -S -flto -m32 -fplugin=dragonegg.so
mv caxpy.s caxpy.ll
llvm-as caxpy.ll -o caxpy.o

Note the "m32" flag which I missed earlier.

Warnings like

warning: unresolved symbol: _gfortran_st_write

can be ignored safely. Emscripten creates empty functions in the JavaScript file with this name so if these functions are not called at all there is no problem. If they get called you can easily substitute them with your own functions; the names are somewhat descriptive. Additionally you can have a look at the libgfortran source code (be aware it is GPL).

With this Emscripten source can be extended by hand to support Fortran files. Someday I may publish this on github!

I actually pulled this off recently (https://github./harveywi/arpack-js). The Github repo is mostly barren except for the output JS files, but I will be uploading source code, Makefiles, and other instructions soon. After wrangling unsuccessfully with dragonegg for a while, I found a different approach (not as great, but sufficient) which did the trick.

Here is roughly how I did it:

1. Download the ARPACK source code.
2. Run f2c on all of the Fortran files to convert them to C.
3. (This might be the trickiest part): Modify the Makefiles to use Emscripten and LLVM toolchains.
4. Make the project to produce an LLVM binary.
5. Use Emscripten again to transpile the LLVM binary to JS.

本文标签： javascriptUse Emscripten with Fortran LAPACK bindingStack Overflow