F90 – Resumão – IF – UFRGS - Conteúdo de Física na Net Portal Física

Resumo adaptado do original: (IF-UFRGS – Origem: Wikipédia, a enciclopédia livre.).

A linguagem de programação Fortran foi desenvolvida na década de 1950 e continua a ser usada hoje em dia. O nome tem como origens a expressão Formula Translator (ou Translation). A linguagem Fortran é principalmente usada em ciências da computação e análise numérica.

O primeiro compilador de Fortran foi desenvolvido para o IBM 704 em 1954-57. O compilador era optimizado, pois os autores acreditavam que ninguém iria usar essa linguagem se a sua prestação não fosse comparável com a da linguagem assembly A linguagem Fortran foi largamente adoptada por cientistas para a escrita de programas numericamente intensivos, o que encorajou os produtores de compiladores a escrever compiladores que gerassem código mais rápido. A inclusão de um tipo de dados de número complexo na linguagem tornou a linguagem Fortran particularmente apta para a computação científica. Ao longo dos tempos foram-se dando algumas revisões da linguagem. Entre elas encontram-se:

FORTRAN IV (também conhecida como FORTRAN 66)
FORTRAN 77,
Fortran 90, 95, 2003 e 2008

Olá Mundo

Exemplos de código típico e mínimo para testar se o Fortran 90 está instalado (há um equivalente para cada linguagem)

 print*,"olá Mundo" 
 end

Editar

Essas duas linhas devem estar dentro de um arquivo (ex: ola.f90), isto pode ser feito, no terminal linux por exemplo, com cat linha por linha ou com um editor de texto qualquer, por exemplo no linux (vi, nano, emacs, etc…)

Compilar

ola.f90 deve ser compilado, isto é traduzido de linguagem de Fortran (código fonte) para Assembler (código objeto ou de maquina) com a seguinte instrução:

> f90 ola.f90 -o ola

f90 invoca o compilador que faz a dita conversão

-o nome é uma opção para nomear o código objeto (o programa executável) que no exemplo terá o nome ola

sem essa opção o executável criado terá o nome a.out

Rodar (executar)

ola.f90 não é executável, o computador não sabe o que fazer com ele, para isso o compilador traduz para una série de instruções que o computador sim sabe executar)

Para rodar fazemos:

> ./ola

./ é uma espécie de “comando” que deve ir na frente do nome do programa, necessário para que o intérprete de comandos (shell) do terminal (ou sistema) linux saiba que o programa está no diretório de trabalho. Ao ser executado será exibido na tela do computador o seguinte:

olá Mundo

Desta forma completamos o pequeno ciclo de teste de existência de compilador Fortran 90 no computador que estamos usando e também de ter criado nosso primeiro código Fortran, compilado e executado o mesmo.

Código fonte do programa Olá Mundo na versão padrão de formato-sintaxe: program ... end program:

program ola
print*,"olá Mundo"
end program ola

Esta versão está no estilo Fortran 90, mesmo que com declarações opcionais como program nome (terminada com end program) é útil quando temos códigos longos com subrotinas (programas chamados pelo programa principal)

NOTA: A título de curiosidade, a única declaração obrigatória em Fortran é o comando end.

Variáveis

Em linguagem de programação Variáveis são nomes definidos pelo programador onde podem ser armazenados números ou texto para uso posterior, sobre os quais podem ser feitas operações matemáticas, lógicas ou de manipulação de caracteres, e que também podem ser alteradas ou substituídas por novos conteúdos dependendo das instruções do programa.

Um programa é feito de constantes, variáveis e instruções de operações entre elas que definem uma tarefa a ser desenvolvida automaticamente pelo computador.

Em Fortran existem cinco tipos primitivos de variáveis:

Tipo Primitivo	Descrição
INTEGER	inteira, exemplo: 0, 1, 100 -34
REAL	real ou de ponto flutuante, exemplo: 3.14, -0.003, 3.2E-10, -5.02E100 (notação científica)
COMPLEX	complexa, exemplo: \[(2,3 + 0i):\] (2.3, 0) \[(-i):\] (0, -1.0) \[(1,2 \times 10^{-5} – 1 \times 10^{-3}i):\] (1.2E-5, -0.001)
LOGICAL	logica, exemplo: .TRUE. .FALSE.
CHARACTER	caractere, exemplo: “ola”, “Entre seu nome!”

