Hey, what's going on?

Archive for January, 2011

Installing LAPACK/BLAS on Mandriva 2010

Posted by Syeilendra Pramuditya on January 19, 2011

Download the source file from the web. I downloaded version 3.3.0 (lapack-3.3.0.tgz). Save it to /usr/local/ (you need to use root account). Make sure (change if necessary) that the “owner” and the “group” are your linux user name.

  1. By using root account, extract the file: gunzip -c lapack-3.30.tgz | tar xvf –
  2. Type: cd lapack-3.3.0
  3. Type: chown yourusername *.* * -R
  4. Type: chgrp yourusername *.* * -R
  5. Exit the root account and just use your user account
  6. Type: cp make.inc.example make.inc (by default, gfortran compiler will be used, otherwise you need to edit this make.inc file)
  7. Type: make blaslib
  8. Type: make lapacklib
  9. Type: cp blas_LINUX.a libblas.a
  10. Type: cp lapack_LINUX.a liblapack.a
  11. Finish

Now you need to run a sample case, I took the following code from Ref. 3:

c This code is written in Fortran
      Program LinearEquations
c solving the matrix equation A*x=b using LAPACK
      Implicit none
c declarations, notice single precision
      Real*4 A(3,3), b(3)
      integer i, j, pivot(3), ok
c define matrix A
      A(1,1)=3.1
      A(1,2)=1.3
      A(1,3)=-5.7
      A(2,1)=1.0
      A(2,2)=-6.9
      A(2,3)=5.8
      A(3,1)=3.4
      A(3,2)=7.2
      A(3,3)=-8.8
c define vector b, make b a matrix and you can solve multiple
c equations with the same A but different b
      b(1)=-1.3
      b(2)=-0.1
      b(3)=1.8
c find the solution using the LAPACK routine SGESV
      call SGESV(3, 1, A, 3, pivot, b, 3, ok)
c
c parameters in the order as they appear in the function call
c order of matrix A, number of right hand sides (b), matrix A,
c leading dimension of A, array that records pivoting,
c result vector b on entry, x on exit, leading dimension of b
c return value
c
c print the vector x
      do i=1, 3
      write(*,*) b(i)
       end do
      end

To compile and link, type: gfortran -o main.unix main.f -L/usr/local/lapack-3.3.0 -llapack -lblas

To run type: ./main.unix

The answer should be [1.0 , 1.0 , 1.0]

References:

  1. NETLIB
  2. LLNL
  3. ORST
  4. UTAH

Posted in software & simulation | Tagged: , , , , , | Leave a Comment »

NECTAR Code – A Steady State Thermal Hydraulics Solver

Posted by Syeilendra Pramuditya on January 4, 2011

Shortlink: http://wp.me/p61TQ-zt

Software License Agreement

This software is given “as is”. The author holds its copy rights, one can download, copy, use, modify, and spread it free of charge, provided that the original sources (the author and this webpage) are properly cited. The author assumes no responsibility, implicit or explicit, for any possible loss or damage due to the use of this software.

By downloading the software, you agree to accept this agreement.


Download code package:

Preview:

Posted in nuclear engineering, software & simulation | Leave a Comment »