TAPENADE Example

TAPENADE is a Automatic Differentiation (AD) tool for FORTRAN. It generates differentiation code by source code transformation.

For illustration of TAPENADE usage, following test expression is selected

This expression/function is coded in FORTRAN as a subroutine which takes 3 arguments whose in/out behaviour is stated by intent directives

! FILE : test1.f90
subroutine test1(x,y,z)
	real, intent(in) :: x,y
	real, intent(out) :: z

	z = sin(x+2*y)*cos(x-sin(y));
end subroutine test1

TAPENADE executed to differentiate test1 subroutine (in tangent mode) generates following differentiation code

!        Generated by TAPENADE     (INRIA, Tropics team)
!  Tapenade 3.3 (r3163) - 09/25/2009 09:03
!
!  Differentiation of test1 in forward (tangent) mode:
!   variations  of output variables: z
!   with respect to input variables: x y
SUBROUTINE TEST1_D(x, xd, y, yd, z, zd)
  IMPLICIT NONE
  REAL, INTENT(IN) :: x, y
  REAL, INTENT(IN) :: xd, yd
  REAL, INTENT(OUT) :: z
  REAL, INTENT(OUT) :: zd
  REAL :: arg1
  REAL :: arg1d
  INTRINSIC COS
  INTRINSIC SIN
  arg1d = xd - yd*COS(y)
  arg1 = x - SIN(y)
  zd = (xd+2*yd)*COS(x+2*y)*COS(arg1) - SIN(x+2*y)*arg1d*SIN(arg1)
  z = SIN(x+2*y)*COS(arg1)
END SUBROUTINE TEST1_D

Note that each input and output variable that user specified caused insertion of “differentials” of those variables, that is, the input code

subroutine test1(x,y,z)

results in

SUBROUTINE TEST1_D(x, xd, y, yd, z, zd)

In this example,

are inserted just after the original argument.

Actual differential of test1 function is (total differential)

with

Note that in order to compute partial derivatives and , the generated derivative routine test1_d must be called twice

! derivative computation point
x = ...
y = ...
xd = 1.0
yd = 0.0
CALL TEST1_D(x, xd, y, yd, z, zd) ! z and zd ~ dz/dx is obtained
xd = 0.0
yd = 1.0
CALL TEST1_D(x, xd, y, yd, z, zd) ! z and zd ~ dz/dy is obtained

In order to avoid such multiple calls, TAPENADE has an alternate differentiation mode called “Tangent Vectorial Mode”

If TAPENADE code is generated with this option then following differentiation code is obtained

!        Generated by TAPENADE     (INRIA, Tropics team)
!  Tapenade 3.3 (r3163) - 09/25/2009 09:03
!
!  Differentiation of test1 in forward (tangent) mode: (multi-directional mode)
!   variations  of output variables: z
!   with respect to input variables: x y
SUBROUTINE TEST1_DV(x, xd, y, yd, z, zd, nbdirs)
  USE DIFFSIZES
!  Hint: nbdirsmax should be the maximum number of differentiation directions
  IMPLICIT NONE
  REAL, INTENT(IN) :: x, y
  REAL, INTENT(IN) :: xd(nbdirsmax), yd(nbdirsmax)
  REAL, INTENT(OUT) :: z
  REAL, INTENT(OUT) :: zd(nbdirsmax)
  REAL :: arg1
  REAL :: arg1d(nbdirsmax)
  REAL :: arg2
  REAL :: arg2d(nbdirsmax)
  INTEGER :: nd
  INTEGER :: nbdirs
  INTRINSIC COS
  INTRINSIC SIN
  arg1 = x + 2*y
  arg2 = x - SIN(y)
  DO nd=1,nbdirs
    arg1d(nd) = xd(nd) + 2*yd(nd)
    arg2d(nd) = xd(nd) - yd(nd)*COS(y)
    zd(nd) = arg1d(nd)*COS(arg1)*COS(arg2) - SIN(arg1)*arg2d(nd)*SIN(&
&      arg2)
  END DO
  z = SIN(arg1)*COS(arg2)
END SUBROUTINE TEST1_DV