No programa a definição de variáveis, ou declaração de variáveis, é feita no início e depois seguem as instruções de execução, veja a forma no seguinte exemplo:

 Program test
 Integer I, J, K
 Real  a, x, velocidade, tempo_0, Temp
 Logical prova
 Character*10 Titulo, nome
 ...

Bytes

Em computação byte é uma unidade de medida de armazenagem de informação, tipicamente de oito bits.

O bit, pela sua vez, é a mínima quantidade de informação possível, ou seja a menor unidade de medida de armazenagem de informação.

É onde podemos guardar apenas dois estados: 0 ou 1, sim ou não, etc.

Menos do que isso não é informação nenhuma, pois de fato, menos do que isso é um estado único que só pode ser sempre 0 ou 1 ou 2, mais é sempre igual, ele só pode tomar esse único valor e não mudar. Para ter informação então, no mínimo precisamos de dois estados: ora é um, ora é zero.

Por isso alguém chamou o bit de

a mínima diferença que faz a diferença.

A origem do bit é anterior ao computador, provavelmente de 1936, e o nome vem da contração de binary digit em inglês.

Os dígitos binários eram usados anteriormente e são apenas uma representação dos números como é a decimal, porem com base nas potências de 2. Assim, um número binário só pode estar composto de 0 e 1 (como um decimal esta só composto dos dígitos de 0 a 9).

Como exemplo, o numero binário (7 bits) 1001100 representa o número decimal 76.

1	0	0	1	1	0	0
6	5	4	3	2	1	0
\[1 \times 2^6\]	\[0 \times 2^5\]	\[0 \times 2^4\]	\[1 \times 2^3\]	\[1 \times 2^2\]	\[0 \times 2^1\]	\[0 \times 2^0\]
64	0	0	8	4	0	0

\[64+8+4 = 76\]

(Mais informações: Byte, Bit, Sistema binário )

Entrada e Saída

São as intruções READ e WRITE (ou PRINT)

O READ permite entrar dados antes ou durante a execução do programa.

Com o WRITE podemos ver os resultado final da execução do código ou

acompanhar o desenvolvimento do mesmo.

Exemplos:

 READ (*,*) I, Velocidade
 ...
 Write (*,*) Titulo, Temp

Esta última é equivalente a

 PRINT*, Titulo, Temp

Todas essa representam entrada/saída pelo terminal. Entrada pelo Teclado. Saída na Tela.

São as chamadas entrada/saida padrão (standard)

Se quisermos ler ou gravar num arquivo:

OPEN (1, FILE="entrada.dat")
READ(1,*) I, Velocidade

OPEN (2, FILE="saida.dat")
WRITE (2,*) Titulo, Temp

 CLOSE(1)
 CLOSE(2)

Operadores matemáticos

Operador	Função
`.`	ponto decimal
`+`	soma
`-`	subtração
`*`	multiplicação
`/`	divisão
`**`	potência
`( )`	parênteses

Exemplo usando todas as operações

Equação:

\[z = 3,14 + \frac{\sqrt(x)}{3} – 4,3\times10^{-6}\]

z = 3.14+x**(1./2)/3 - 4.3E-6

Funções intrínsecas

Estas são as funções pré-definidas no compilador, entre as quais encontramos a mais comummente utilizadas como

raiz quadrada, trigonométricas, exponencial, etc:

A sintaxe de uso é o nome da função (geralmente três letras) e o argumento entre parêntesis:

raiz quadrada: sqrt
trigonométricas: sin, cos, tan, asin
exponencial: exp
logaritmo natural: log
módulo: abs

Exemplos:

a = sqrt(2.)
b = sin(x)
Pi = 2.*asin(1.) ! uma forma pratica de ter Pi como constante dentro de um programa
y = (x - exp(-lam*x**2))/sqrt(x)
z = x**(1./4.) ! Cuidado Teste com 1/4 em lugar de 1./4. e verá a grande diferença

Lista das funções intrínsecas mais comuns

Function        Generic  Specific Data type
                name     name     Arg   Res

Square root     SQRT     SQRT     R     R
                         DSQRT    D     D
                         CSQRT    C     C

Exponential     EXP      EXP      R     R
                         DEXP     D     D
                         CEXP     C     C

Natural         LOG      ALOG     R     R
 logarithm               DLOG     D     D
                         CLOG     C     C

Common          LOG10    ALOG10   R     R
 logarithm               DLOG10   D     D

Sine            SIN      SIN      R     R
                         DSIN     D     D
                         CSIN     C     C

