HOME
PDF (letter size)
PDF (legal size)

my Ada corner

Nasser M. Abbasi

August 28, 2015   Compiled on May 23, 2020 at 8:57am

Contents

1 Original Ada web site and Ada links
2 how to call Lapack and Blas directly from Ada
 2.1 Introduction
 2.2 Review of the LAPACK and BLAS Ada binding
 2.3 Installation instruction
 2.4 source code
3 How to compile GTK Ada program
4 How to make Ada generate an exception on some floating points operations?
5 How to use Ada 2005 OO?
6 how to make simple Ada program
7 Ada implementation of decimal representation of exp()
8 Ada implementation of getopt()
9 GNAT 2012 installation log file

1 Original Ada web site and Ada links

My original Ada web site with links is here

2 how to call Lapack and Blas directly from Ada

Important note added June 2013:

This page is obsolete now, left here for archive and information only.

The Ada Lapack code is now housed at http://sourceforge.net/projects/ada-lapack/

2.1 Introduction

I downloaded the original BLAS Ada binding written by Duncan Sands from http://topo.math.u-psud.fr/~sands/Programs/BLAS/index.html and the LAPACK Ada binding written by Wasu Chaopanon from ftp://ftp.cs.kuleuven.be/pub/Ada-Belgium/mirrors/gnu-ada/OLD/contrib/lapack-ada/

And made some minor improvments to the bindings.

This page describes the minor changes made and instructions how to use these bindings from Ada in order call LAPACK and BLAS Fortran functions.

A new tar file for LAPACK and for BLAS with all the changes can be downloaded from the link below.

The changes made to LAPACK binding involve streamlining the source tree structure, writing new Makefiles, simplify the binding to use one package called lapack and also adding the documentation shown below.

Changes for the BLAS binding were minimal. It involved changes to the source tree structure and writing Makefiles and adding the documentation shown below.

2.2 Review of the LAPACK and BLAS Ada binding

The Ada binding to LAPACK and BLAS is an Ada package which allows one to interface to the native lapack and blas libraries.

The native lapack and blas libraries need to be first installed on the system (on Linux, these libraries will normally be found in /usr/lib/liblapack.so and /usr/lib/libblas.so )

To use LAPACK from Ada, one needs to install both the native LAPACK and BLAS libraries since LAPACK depends on BLAS.

The Ada binding is a thin binding, meaning there is 1-1 mapping between the call to the Ada routine and the corresponding Fortran routine using the same function name in the Fortran libraries.

The following diagram illustrates the use of the LAPACK binding with the needed gnatmake command to compile and link the client Ada program.

pict
Figure 1:high level1 lapack

And a similar diagram for the BLAS binding interface

pict
Figure 2:high level1 blas

The source tree structure for LAPACK is described in this diagram

pict
Figure 3:tree structure for lapack

The full content of the LAPACK tree is listed here lapack_tree_listing.txt

The source tree structure for the BLAS binding is described in this diagram

pict
Figure 4:tree structure for blas

2.3 Installation instruction

These instructions explain how to use the Ada LAPACK and BLAS binding in the updated snapshot tar file.

  1. install native lapack 77 and blas libraries on your system. These will normally be installed in /usr/lib/.
  2. download the updated Ada binding in the zip files given in the links at the bottom of this page. They are ada_lapack_073012.zip and ada_blas_073012.zip
  3. Extract the zip files to some location. This will create 2 source trees as shown in the diagrams above.
  4. At the top of the each source tree, there is a file called common.mk where a Makefile variable is set to point to the directory that contains the native LAPACK and BLAS libraries. This is currently defined to point to /usr/lib. Edit this line to change this location only if the location is different in your system.
  5. To build the binding, just type make from the top of each tree. Make will build the whole tree, including the bindings packages and the test programs.
  6. To run the test program, type make testing from the top of tree for LAPACK and for BLAS.
  7. Examples of clients using the bindings are found in the tests/ directory of each tree.
  8. The binding is in the binding/ directory in each tree. This is the directory that you need to set the -I option to point to when using gnatmake as shown in the diagrams above.
  9. The following is a simple example of using the Ada LAPACK binding to solve \(Ax=b\) mysolve.adb.txt which can also be found in the ada_lapack/tests/pilot/ directory
  10. That is all! Have fun using Ada with LAPACK and BLAS.