Note that generated code executes differentiation code for all given differentiation directions inside a loop. This, in effect, corresponds to calling “single-direction” mode successively, but with the direction loop inside the differentiation code, subroutine arguments are not passed to calls again and again, increasing the efficiency of code. Another benefit of this mode is hiding details of successive calls.

Following (incomplete) program illustrates use of “multi-directional” mode.

integer, parameter :: nbdirs = 8
real*8 xd(nbdirs),yd(nbdirs)
xd = (/1.0, 0.0/)
yd = (/0.0, 1.0/)
call TEST1_DV(x, xd, y, yd, z, zd, nbdirs) ! zd contains all partial derivatives

TAPENADE GUI Usage

First click add button

Select the source file(s) test.f90

In listbox component all source code added is shown

Click Parse Button. On the left drop-down list box, subroutines/functions available in parsed source codes will appear

In input variables, click Add and select x,y

Pressing OK will add selected variable(s) to input variables listbox component

In output variable, click Add and select z

Pressing OK will add selected variable(s) to output variables listbox component

Finally select dierentiation mode and click dierentiate. This will generate code and automatically open html report page (if browser is properly dened in environment)

The Levenshtein distance function returns how many edits (deletions, insertions, transposition) are required to turn one string into another.  Here’s my attempt to write that function in Fortran.  It seems to work.  Please chime in with improvements.

integer function distance (a,b)

implicit none

! define variables

integer :: d, len_a, len_b, i, j, cost
character(len=*), intent(in) :: a, b

! matrix for calculating levenshtein distance
integer, dimension(0:len_trim(a), 0:len_trim(b)) :: leven_mat

integer, dimension(3) :: three_vals

len_a = len_trim(a)
len_b = len_trim(b)

do i = 0,len_a
leven_mat(i,0) = i
end do

do j = 0, len_b
leven_mat(0,j) = j
end do

do i = 1, len_a
do j = 1, len_b
if (a(i:i) == b(j:j)) then
cost = 0
else
cost = 1
end if
three_vals(1) = leven_mat(i-1,j) + 1
three_vals(2) = leven_mat(i, j-1) + 1
three_vals(3) = leven_mat(i-1,j-1) + cost
leven_mat(i,j) = minval(three_vals)
end do
end do

distance = int(leven_mat(len_a,len_b))

end function distance

To install the g77 Fortran compiler in Ubuntu 9.10 (Karmic Koala) just follow these simple steps.

Open /etc/apt/sources.list and add the following Ubuntu 8.04 Hardy repositories:

deb http://hu.archive.ubuntu.com/ubuntu/ hardy universe
deb-src http://hu.archive.ubuntu.com/ubuntu/ hardy universe
deb http://hu.archive.ubuntu.com/ubuntu/ hardy-updates universe
deb-src http://hu.archive.ubuntu.com/ubuntu/ hardy-updates universe

Once you’ve saved the file open a terminal and run the commands:

sudo aptitude update
sudo aptitude install g77

After g77 is installed, comment out or remove the lines above from /etc/apt/sources.list and rerun:

sudo aptitude update

The instructions in this post were largely taken from install g77 in Ubuntu 9.04 (Jaunty Jackalope)

418214     การเขียนโปรแกรมภาษาฟอร์แทรน     FORTRAN Programming

โครงสร้างและองค์ประกอบของภาษาฟอร์แทรน หลักการเขียนโปรแกรมภาษาฟอร์แทรนและการประยุกต์สำหรับงานด้านต่างๆ

(Structure and elements of FORTRAN. Principles of programming in FORTRAN and its applications.)

(418214 มหาวิทยาลัยเกษตรศาสตร์)

to get the current time in Fortran, we can use “system_clock” built-in function. To do this, take a look at the code below.

integer::count_0,count_rate,count_max, walltime, tstart,tend
call system_clock(count_0,count_rate,count_max)
tstart=count_0/count_rate
write(*,*) 'begin (system_clock)',tstart


To find out how long your program has run, add a similar line and save the time to, for example “t_end” variable, and get the difference. It will give the duration time in seconds.

I’ve uploaded the previous post’s source code to github.
Get it at  http://gist.github.com/226473