Cosine          COS      COS      R     R
                         DCOS     D     D
                         CCOS     C     C

Tangent         TAN      TAN      R     R
                         DTAN     D     D

Arcsine         ASIN     ASIN     R     R
                         DASIN    D     D

Arccosine       ACOS     ACOS     R     R
                         DCOS     D     D

Arctangent      ATAN     ATAN     R     R
                         DATAN    D     D
                ATAN2    ATAN2   2R     R
                         DATAN2  2D     D

Hyperbolic      SINH     SINH     R     R
 sine                    DSINH    D     D

Hyperbolic      COSH     COSH     R     R
 cosine                  DCOSH    D     D

Hyperbolic       TANH    TANH     R     R
 tangent                 DTANH    D     D

Ciclos ou laços (loop)

O grande lance dos programas (e do computador claro) é poder repetir operações sem ter que fazer uma de cada vez. Ou seja, utilizando a lógica para repetir tarefas que podem ter variações especificadas pelo programador dentro de um ciclo. Vejamos um exemplo usando a instrução do i =... end do:

...
do i = 1, 100
   x = i * 0.1;   y = x**2
   print*, x,y
end do
...

Nota

há diferenciação no shell onde maiúsculas e minúsculas são diferentes (A é diferente de a) no Fortran elas são interpretadas da mesma forma (print = PRINT = PrInT=…).

Porem, é recomendado, por uma questão de estilo, usar preferencialmente minúsculas, reservando maiúsculas para aquilo que se quer chamar a atenção. O Ponto e vírgula `;` permite separar instruções diferentes em uma mesma linha

a interpretação do do é a seguinte:

* comece com i=1
* faca o que esta entre `do ... end do`
* incremente i em 1 e veja se ultrapassou o valor máximo de i (no exemplo 100)
* se a resposta e não faca de novo  com o novo valor de i
* o ciclo se repete até o valor máximo de i (= 100 no exemplo)
* quando i ultrapassa esse valor o laço é interrompido e o programa continua com a seguinte instrução depois de `end do`

Pergunta: Qual será no exemplo acima o valor da variável i finalizado o laço?

Exemplo: MRUA

Neste exemplo criamos um programa que gera a trajetória x(t) (unidimensional) de um objeto com aceleração constante:

\[x(t) = x_0 + v_0t + \frac{at^2}{2}\]

 x(t) = x(0) + v(0)*t + 0.5*a*t²

program mrua
implicit none ! isto nos obriga a declarar todas as variáveis que usaremos (recomendável)
real :: x0, x, v0, a, t, dt  ! variáveis reais (ponto flutuante)
integer :: I, N

print*, 'entre x0, v0, a, t, N'
read*, x0, v0, a, t, N
dt = t/N  ! passo de tempo: definido pelo tempo máximo e o numero de pontos da trajetória
t = 0

do i = 1, N
   x = x0 + v*t + 0.5*a*t**2
   print*, t, x
   t = t + dt
end do

end program mrua

Desafio: este programa contém um erro que você deve consegue descobrir? Se sim comente no formulário no final do post. :o)

Controle de Fluxo (IF e CASE)

Possibilitam a escolha de rumos diferentes no programa dependendo de valores de variáveis ou condições entre elas.

IF … THEN …

A sintaxe é:

if (condição) then
     instruções executadas quando condição verdadeira
else
     instruções executadas quando condição falsa
end if

Para apreciar o poder dela primeiro vejamos:

Operadores relacionais e lógicos

Operador	análogo	Função
`==`	`.eq.`	igual
`/=`	`.ne.`	diferente
`>`	`.gt.`	maior
`>=`	`.ge.`	maior ou igual
`<`	`.lt.`	menor
`<=`	`.le.`	menor ou igual
`.and.`	`.and.`	E lógico
`.or.`	`.or.`	OU lógico
`.not.`	`.not.`	NÃO lógico

Eles são usados junto com as instruções de controlo como o IF (se), exemplo:

! exemplo de divisão a x b:
...
read*, a, b

if (b /= 0) then
    c = a/b
else
     stop 'b=0'
end if
print*, 'a/b=', c
...

CASE

Este serve para quando temos muitas opções ou escolhas para tomar rumos diferentes no programa. Vejamos a sintaxe diretamente no exemplo:

Select Case (J)
    Case(1)
      print*, "caso 1"
    Case(2)
      print*, "caso 2"
    Case (4:10)
      print*, "caso 4-10"
    ....
    Case Default
      print*, "não foi nenhum dos casos contemplados"   
