Tópico: Corujão do Clube Cético

[Andre]

Desde que vim pra Blumenau, só acredito que está calor mesmo quando alguém daqui ou de Joinville fala. Pessoal de SP fica reclamando, mas quando fui pra lá, achei um alívio (e fiquei um bom tempo).

[SnowRaptor]

E seu parâmetro de frio? Quanto dá 25⁰C tem gente em Salvador com tr6es blusas, tremendo e xingando.

[Andre]

Ah, frio eu gosto, então não fico reclamando.

Aqui tem madrugada que faz 6°C. Geralmente de manhã no inverno faz uns 12°C. Eu acho bom assim.

[SnowRaptor]

[Nina]

Hum... friozinho chegando, hora de dormir!

[Diegojaf]

Calor do inferno... todo dia. E esse mês de Fevereiro será ainda mais quente segundo disseram.

[Andre]

Hoje os termômetros do centro daqui estavam marcando 45°C às 14h50min, às 16h o google dizia que estava 38°C. Uma amiga saiu agora pouco (às 21h30min) e disse que os termômetros estão marcando 32°C. Durante o dia, a RBS (filial da Globo) disse que a sensação térmica estava de 56°C.

Vou morrer assim...

[SnowRaptor]

Vem pra SP. Toda tarde tem uma chuvinha pra refrescar.

[Fabi]

Ah, frio eu gosto, então não fico reclamando.

Aqui tem madrugada que faz 6°C. Geralmente de manhã no inverno faz uns 12°C. Eu acho bom assim.

Eu só reclamo quando faz -2°C ou se os meus dedos ficam congelados.

Eu já tive em blumenau e joinville. Aqui tá parecendo joinville... tá muito calor.

Não sei o que é pior, o meu ventilador ter queimado ou não ter ventilador pra comprar. 

[Barata Tenno]

Vem pra SP. Toda tarde tem uma chuvinha pra refrescar.

Dependendo de onde morar depois tem um desbarrancamento pra animar e a enchente pra exercitar o nado....

[Unknown]

Vem pra SP. Toda tarde tem uma chuvinha pra refrescar.

Depende do lugar. Terça-feira rolou um toró monstruoso na Vital Brasil e na minha casa não caiu uma gota de chuva. Nos últimos dias o tempo sempre fecha perto de casa, mas chover que é bom...

[SnowRaptor]

Já na USP, todo dia tem pé d'água.

E Joinville é insuportavelmente quente!

[SnowRaptor]

Passei a noite em claro.

Enquanto meu computador (quase) incansavelmente mastigava números do meu doutorado, eu tentei ir dormir e quando encostei a cabeça no meu travesseiro, comecei a bolar algoritmos pra pegar os dados que importam do resultado das contas que tavam rodando.

Provavelmente seria mais fácil e rápido fazer em perl, mas como dizem que a melhor linguagem é a que a gente conhece,  aí está:

Código: [Selecionar]

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>

void conta_atomos(FILE* logfile, int* nbasis, int* base, int* elemento);
void le_log(FILE* logfile, int* base, int *elemento, double** autovetores, double* autovalores, int nbasis);
void compara_pi(int* base, int* elemento, double** autovetores, double* autovalores, int nbasis, double threshold);
void imprime_matriz(int* base, int *elemento, double** autovetores, double* autovalores, int nbasis);
void libera(FILE* logfile, double** autovetores, double* autovalores, int nbasis);
void *mallocX (unsigned int nbytes);