Genap dua kali sudah, saya sedikit berkecimpung dalam membantu mata kuliah Bahasa Pemrograman di Teknik Sipil UGM. Dua periode itu, membuat saya sedikit banyak mengetahui permasalahan yang kerap dihadapi mahasiswa baru di semester I, yang begitu beberapa bulan masuk kuliah langsung dijejali dengan logika-logika pemrograman. Bagi mahasiswa yang telah mendapatkan pelajaran TIK (Teknologi Informasi dan Komunikasi) semasa SMA, hampir dipastikan akan segera menyerap logika-logika tersebut dengan lancar.

Sementara bagi mereka yang kurang beruntung, dan tidak mendapatkan dasar-dasar materi Pemrograman semasa SMA, harus berjuang sedikit lebih keras untuk bisa memahami logika-logika percabangan, perulangan (looping), algoritma, flowchart, dan tetek bengek pemrograman lainnya.

Bagi mereka yang kritis, biasanya ada dua pertanyaan seputar mata kuliah Bahasa Pemrograman yang kerap ditanyakan :

1. Mengapa di jurusan seperti Teknik Sipil harus belajar Bahasa Pemrograman ?
2. Mengapa harus mempelajari bahasa FORTRAN ?

Untuk pertanyaan yang pertama, biasanya akan segera terjawab dengan sendirinya jika telah memasuki perkuliahan di semester-semester 4 dan 5. Terutama ketika berhadapan dengan kasus-kasus Hidrolika, peninggian muka air yang harus diselesaikan dengan selembar penuh hitungan iterasi, penyelesaian debit perpipaan dengan sistem looping, perancangan dimensi saluran drainase perkotaan dengan software DUFLOW, dan juga yang paling legendaris dalam hal struktur, penyelesaian dengan Finite Element Method dan juga penyelesaian dengan Metode Matriks. Ruwetnya hitungan-hitungan tersebut akan dapat diselesaikan dengan program sederhana, yang dapat dibuat di kalkulator programmable sekalipun. Belum jika kita menemukan kasus-kasus lokal dalam proyek, yang memaksa kita untuk membuang software yang telah ada, dan membuat sendiri software yang sesuai dengan kasus tersebut.

Sementara untuk pertanyaan yang kedua, mengapa harus FORTRAN ? Perlu diketahui, biasanya dalam tengah semester awal, mata kuliah Bahasa Pemrograman murni mengajarkan program FORTRAN. Sebuah bahasa yang boleh dibilang kuno, yang dibuat di periode 1960-an saat awal munculnya komputer IBM. Jenis FORTRAN yang digunakan adalah WATFOR77, versi FORTRAN yang sesuai dengan namanya, dirilis pada tahun 1977. Lebih dari 30 tahun lalu.

Bahasa FORTRAN, ketika pertama kali mempelajarinya, terlihat sebagai bahasa yang lebih dulu diciptakan daripada C++, Pascal, dan bahasa pemrograman lainnya. Hal ini terasa dari beberapa poin :

1. Tidak adanya operator “>” (lebih dari), “<” (kurang dari) dan “=” (sama dengan) pada perintah IF. Sebagai gantinya, digunakan LT (Less Than “<”), GT (Bigger Than “>”), dan EQ (Equal “=”)
2. Sangat Case Sensitive, terutama ketika kita menyimpan file, harus dengan huruf kapital, jika tidak, jangan harap bisa dibaca oleh Compiler.
3. Aturan penggunaan column yang ketat. Semua penulisan isi program harus diawali pada kolom ke 7.  Setahu saya, ini berguna bagi penyimpanan data untuk Punch Card (“Kartu Plong”) yang menjadi cikal bakal disket pipih nan lebar, sebagai media penyimpanan pada jaman doeloe.
4. Pangkat ditulis dalam ** dan bukan ^
5. dan kekunoan-kekunoan lainnya

Semua kekunoan di atas membuat banyak mahasiswa berpikir : Mengapa harus FORTRAN ? Iseng-iseng, hal ini saya tanyakan kepada Pak Wir, sebutan akrab Dr. Ir. Wiryanto Dewobroto, dosen teknik sipil Universitas Pelita Harapan, dan juga alumni Teknik Sipil UGM. Beliau mengasuh sebuah blog yang amat ramai di sini. Saya mengajukan pertanyaan sebagai berikut :