End Select

Ciclo infinito (Do … end Do)

Se não especificamos com um índice o numero de vezes que o laço deve ser executado este será executado infinitamente. Exemplo:

do
   print*, "seja aproximadamente exato, não exatamente errado!"
enddo

o programa ficará rodando infinitamente, fazendo a frase ser escrita na tela em um laço infinito, que só poderá ser cancelado teclando CTRL C no terminal.

Para que um ciclo infinito tenha utilidade prática, em algum momento (lógico) do programa é preciso lançar mão do comando exit, como no exemplo abaixo:

 Do
     Print*, "Digite o número 1 ou 2"
     Read*, num
     Select Case (num)
     Case (1)
        Print*, "Você digitou 1"
        Exit
     Case (2)
        Print*, "Você digitou 2"
        Exit
     Case Default
        Print*, "Opção inválida"
     End Select
 End Do

Se queremos um programa a prova de fogo (que não de erro por entrada de dados) podemos fazer uma pequena modificação:

Na qual foi usado ERR=100 para o código se recuperar de dado inesperado, transferindo o controle do programa para a linha com esse número em caso de erro (texto em lugar de número por exemplo).

Ciclo condicional (Do While)

O do while é executado enquanto uma condição lógica é satisfeita. Exemplo:

  i = 1
  do while ( i .le. 100 )
    i = i + 1; print *, i
  end do

Resulta equivalente a:

  do i + 1, 100
     print*, i
  End do

Porem é de utilidade quando a condição não é trivial. Exemplo:

  fac = 1. + 1./3.;  imax = 1.E6
  I = 1;  print*, I
  do while ( i .lt. imax )
     i = fac * i;  print*, i 
  end do

Arranjos (Arrays)

Tecnicamente, são tipos de dados ou variáveis estruturadas.

De outra forma, um arranjo é um conjunto de dados escalares, todos do mesmo tipo (REAL ou INTEGER, etc), organizados de

forma regular de maneira que é fácil atribuir ou extrair valores deles.

Assim como um escalar real é definido como

 Real A

um array real pode ser definido assim:

 Real A(10)

No primeiro exemplo A pode conter apenas um valor por vez, enquanto no segundo, 10 diferentes valores podem ser armazenados

em A, cada um dos quais pode ser referido via o índice assim

 A(1) = 0;  A(2:7) = 4;  A(8) = -1;  A(9:10) = 2

  Do i =2, 10, 2
     print*, A(i) !imprimindo só as posições pares
  End do

Rodando esse programa teremos na tela:

Arranjos Estáticos

Arranjo estático e aquele cuja dimensão (números de elementos do arranjo) é definido nas declarações (instruções no início do programa e que não são executadas durante o tempo de rodada), é um número fixo que não pode ser mais alterado durante a execução. Exemplos:

 Real A(10), B(30)
 Integer GG(10,10), ax(5,20), FG(6,10,5)

Também pode ser assim:

 Character*10, Dimension (10,10) :: A, B, C, D

A dimensão é definida com Dimension(10,10) para todos os arrays (A,B,C,D), todos eles da mesma forma (shape).

NOTA: quando mais de um atributo é aplicado a uma ou mais variáveis (arrays ou escalares) os dois pares de pontos (:: quadrupleto) entre atributos e variáveis é mandatório sempre que houver atribuição inicial na declaração.

Arranjos Dinâmicos

A diferença com os anteriores é que a dimensão destes é definida durante durante a execução do programa, pudendo até ser modificada durante o mesmo. Isto se consegue com o atributo ALLOCATABLE. Vejamos o exemplo:

 Real, Allocatable :: A(:), B(:,:)
 Integer, Allocatable :: K(:), AJ(:,:,:)

Nesses exemplos definimos A e K como arrays vetoriais, B matricial, e AJ tensorial.

O número de elementos deles é definido durante o tempo de execução com ALLOCATE.

A seguir um exemplo que ilustra a vantagem de usar arranjos dinâmicos

Real, Allocatable :: A(:), B(:,:)
Integer I, N, M

Read*, N

Allocate (A(N), B(N,N))

Do I = 1, N
   A(I) = 2*I
End Do

Do i =1, N
   B(i,:) = A + i
End do

...
if (j .gt. 0) then
  Deallocate (B)
else 
  Deallocate (A)
end if

