I made a new port "lapack++"
This facilitate quite awkward usage of lapack, by
making use of the C++ class. For example, diagonalization of
matrices can be achived in Fortran is quite troublesome:

call dgeev('N','V',a,b,lda,wr,wi,vl,ldvl,vr,ldvr
     1     ,work,lwork,info)
a is dimension of matrix, b is a matrix which diagonalized
'N', 'V' are some spell ;-) that controlls the LAPACK subroutine,
lda, wr, wi, vl... are... 
However, using Lapack++, it can be shortened:
LaEigSolve(A, v, B);
A, B are n times n matrices, and v is n dimension vector.
Very Very easy!!

It is worth to include in FreeBSD ports tree!