Pak Wir, iseng saya ingin bertanya, mengapa di kurikulum teknik sipil kami (UGM), diharuskan mempelajari bahasa FORTRAN? Bukankah beralih ke C++ atau Java akan lebih baik mengingat sebagian besar program dasarnya adalah dua bahasa itu? Sampai sekarang saya masih belum bisa menemukankelebihan FORTRAN dibanding C++

Terima kasih ^^

Keesokan harinya, ternyata pertanyaan saya langsung dijawab :

@yhouga
Lihat penjelasan atau pernyataan dari sdr Jedliem di atas. Saya kira saya sependapat dengannya. Karena bagi mahasiswa teknik sipil, penguasaan bahasa pemrograman hanya untuk melatih logika. Jadi gunakan saja bahasa pemrograman yang paling dikuasai dosennya. Kami di UPH tidak memakai Fortran tetapi memakai Visual Basic, alasan utamanya karena saya sebagai dosennya menikmati memakai bahasa program tersebut. Nggak usah mikir lagi.

Kalau ternyata nanti saya juga menguasai bahasa lain, misalnya C++ atau Java, maka ada kemungkinan saya pakai juga untuk mengajar. Pokoknya nggak ada yang kaku dengan bahasa program tersebut.

O ya, selain pemikiran tersebut ada baiknya suatu saat nanti dikaitkan dengan lisensi. Apakah kita sudah memakai program berlisensi. Jadi ada kemungkinan pakai yang versi “open”.

Sebuah pemikiran yang bagus. Jadi, alih-alih harus belajar dua kali jika ingin mempelajari Visual Basic (dengan graphical user interface yang sangat friendly), mengapa tidak langsung saja diajarkan Visual Basic agar mahasiswa dapat membuat program yang lebih kreatif, dengan tampilan yang menarik pula? Namun, jika memang ingin melatih logika semata, boleh lah dengan FORTRAN. Hanya, meminjam bahasa ibu, itu hal yang mindhogaweni menurut saya, membuat pekerjaan menjadi dobel. Namun, pemikiran untuk menggunakan bahasa yang “open” juga perlu dipertimbangkan. Bukankah saat ini dicanangkan program UGOS (UGM Goes Open Source) ?

Semoga bermanfaat!

-yhougam-

Genap dua kali sudah, saya sedikit berkecimpung dalam membantu mata kuliah Bahasa Pemrograman di Teknik Sipil UGM. Dua periode itu, membuat saya sedikit banyak mengetahui permasalahan yang kerap dihadapi mahasiswa baru di semester I, yang begitu beberapa bulan masuk kuliah langsung dijejali dengan logika-logika pemrograman. Bagi mahasiswa yang telah mendapatkan pelajaran TIK (Teknologi Informasi dan Komunikasi) semasa SMA, hampir dipastikan akan segera menyerap logika-logika tersebut dengan lancar. Termasuk saya yang ketika SMA, sekolah kami mempekerjakan mahasiswa-mahasiswa Teknik Elektro Universitas Brawijaya sebagai pengajar di pelajaran TIK.

Sedikit nostalgia, waktu itu ada 5 sekawan : Pak Ryan, Pak Bison, Pak Galih, Pak Fariz, dan juga Pak Aswin. Sengaja saya panggil Pak, karena saya menghormati mereka, walau keseharian kami biasa memanggil beliau-beliau dengan sebutan Mas saja. Karena yang mengajar mahasiswa, kami langsung saja diberi berbagai materi seputar jaringan, hardware (ketika SMP bahkan ada ujian merakit PC), dan tentu saja pemrograman berbasis C++. Bahkan, waktu itu kurikulum TIK kami disebut-sebut adalah yang paling maju di kota Malang. Dengan dasar seperti itulah, tidak sulit bagi saya untuk memahami 2 SKS mata kuliah Bahasa Pemrograman di Teknik Sipil UGM.

Sementara bagi mereka yang kurang beruntung, dan tidak mendapatkan dasar-dasar materi Pemrograman semasa SMA (termasuk yang mengherankan adalah terkadang dari SMA-SMA favorit di Jogja), harus berjuang sedikit lebih keras untuk bisa memahami logika-logika percabangan, perulangan (looping), algoritma, flowchart, dan tetek bengek pemrograman lainnya.

