Todo

Todo
dia

quarta-feira, 23 de abril de 2014

prece #1

por sentido
outras vezes tanta
vez dez a potência
razão necessariamente
suficiente

queremos inferir a
consequência implicada
figura colimada
forma polarizado feixe
ressonante seu
modo harmônico
deus

e a vossa glória
anunciada, profetizada, proferida,
negociada e preferida à ira da divindade
no inverno do deserto, sob a tenda de tapetes,
ao preço do sexo da posse da pedra;
para a resolução das vontades
para a renovação dos costumes,
para
o esquecimento


quarta-feira, 9 de abril de 2014


#include
#include
#include

#define NC 10000000
#define dx 0.0001
#define ntime 600
#define nsamples 100000
#define nker 10000
#define ntransient 15000
#define IA 16807
#define IM 2147483647
#define AM (1.0/IM)
#define IQ 127773
#define IR 2836
#define NTAB 32
#define NDIV (1+(IM-1)/NTAB)
#define EPS 1.2e-7
#define RNMX (1.0-EPS)
#define TTT 100

float ran1(long *idum)

{
int j;
long k;
static long iy=0;
static long iv[NTAB];
float temp;
if (*idum <= 0 || !iy) { 
if (-(*idum) < 1) *idum=1
else *idum = -(*idum);
for (j=NTAB+7;j>=0;j--) { 
k=(*idum)/IQ;
*idum=IA*(*idum-k*IQ)-IR*k;
if (*idum < 0) *idum += IM;
if (j < NTAB) iv[j] = *idum;
}
iy=iv[0];
}
k=(*idum)/IQ
*idum=IA*(*idum-k*IQ)-IR*k; if (*idum < 0) *idum += IM
j=iy/NDIV;
iy=iv[j]; 
iv[j] = *idum; 
if ((temp=AM*iy) > RNMX) return RNMX; else return temp;
}

float gasdev(long *idum)
{
float ran1(long *idum);
static int iset=0;
static float gset;
float fac,rsq,v1,v2;
if (*idum < 0) iset=0
if (iset == 0) { 
do {
v1=2.0*ran1(idum)-1.0;v2=2.0*ran1(idum)-1.0; rsq=v1*v1+v2*v2; 
} while (rsq >= 1.0 || rsq == 0.0); 
fac=sqrt(-2.0*log(rsq)/rsq);
gset=v1*fac;
iset=1
return v2*fac;
} else
iset=0
return gset; 
}
}


int main (void) {
int i, kn, km, ke, ki, kij;
float p[NC], ps[NC], xmin, xvmin, xmax, xvmax, jmin, jmax, alfa, sv, gamma, we, ca,cn, cj, xll, tj, dt, dt2, dt3, dts, cam, w, idum;
float cnm, gammam, cjm, cd, ncoarse, ntimenoise, xm, xm2, xm3, xm4, xv, xv2, xv3,xv4, jam, jam2, jam3, jam4, t1, t2, t3, t4, tk;
float xx1, xx10, xxv1, xxv10, v1, v10, v101;
float jitot1[ntime], jdtot1[ntime], energy [ntime], correlacaov1[nsamples][1000], correlacao;
float xv1[nker];
for(i=0;i<NC;i++)
{
p[i]=0;
ps[i]=0;
}
xmin = - NC * dx /2.0;
xvmin = - NC * dx /2.0;

xmax = NC * dx /2.0;
xvmax = NC * dx /2.0;

jmin = - NC * dx /2.0;
jmax = NC * dx /2.0;
  
idum = -798764321;
alfa = 5.0;
      sv = 1.0;

  we = 1.0;
  gamma = 1.0;
  ca = 1.0;
  cn = 0.0;
          cj = 10.0;


      xll = 20.0;
      tj =  25.0;
 
  dt = 0.0001;
  dt2 = dt * dt;
      dt3 = dt2 * dt;
      dts = sqrt(dt);
 
  cam = ca / we;
      cnm = cn / we;
      gammam = gamma / we;
      cjm = cj / we;

      cd = (1.0 + gammam * dt / 2.0);       

      ncoarse = floor(0.005/dt);

      ntimenoise = 10 * floor((1.0 / alfa) / 0.005);

  xm = 0.0;
  xm2 = 0.0;
  xm3 = 0.0;
  xm4 = 0.0;

          xv = 0.0;
          xv2 = 0.0;
          xv3 = 0.0;
          xv4 = 0.0;

          jam  = 0.0;
          jam2 = 0.0;
          jam3 = 0.0;
          jam4 = 0.0;
 
  for(kn=0;kn<nsamples; kn++)
  {
  for(km=0;km<nsamples;km++)
  {
  correlacaov1[kn][km]=0;
  }
  }

  for(i=0;i<ntime;i++)
  {
jitot1[i] = 0.0;
jdtot1[i] = 0.0;
                energy[i] = 0.0;        
}
for(i=0; i<nsamples; i++)
{
v1 = 0.0;
         for (kn = 0; kn
{
           for(km = 0;km
  {
               w = gasdev(&idum);
               v1 = v1 - alfa * v1 * dt + sv * w * dts;
   
            }
         }
}
  xxv1 = 0.000;
      xxv10 = 0.0;

      xx1 =  0.0;
      xxv1 = 0.0;

      for(ke=0; ke<nker; ke++)
         xv1[ke] = 0.0;
for(ki = 0; ki< ntransient;ki++)
{
for(kij = 0; kij< ncoarse; kij++)
{
correlacaov1[ki][kij]=v1;
           v10 = v1;
           w = gasdev(&idum);
           v1 = v1 - alfa * v1 * dt + sv * w * dts;
           v101 = sqrt(dt)*(v10+v1)/(2.0)/we;

           xxv10 = xxv1;
           xxv10 = xv1[0];
  xx10 = xx1;



           t1 = 0.0;
           t2 = 0.0;
           t3 = 0.0;
           t4 = 0.0;
           tk = 0.0;

           for(ke = 0; ke< nker; ke++)
              tk = tk + exp(-(ke-1) * dt * alfa) * xv1[ke];

           t1 = (1.0 - gammam * dt) * xxv10;
           t1 = xxv10 - gammam * dt * tk;
           t2 = - (cam + cjm) * xx10 * dt;
           t3 = - cnm * xx10 * xx10 * xx10 * dt;          
           t4 = v101;            

           xxv1 = t1 + t2 + t3 + t4;

           for(ke = 0; ke< nker-1;ke++)
              xv1[ke+1] = xv1[ke];

           xv1[1] = xxv1;              

           xx1 = xx10 + (xxv1+xxv10)*dt/2.0;

}
}
correlacao=0;

for(ki = 0; ki< ntransient;ki++) /*dado TTT */
{
for(kij = 0; kij< ncoarse; kij++)
{
correlacao= correlacao + correlacaov1[ki][kij+TTT]*correlacaov1[ki][kij];
}
correlacao= correlacao/ncoarse;
}
  printf("AAAAAAAA\n\n\n\n %f\n\n\n\n", correlacao);
    return 0;
}