int main(int argc, char** argv){
FILE *logfile;
double **autovetores;
double *autovalores;
double threshold;
int base[2000], elemento[2000];
int nbasis, i;

if (argc != 3 ){
printf("Erro! Numero invalido de argumentos!\nUso: %s logfile threshold\n", argv[0]);
return EXIT_FAILURE;
}



/* Abertura do arquivo */
printf("Abrindo. . .\n");

logfile = fopen(argv[1], "r");
if (logfile == NULL){
printf("Erro! Nao consegui abrir arquivo %s para leitura.\n", argv[1]);
libera(logfile, autovetores, autovalores, nbasis);
return EXIT_FAILURE;
}

printf("Contando. . .\n");
conta_atomos(logfile, &nbasis, base, elemento);

/* Alocacao */

printf("Alocando. . . nbasis=%d\n", nbasis);
autovetores = mallocX (nbasis * sizeof (double *));
for (i = 0; i < nbasis; i++)
autovetores[i] = mallocX (nbasis * sizeof (double));

autovalores = mallocX (nbasis * sizeof (double *));


/* Lendo o log e gravando matrizes e vetores */

le_log(logfile, base, elemento, autovetores, autovalores, nbasis);

/* Agora vamos imprimir!
  Debug de printf FTW! \o/
imprime_matriz(base, elemento, autovetores, autovalores, nbasis);
*/


threshold = atof(argv[2]);
compara_pi(base, elemento, autovetores, autovalores, nbasis, threshold);

/* FREE!! */
libera(logfile, autovetores, autovalores, nbasis);

return 0;

}

void conta_atomos(FILE* logfile, int* nbasis, int* base, int* elemento){
char string[80];
int i, bases=0;

/* Pula o arquivo até achar os orbitais */
do
fgets(string,80,logfile);
while (strncmp(string," Charge =  0 Multiplicity = 1",29) && !feof(logfile));

for(fgets(string,80,logfile);strncmp(string," Stoichiometry",14) && !feof(logfile); fgets(string,80,logfile)){
if(!strncmp(string," H", 2)){
elemento[bases]=1;
base[bases]=0; /* 0 == S */
bases++;
}
else if(!strncmp(string," C", 2)){
for (i=0; i < 4; i++){
elemento[bases+i]=6;
base[bases+i]=i; /* 0 == S; 1 == PX; 2 ==PY; 3 == PZ */
}
bases+=4;
}
}

*nbasis = bases;

}

void le_log(FILE* logfile, int* base, int* elemento, double** autovetores, double* autovalores, int nbasis){
char string[80];
int resto, nbasis_5, orbital, componente, i;
printf("Lendo. . . \n");
/* Pula o arquivo até achar os orbitais */
do
fgets(string,80,logfile);
while (strncmp(string,"     Molecular Orbital Coefficients",35) && !feof(logfile));

/* Comecando a ler orbitais. */

resto = nbasis % 5;
nbasis_5 = nbasis - resto;
for(orbital=0; orbital < nbasis_5; orbital +=5){
/* Pula duas linhas:
  1         2    ...
(A)--O    (B)--O  ...
*/

fgets(string,80,logfile);
fgets(string,80,logfile);

/* Primeira linha: autovalores: Joga fora primeiras 20 colunas: */
fgets(string,20,logfile);
fscanf (logfile, "%lf %lf %lf %lf %lf\n",
&autovalores[orbital],
&autovalores[orbital+1],
&autovalores[orbital+2],
&autovalores[orbital+3],
&autovalores[orbital+4]);

/* Guardando as componentes dos autovetores */
for (componente=0; componente < nbasis; componente++){
/* Primeiro, o base
  Eu sei que é redundante pôr isso no loop dos orbitais,
  mas isso implica em menos IFs.
*/
fgets(string, 20, logfile);

/* Agora, as componentes em si */
fscanf (logfile, "%lf %lf %lf %lf %lf\n",
&autovetores[orbital][componente],
&autovetores[orbital+1][componente],
&autovetores[orbital+2][componente],
&autovetores[orbital+3][componente],
&autovetores[orbital+4][componente]);
}
}

/* As últimas colunas tem que ser diferentes. Damn you, FORTRAN! */
/* Pula duas linhas:
  1         2    ...
(A)--O    (B)--O  ...
*/

fgets(string,80,logfile);
fgets(string,80,logfile);

/* Primeira linha: autovalores: Joga fora primeiras 20 colunas: */
fgets(string,20,logfile);

for (i=0; i < resto; i++)
fscanf (logfile, "%lf", &autovalores[orbital+i]);


/* Guardando as componentes dos autovetores */
for (componente=0; componente < nbasis; componente++){
/* Primeiro, o base
  Eu sei que é redundante pôr isso no loop dos orbitais,
  mas isso implica em menos IFs.
*/
fgets(string, 20, logfile);

/* Agora, as componentes em si */
for (i=0; i < resto; i++)
fscanf (logfile, "%lf\n",&autovetores[orbital+i][componente]);
}
}