Bagi mereka yang kritis, biasanya ada dua pertanyaan seputar mata kuliah Bahasa Pemrograman yang kerap ditanyakan :

1. Mengapa di jurusan seperti Teknik Sipil harus belajar Bahasa Pemrograman ?
2. Mengapa harus mempelajari bahasa FORTRAN ?

Untuk pertanyaan yang pertama, biasanya akan segera terjawab dengan sendirinya jika telah memasuki perkuliahan di semester-semester 4 dan 5. Terutama ketika berhadapan dengan kasus-kasus Hidrolika, peninggian muka air yang harus diselesaikan dengan selembar penuh hitungan iterasi, penyelesaian debit perpipaan dengan sistem looping, perancangan dimensi saluran drainase perkotaan dengan software DUFLOW, dan juga yang paling legendaris dalam hal struktur, penyelesaian dengan Finite Element Method dan juga penyelesaian dengan Metode Matriks. Ruwetnya hitungan-hitungan tersebut akan dapat diselesaikan dengan program sederhana, yang dapat dibuat di kalkulator programmable sekalipun. Belum jika kita menemukan kasus-kasus lokal dalam proyek, yang memaksa kita untuk membuang software yang telah ada, dan membuat sendiri software yang sesuai dengan kasus tersebut.

Sementara untuk pertanyaan yang kedua, mengapa harus FORTRAN ? Perlu diketahui, biasanya dalam tengah semester awal, mata kuliah Bahasa Pemrograman murni mengajarkan program FORTRAN. Sebuah bahasa yang boleh dibilang kuno, yang dibuat di periode 1960-an saat awal munculnya komputer IBM. Jenis FORTRAN yang digunakan adalah WATFOR77, versi FORTRAN yang sesuai dengan namanya, dirilis pada tahun 1977. Lebih dari 30 tahun lalu.

Bahasa FORTRAN, ketika pertama kali mempelajarinya, terlihat sebagai bahasa yang lebih dulu diciptakan daripada C++, Pascal, dan bahasa pemrograman lainnya. Hal ini terasa dari beberapa poin :

1. Tidak adanya operator “>” (lebih dari), “<” (kurang dari) dan “=” (sama dengan) pada perintah IF. Sebagai gantinya, digunakan LT (Less Than “<”), GT (Bigger Than “>”), dan EQ (Equal “=”)
2. Sangat Case Sensitive, terutama ketika kita menyimpan file, harus dengan huruf kapital, jika tidak, jangan harap bisa dibaca oleh Compiler.
3. Aturan penggunaan column yang ketat. Semua penulisan isi program harus diawali pada kolom ke 7.  Setahu saya, ini berguna bagi penyimpanan data untuk Punch Card (“Kartu Plong”) yang menjadi cikal bakal disket pipih nan lebar, sebagai media penyimpanan pada jaman doeloe.
4. Pangkat ditulis dalam ** dan bukan ^
5. dan kekunoan-kekunoan lainnya

Semua kekunoan di atas membuat banyak mahasiswa berpikir : Mengapa harus FORTRAN ? Iseng-iseng, hal ini saya tanyakan kepada Pak Wir, sebutan akrab Dr. Ir. Wiryanto Dewobroto, dosen teknik sipil Universitas Pelita Harapan, dan juga alumni Teknik Sipil UGM. Beliau mengasuh sebuah blog yang amat ramai di http://wiryanto.wordpress.com Saya mengajukan pertanyaan sebagai berikut :

Pak Wir, iseng saya ingin bertanya, mengapa di kurikulum teknik sipil kami (UGM), diharuskan mempelajari bahasa FORTRAN? Bukankah beralih ke C++ atau Java akan lebih baik mengingat sebagian besar program dasarnya adalah dua bahasa itu? Sampai sekarang saya masih belum bisa menemukankelebihan FORTRAN dibanding C++

Terima kasih ^^

Keesokan harinya, ternyata pertanyaan saya langsung dijawab :