2.4 source code

  1. ada_lapack_073012.zip
  2. ada_blas_073012.zip

3 How to compile GTK Ada program

$ gnatmake -I../pragmarc mine_detector.adb `gtkada-config`
gcc-4.4 -c -I../pragmarc -I/usr/share/ada/adainclude/gtkada2 mine_detector.adb
gcc-4.4 -c -I../pragmarc -I/usr/share/ada/adainclude/gtkada2 user_if.adb
gcc-4.4 -c -I../pragmarc -I/usr/share/ada/adainclude/gtkada2 field.ads
gcc-4.4 -c -I../pragmarc -I/usr/share/ada/adainclude/gtkada2 field-operations.adb

gnatbind -I../pragmarc -aI/usr/share/ada/adainclude/gtkada2 -aO/usr/lib/ada/adalib/gtkada2 -x mine_detector.ali

gnatlink mine_detector.ali -L/usr/lib -lgtkada2 -pthread -lgtk-x11-2.0 -lgdk-x11-2.0 -latk-1.0 -lgio-2.0
-lpangoft2-1.0 -lpangocairo-1.0 -lgdk_pixbuf-2.0 -lm -lcairo -lpango-1.0 -lfreetype
-lfontconfig -lgobject-2.0 -lgmodule-2.0 -lgthread-2.0 -lrt -lglib-2.0
$

4 How to make Ada generate an exception on some floating points operations?

On Thu, 20 Nov 2008 12:09:41 +0100, Markus Schoepflin wrote:

> is it possible to influence the behaviour of GNAT regarding the handling of
> NANs? (Most importantly in the special case of division by zero.)
>
> We need to get exceptions whenever a NAN is generated, is this possible
> somehow? (For example by setting Machine_Overflow to True and recompiling
> the compiler itself.)

You can scrap IEEE stuff in favor of Ada semantics by declaring your own
floating-point [sub]type with a range specified. The compiler will be
forced to check values:

   type Safe_Float is digits 6 range -10.0E10..+10.0E10;

or

   subtype Safe_Float is Float range Float'Range;

then

   X : Safe_Float := 1.0;
   Y : Safe_Float := 0.0;
begin
   Y := X / Y;
exception
   when Error : others =>  -- Should print "range check failed"
      Put_Line (Exception_Message (Error));
end;

--
Regards,
Dmitry A. Kazakov
http://www.dmitry-kazakov.de

5 How to use Ada 2005 OO?

from http://en.wikibooks.org/wiki/Ada_Programming/Object_Orientation

6 how to make simple Ada program

If gnat is not installed, install it (on linux) using something similar to

sudo apt-get install gnat-4.6

write the following code in file called hello_world.adb

compile using

gnatmake hello_world.adb

         gcc-4.6 -c hello_world.adb
         gnatbind -x hello_world.ali
         gnatlink hello_world.ali

Run it using ./hello_world

7 Ada implementation of decimal representation of exp()

This is an Ada implementation of decimal representation of \(e\) based on SPIGOT algorithm for \(\pi \) by S. Rabinowitz & S. Wagon, The American Mathematical Monthly, March 1995.

This source file is here

8 Ada implementation of getopt()

This package is an Ada implementation of getopt() as specified by the document ”The Single UNIX Specification, Version 2”, Copyright 1997 The Open Group

This source file is here

This is a test program of the above package. source code

9 GNAT 2012 installation log file

gnat2012_installation_log_file.txt