void compara_pi(int* base, int* elemento, double** autovetores, double* autovalores, int nbasis, double threshold){
double sp[4], razao;
int orbital, componente;

printf("Comparando. . .\n");
printf("\n\n Orbital \tAutovalor\tpi/sigma\n");
for (orbital=0; orbital < nbasis; orbital++){
sp[0] = 0;
sp[1] = 0;
sp[2] = 0;
sp[3] = 0;
for (componente=0; componente < nbasis; componente++)
if(elemento[componente] == 6)
sp[(base[componente])] += autovetores[orbital][componente] * autovetores[orbital][componente];

razao = sqrt(sp[base[1]] + sp[base[2]] + sp[base[3]])/ sqrt(sp[base[0]]);
if ( razao >  threshold ){
printf("%8d\t%9g\t% 8g\n", orbital +1, autovalores[orbital], razao);
}
}
}


void imprime_matriz(int* base, int *elemento, double** autovetores, double* autovalores, int nbasis){
int orbital, componente;

for(componente=0; componente < nbasis; componente++){
printf("%d\t%d", elemento[componente], base[componente]);
for(orbital=0; orbital < nbasis; orbital++)
printf("\t% 8.6g", autovetores[orbital][componente]);
printf("\n");
}
}




void libera(FILE* logfile, double** autovetores, double* autovalores, int nbasis){
int i;

if (logfile != NULL)
fclose(logfile);
if(autovetores != NULL)
for (i = 0; i < nbasis; i++)
if (autovetores[i] != NULL)
free(autovetores[i]);
free(autovetores);

if (autovalores != NULL)
free(autovalores);
}
void *mallocX (unsigned int nbytes)
{
   void *ptr;
   ptr = malloc (nbytes);
   if (ptr == NULL) {
      printf ("Socorro! malloc devolveu NULL!\n");
      exit (EXIT_FAILURE);
   }
   return ptr;
}



[Penny Lane]

Joinville tem duas coisas, calor e chuva, quando não é um é outro... 

[FZapp]

Passei a noite em claro.

Enquanto meu computador (quase) incansavelmente mastigava números do meu doutorado, eu tentei ir dormir e quando encostei a cabeça no meu travesseiro, comecei a bolar algoritmos pra pegar os dados que importam do resultado das contas que tavam rodando.

Provavelmente seria mais fácil e rápido fazer em perl, mas como dizem que a melhor linguagem é a que a gente conhece,  aí está:

Muito bom seu programa, mas realmente pela manipulação de printf e fgets, seria algo ideal para o perl, um dia temos que nos juntar e te mostrar umas coisas...

[Fabi]

Que chato ter visto todos os episódios de lost e ter que ficar esperando até terça pra ver o episódio novo...

[Unknown]

Tá duro dormir com esse calor. Pelo menos não estou no RJ, que hoje teve a segunda maior sensação térmica do planeta, batendo até o deserto do Saara.

[Mr."A"]

madrugada de transição da quarta de cinzas para a quinta.... e eu acordado, vidrado no CC!

[Flavia]

2. Perdi o sono e to afundada em tédio.

[Mr."A"]

idem - pior que amanhã vou ter que acordar cedo. Opções:
I. durmo a hora que o sono chegar durante uns 30 minutos (é o que eu suponho que será o horário)
II. fico acordado direto e vou fazer o que eu tenho pra fazer.

[Sr. Bradley P. Richfield]

Estou usando as madrugadas pra ver as Olimpiadas de Inverno.

[FZapp]

Eu estou usando as madrugadas para que o emprego não fique impossível.


Sério.

[Fabi]

E o auge da noite foi o carro do sonho passando 2:30 da manhã...

[Fabi]

To brincando de controlar o tempo no BBB http://bbb.globo.com/

[Mussain!]

eu uma hora dessa......... por algum motivo sem sono.........

quem ta na mesma?