@yhouga
Lihat penjelasan atau pernyataan dari sdr Jedliem di atas. Saya kira saya sependapat dengannya. Karena bagi mahasiswa teknik sipil, penguasaan bahasa pemrograman hanya untuk melatih logika. Jadi gunakan saja bahasa pemrograman yang paling dikuasai dosennya. Kami di UPH tidak memakai Fortran tetapi memakai Visual Basic, alasan utamanya karena saya sebagai dosennya menikmati memakai bahasa program tersebut. Nggak usah mikir lagi.

Kalau ternyata nanti saya juga menguasai bahasa lain, misalnya C++ atau Java, maka ada kemungkinan saya pakai juga untuk mengajar. Pokoknya nggak ada yang kaku dengan bahasa program tersebut.

O ya, selain pemikiran tersebut ada baiknya suatu saat nanti dikaitkan dengan lisensi. Apakah kita sudah memakai program berlisensi. Jadi ada kemungkinan pakai yang versi “open”.

Sebuah pemikiran yang bagus. Jadi, alih-alih harus belajar dua kali jika ingin mempelajari Visual Basic (dengan graphical user interface yang sangat friendly), mengapa tidak langsung saja diajarkan Visual Basic agar mahasiswa dapat membuat program yang lebih kreatif, dengan tampilan yang menarik pula? Namun, jika memang ingin melatih logika semata, boleh lah dengan FORTRAN. Hanya, meminjam bahasa ibu, itu hal yang mindhogaweni menurut saya, membuat pekerjaan menjadi dobel. Namun, pemikiran untuk menggunakan bahasa yang “open” juga perlu dipertimbangkan. Bukankah saat ini dicanangkan program UGOS (UGM Goes Open Source) ?

Semoga bermanfaat!

As the code becomes parallel, the complexity increases. One of the common problem in parallel code is data race. To debug this problem, I encounter that Intel’s Thread Checker to be a good software. You can download from the Intel’s website for evaluation. I found about this when I read this article

When I run tcheck_cl, I realized that some parts has a data race. So I need to eliminate this. I couldn’t make it private since it will give segmentation fault (refer to this bug). So I simply add the !$omp critical to the section that has the data race. And it solves the problem.

Anyway, Intel’s Thread Checker really helped to identify which section of the code has the problem.

Toda Sexta postarei um artigo mais descontraído.

Como os programadores matam os dragões…

Java
Chega, encontra o dragão. Desenvolve um framework para aniquilamento de dragões em múltiplas camadas.
Escreve vários artigos sobre o framework, mas não mata o dragão.

.NET
Chega, olha a idéia do Javanês e a copia, tenta matar o dragão, mas é comido pelo réptil.

ASP
Os componentes necessários para levantar a espada são proprietários e caros. Outros tantos componentes proprietários para achar a localização do dragão, e mais outros tantos a localização da princesa. Chama então seu amigo programador de PHP.

C
Chega, olha para o dragão com olhar de desprezo, puxa seu canivete, degola o dragão. Encontra a princesa, mas a ignora para ver os últimos checkins no cvs do
kernel do linux.

C++
Cria um canivete básico e vai juntando funcionalidades até ter uma espada complexa que apenas ele consegue entender … Mata o dragão, mas trava no meio da ponte por causa dos memory leaks.

COBOL
Chega, olha o dragão, pensa que tá velho demais para conseguir matar um bicho daquele tamanho e pegar a princesa e, então, vai embora de volta ao seu
mundinho.

Pascal
Se prepara durante 10 anos para criar um sistema de aniquilamento de dragão… Chegando lá descobre que o programa só aceita lagartixas como entrada.

VB
Monta uma arma de destruição de dragões a partir de vários componentes, parte pro pau pra cima do dragão e, na hora H, descobre que a espada só funciona
durante noites chuvosas…

PL/SQL
Coleta dados de outros matadores de dragão, cria tabelas com N relacionamentos de complexidade ternária, dados em 3 dimensões, OLAP, demora 15 anos para processar a informação. Enquanto isso a princesa virou lésbica.

PHP
Pesquisa bancos de scripts e acha as classes de construção de espada, manuseio da espada, localização da princesa e dragão. Remenda tudo e coloca umas firúlas próprias.
Mata o dragão e casa com a princesa. Como tudo foi feito com gambiarras, o dragão um dia vai ressuscitar e comer os dois.

