%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  MAP05 - Metodos Matriciais Computacionais - PPGMAp UFRGS
%%  Mat01050 - Algebra Matricial Computacional - DMPA UFRGS

Tecnicamente, esse nao eh um script, pois nao pode ser executado.
Entenda como um conjunto de instrucoes que facilmente podem
ser carregados com o mouse.

  Estudo de caso: usando gcc para chamar uma rotina CBLAS na bulldog
(CBLAS eh uma versao de BLAS que foi convertida com f2c)

  whereis libcblas    ( /usr/lib/libcblas.a  )

edite arquivo dnrm2_cd.c:   (nano dnrm2_cd.c)

#include <stdio.h>
#include <cblas.h>
main() {
  // serah um driver para rotina (funcao) dnrm2 da biblioteca CBLAS
  double x[2000],norma;
  float xi; 
  int i,n;
  printf("Qual o valor de n, o numero de componentes ? ");
  scanf("%i",&n);
  printf("Qual o valor das componentes, linha a linha ?\n");
  for (i=0;i < n;i++){ scanf("%f",&xi); x[i]=xi;}
  norma = cblas_dnrm2(n,x,1); 
  // normalmente, usamos cblas_ antes no nome da respectiva rotina BLAS
  printf("Valor da norma, via (cblas_dnrm2) : %f \n",norma);
}

no servidor bulldog, compile com  gcc via
      gcc dnrm2_cd.c -lcblas -o dnrm2_cd
execute com ./dnrm2_cd
%%
%% Feito    JBC 10Set2012
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%