...
Select Case (J)
 Case(-1)
   Allocate (A(N/2)
 Case(2) 
   Allocate (B(N,N/2))
End Select case

Funções e Subrotinas

São módulos independentes para fazer tarefas especificas definidas pelo programador.

Com a declaração FUNCTION é possível definir funções matemáticas diferentes ou além do conjunto de funções intrínsecas do compilador. Elas podem ser de uma ou mais variáveis, porem devolvem um valor só.

ATENÇÃO: a sintaxe na declaração de função se escreve g(x) (primeira linha nos exemplos onde se declara o nome da função e das variáveis). Na definição da função, implementada em Fortran, se escreve:

 g = ...

não se colocam os argumentos ou variáveis a esquerda, elas devem aparecer só a direita da instrução.

Quanto a subrotina ela e um procedimento mais geral que pode operar em um conjunto ou todas as variáveis que estão no programa principal, e no retorno ao programa não há limite nas variáveis que pode devolver. Vejamos a sintaxe de cada uma:

Function

A sintaxe é:

 REAL FUNCTION NOME(VARIAVEL)
 ...
 NOME = operações com VARIAVEL
 END FUNCTION NOME

Vamos para o exemplo:

\[g(x)=x^{2} sin(x)\]

Real Function g(x)
Implicit None
Real x
 g = x**2 * sin(x)
End Function g

Exemplo de duas variáveis:

Real Function paraboloide(z,w)
Implicit None
Real z,w

  paraboloide = z**2 + w**2

End Function paraboloide

É possível também definir funções não analíticas, como por exemplo:

...
if (x .le. 0) then
    g = 0.
else
    g = 1.
end if
...

Logo na execução, dentro do programa principal, se chama como qualquer função intrínseca do Fortran. Vejamos um exemplo completo de programa com definição de função.

! este e o program principal que usa a função h(x)
program exemplo_uso_funcao
implicit None
real, external :: h(x) ! external pq está depois do end program
real :: x, dx, y
integer :: i,N

Print*, "programa para construir tabela de h(x)"
Print*, "entre o primeiro e ultimo x e o numero de pontos a graficar"
read*, x0, xf, N

dx = (xf - x0)/N

do i = 0, N
   x = x0 + i*dx
   print*, x, h(x)
enddo

end program exemplo_uso_funcao

! aqui a definicao de funcao h(x)
real function h(x)
real x
   h = exp(-x**2)/2.0
end function h

Subrotina (subroutine)

A sintaxe de subrotina e:

 SUBROUTINE EXEMPLO (Argumentos)
 ...
 END SUBROUTINE

No meio todas as instruções Fortran são validas, pois a subrotina é um programa. A idéia é que se um conjunto de operações se repete muito, então resulta conveniente colocá-las num módulo e chamá-lo cada vez que for necessário, assim:

PROGRAM Principal
...
  CALL EXEMPLO(Argumentos)
  ...
  DO
     ...
     CALL EXEMPLO(Argumentos)
     ...  
  END DO
 ...
END PROGRAM Principal

Exemplo:

program main
implicit None
real, dimension(3,3) :: A=0., B=0.

  call matriz(A,B)
  print*,"Matriz A"
  print'(3(f5.3,1x))',A
  
  print*,"Matriz B"
  print 100,B

100 format (3(f5.3,1x))

stop("Programa finalizado!")
end program main

subroutine matriz(A,B)
real, dimension(3,3) :: A,B

    A = B + 1

end subroutine matriz

Resultado:

Matriz A                                                  
1.000 1.000 1.000
1.000 1.000 1.000
1.000 1.000 1.000   
Matriz B                   
0.000 0.000 0.000
0.000 0.000 0.000
0.000 0.000 0.000

Module

Esta estrutura permite definir variáveis de uso comum entre o programa principal e as subrotinas sem se preocupar de passar elas como argumentos das mesmas. Vejamos um exemplo:

module variaveis
  real a, b
  integer N
end module

Depois no programa e subrotina

program principal
use variaveis           !o main sabe agora de a,b,N
implicit None
...
  read (*,*) a, b, N    !estes valores de a,b,N serão passados para sub
  call Sub
...
end program

subroutine sub
use variaveis           !a Subrotina sabe de a,b,N
implicit none
...
  h = (a-b)/N           !e pode usá-las com os valores que tinham quando foi chamada
...
end subroutine

Links

—

IDPGPT-1890

F90 – Resumão – IF – UFRGS