Ruby
Chega com uma baita fama, falando que é o melhor faz tudo, quando vai enfrentar o dragão mostra um videozinho dele matando um dragão … O dragão come ele de tédio.

Smalltalk
Chega, analisa o dragão e a princesa, vira as costas e vai embora, pois eles são muito inferiores.

ASSEMBLY
Acha que está fazendo o mais certo e enxuto, porém troca um A por D, mata a princesa e transa com o dragão.

Shell
Cria uma arma poderosa para matar os dragões, mas na hora H, não se lembra como usá-la.

Shell (2)
O cara chega no dragão com um script de 2 linhas que mata, corta, stripa, pica em pedacinhos e empalha o bicho, mas na hora que ele roda, o script aumenta,
engorda, enfurece e coloca álcool no fogo do dragão.

Fortran
Chega desenvolve uma solução com 45000 linhas de código, mata o dragão vai ao encontro da princesa …
mas esta o chama de tiuzinho e sai correndo atrás do programador java que era elegante e ficou rico.

FOX PRO
Desenvolve um sistema para matar o dragão, por fora é bonitinho e funciona, mas por dentro está tudo remendado. Quando ele vai executar o aniquilador de
dragões lembra que esqueceu de indexar os DBF’s.

CLIPPER
Monta uma rotina que carrega um array de codeblocks para insultar o dragão, cantar a princesa, carregar a espada para memória, moer o dragão, limpar a sujeira, lascar leite condensado com morangos na princesa gostosa, transar com a princesa, tomar banho, ligar o carro, colocar gasolina e voltar pra casa. Na hora de
rodar recebe um “Bound Error: Array Access” e o dragão come ele com farinha.

ANALISTA DE PROCESSOS
Chega ao dragão com duas toneladas de documentação desenvolvida sobre o processo de se matar um dragão genérico, desenvolve um fluxograma super complexo para libertar a princesa e se casar com ela, convence o dragão que aquilo vai ser bom pra ele e que não será doloroso. Ao executar o processo ele estima o esforço e o tamanho do estrago que isso vai causar, consegue o aval do papa, do Buda e do Raul Seixas para o plano, e então compra 2 bombas nucleares, 45 canhões, 1 porta aviões, contrata 300 homens armados até os dentes, quando
na verdade necessitaria apenas da espada que estava na sua mão o tempo todo.

I write a fortran code and save the data in VTK legacy format with structured points. However, for certain data somehow Paraview and Mayavi are unable to read the VTK file. It gives a “Generic Warning”, and says that
Error reading ascii data!

After days of debugging, it turns out that the cause is my precision, I used double precision for my data, and there are some data with 3 digit for the exponents, and I write to the file using ES format. This turns out cannot be read by Paraview and Mayavi. So when I changed the format to F which is just a floating point with decimal points (no exponents), then paraview and mayavi2 can read my data. But this is of course is not a good format to store. And after looking for more detail in Fortran languange, I can set my digit in the exponent to only 2. This is the format in fortran

ESw.dEe

where “e” is the digit for the exponent.

The other thing is that In VTK, I set the data type as “float”, I need to change this to “double”.

The days of C++ as a general purpose programming language are quickly ending for most developers. There are still lots of great uses for C++, particularly for OS-level work, low-level work (embedded devices, device drivers, etc.), certain high-performance applications, and applications where the overhead of a system like .NET or Java would be too heavy (like an office suite). Some developers will continue to use C++ for applications that other, less complex languages can handle as well. But for the typical developer, C++ is a big headache for minimal gain.

The performance issues that most developers face are not the kinds of issues that moving to native code will resolve; once you take performance out of the equation, C++ is a fairly unattractive option for application development in most cases.

The handful of C++ developers that I’ve talked to say using C++ in the .NET managed environment is not particularly attractive to them; this takes away much of the opportunity to use it in a Web development capacity, unless you want to use it in the CGI model. There are good things about CGI (less overhead, simple conceptual model) and bad things about CGI (your application has to be “aware” of many more low-level tasks). From what I’ve heard, under the .NET CLR, C++ loses its speed, as well as many of the things that make C++ useful.

This is not to say that C++ is going away any time soon. I see C++ joining the ranks of COBOL and FORTRAN as a legacy language with a massive installation base and a need for people to maintain/extend existing applications for more than 50 years. In addition, a number of new development projects will be started in it for a variety of reasons (familiarity, library support, tradition/habit, cultural, etc.). I also suspect that it will pick up a reputation as a “dead” language (again, like COBOL), due more to a lack of buzz and hype than actual non-usage (also like COBOL).

I don’t want to make the future for C++ sound dismal; if anything, I think there is great potential for C++ developers to do quite nicely for themselves. If you’re a C++ developer, I suggest that you stick with the language. Are the things you’re working on flashy or get the same attention as Web applications in the mainstream publications? No. But with the current salary structures, I feel that experienced C++ developers will see very nice paychecks for some time. In addition, as the remaining C++ work is of higher difficulty and fewer people learn C++ (it isn’t taught as frequently in colleges these days), I expect C++ developers to have more job security and better compensation than .NET or Java developers over the long run.

While C++ in Web development is not likely to become mainstream any time soon and desktop application development in C++ becomes less common, I think there is a lot of upside opportunity for C++ in certain aspects of cloud computing. For some projects (think of ones that are well suited to supercomputers), the cloud offers C++ developers a way to get the same benefits of grid computing but with much more flexibility. There is a lot of overlap between those projects and the kinds of projects for which developers regularly use C++. As a result, I think cloud computing will replace or supplement grids and supercomputers in many projects and will provide an excellent opportunity to see C++ used in new and innovative ways.

Keep in mind that much of this analysis is focused on the Windows world. From what I see, the *Nix development community is still very C/C++ oriented. C++ developers who are concerned about dwindling opportunities in Windows should definitely take a look at *Nix development.

I believe that C++ will slowly fade into the background, but it will neither die nor will it ever become unimportant. While most developers I know have never touched C++ in a real-world environment, many developers would benefit from learning it, if only to gain some appreciation of various languages, including Java, .NET, and Ruby.

Intel firması derleyici konusunda yaptığı ataklarda Fortran’ı da unutmadı. Bunun sebebi ise dünyada vazgeçilemeyecek kadar büyük maliyette Fortran kodunun bulunması. Zaten son dönemde eklenen özellikler ile dilin yapabileceklerinin de sınırı epeyi genişletildiğinden Fortran tahmin edildiği gibi ölmedi. Pek çok Fortran kullanıcısı gibi ben de onlarca programlama dili öğrenip çeşitli yazılımlar geliştirmiş olsam da iş mühendislik hesabına gelince Fortran’ı -alışkanlıklar nedeniyle C’yi göz ardı edersem- bu ayakta tek geçiyorum.

Yıllarca kendimi Windows’a hapsettikten sonra Vista macerasıyla zevki kırılan birisi olarak Linux bu konuda da bir yenilikler zinciri başlattı. Çünkü Intel ticari olmayan amaçlar için Fortran derleyicisinin son sürümünü ücretsiz dağıtıyor. Hatta Akademik programı var. Eğer akademisyenseniz 370 dolara öğrenciyseniz 130 dolara bu derleyicinin son sürümüne ekleriyle birlikte sahip oluyorsunuz. Koduna güvenen herkes için değecek bir fiyat. Bu ülke araştırma projelerinin bütçelerinde onbinlerce dolarlık yabancı yazılımların satın aldığını da görüyor olmasına rağmen 3-4 yüz dolara derleyici alıp kod geliştiren araştırmacıları da görüyor.

Sonra henüz bütün seti satın almadıysanız ücretsiz takılıyorsanız ve Linux’un engin denizinde yol almaya başlamışsanız sizi dünyanın gelmiş geçmiş en iyi kod geliştirme ortamı olan Eclipse ve ona Fortran için entegre olmuş Photran sizi ücretsiz olarak beklemekte.

Bütün bunlar size çok ciddi bir özgürlük kazandırıyor. Kendinizi sonuna kadar rahat hissediyorsunuz. Bir gün geliştirdiğiniz kod paraya döner hale gelirse de Intel’e vereceğiniz para hazır kodların ülkeye maliyeti yanında solda sıfır kalır.

Bu sebeple herkese Linx altında Fortran derleyicisini Photran ile entegre ederek kullanmaları şiddetle tavsiye olunur.

Selam, sevgi, saygı ve dostluk…