List of im() Examples

Examples of im()

com.nr.Complex.im()

Examples of com.nr.Complex.im()

      if (localflag) {
        fail("*** cisi: Incorrect symmetry for cosine integral");
        
      }


      localflag = localflag || zz1.im() != -zz2.im();
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** cisi: Incorrect symmetry for sine integral");
        
      }

View Full Code Here

Examples of com.nr.Complex.im()

      Complex r1 = Complex.I.mul(2.0*pi*5.0*n0/nn[0]);
      Complex r2 = Complex.I.mul(2.0*pi*7.0*n1/nn[1]);
      Complex r3 = Complex.I.mul(2.0*pi*9.0*n2/nn[2]);
      Complex r = r1.exp().mul(r2.exp()).mul(r3.exp());          
      data2[2*i]=r.re();
      data2[2*i+1]=r.im();
    }
//    System.out.printf(maxel(vecsub(data1,data2)));
    localflag = localflag || maxel(vecsub(data1,data2)) > sbeps;
    globalflag = globalflag || localflag;
    if (localflag) {

View Full Code Here

Examples of com.nr.Complex.im()

      for (j=0;j<nn[1]/2;j++) {
        Complex r1 = Complex.I.mul(2.0*pi*5.0*i/nn[0]);
        Complex r2 = Complex.I.mul(2.0*pi*7.0*j/nn[1]);
        Complex r = r1.exp().mul(r2.exp());          
        data2[i][2*j] = r.re();
        data2[i][2*j+1]= r.im();
        /*
        data2[i][2*j]= real(exp(2.0*pi*Complex(0.0,1.0)*5.0*double(i)/double(nn[0]))
          *exp(2.0*pi*Complex(0.0,1.0)*7.0*double(j)/double(nn[1])));
        data2[i][2*j+1]= imag(exp(2.0*pi*Complex(0.0,1.0)*5.0*double(i)/double(nn[0]))
          *exp(2.0*pi*Complex(0.0,1.0)*7.0*double(j)/double(nn[1])));

View Full Code Here

Examples of com.nr.Complex.im()

      }
      Complex r1 = Complex.I.mul(2.0*pi*5.0*i/nn[0]);
      Complex r2 = Complex.I.mul(pi*7.0);
      Complex r = r1.exp().mul(r2.exp());
      speq2[2*i]=r.re();
      speq2[2*i+1]=r.im();
      /*
      speq2[2*i]=real(exp(2.0*pi*Complex(0.0,1.0)*5.0*double(i)/double(nn[0]))*exp(pi*Complex(0.0,1.0)*7.0));
      speq2[2*i+1]=imag(exp(2.0*pi*Complex(0.0,1.0)*5.0*double(i)/double(nn[0]))*exp(pi*Complex(0.0,1.0)*7.0));
      */
    }

View Full Code Here

Examples of com.nr.Complex.im()

          Complex r1 = Complex.I.mul(2.0*pi*5.0*i/nn[0]);
          Complex r2 = Complex.I.mul(2.0*pi*7.0*j/nn[1]);
          Complex r3 = Complex.I.mul(2.0*pi*9.0*k/nn[2]);
          Complex r = r1.exp().mul(r2.exp()).mul(r3.exp());
          data4[i][j][2*k] = r.re();
          data4[i][j][2*k+1] = r.im();
          /*
          data4[i][j][2*k]=real(exp(2.0*Complex(0.0,1.0)*pi*5.0*double(i)/double(nn[0]))
            *exp(2.0*Complex(0.0,1.0)*pi*7.0*double(j)/double(nn[1]))
            *exp(2.0*Complex(0.0,1.0)*pi*9.0*double(k)/double(nn[2])));
          data4[i][j][2*k+1]=imag(exp(2.0*Complex(0.0,1.0)*pi*5.0*double(i)/double(nn[0]))

View Full Code Here

Examples of com.nr.Complex.im()

        Complex r1 = Complex.I.mul(2.0*pi*5.0*i/nn[0]);
        Complex r2 = Complex.I.mul(2.0*pi*7.0*j/nn[1]);
        Complex r3 = Complex.I.mul(pi*9.0);
        Complex r = r1.exp().mul(r2.exp()).mul(r3.exp());
        speq4[i][2*j]=r.re();
        speq4[i][2*j+1]=r.im();
        /*
        speq4[i][2*j]=real(exp(2.0*Complex(0.0,1.0)*pi*5.0*double(i)/double(nn[0]))
            *exp(2.0*Complex(0.0,1.0)*pi*7.0*double(j)/double(nn[1]))
            *exp(Complex(0.0,1.0)*pi*9.0));
        speq4[i][2*j+1]=imag(exp(2.0*Complex(0.0,1.0)*pi*5.0*double(i)/double(nn[0]))

View Full Code Here

Examples of com.nr.Complex.im()

            a[na][na]=q/a[nn][na];
            a[na][nn]=-(a[nn][nn]-p)/a[nn][na];
          } else {
            Complex temp=new Complex(0.0,-a[na][nn]).div(new Complex(a[na][na]-p,q));
            a[na][na]=temp.re();
            a[na][nn]=temp.im();
          }
          a[nn][na]=0.0;
          a[nn][nn]=1.0;
          for (i=nn-2;i>=0;i--) {
            w=a[i][i]-p;

View Full Code Here

Examples of com.nr.Complex.im()

            } else {
              m=i;
              if (wri[i].im() == 0.0) {
                Complex temp = new Complex(-ra,-sa).div(new Complex(w,q));
                a[i][na]=temp.re();
                a[i][nn]=temp.im();
              } else {
                x=a[i][i+1];
                y=a[i+1][i];
                vr=SQR(wri[i].re()-p)+SQR(wri[i].im())-q*q;
                vi=2.0*q*(wri[i].re()-p);

View Full Code Here

Examples of com.nr.Complex.im()

                vi=2.0*q*(wri[i].re()-p);
                if (vr == 0.0 && vi == 0.0)
                  vr=EPS*anorm*(abs(w)+abs(q)+abs(x)+abs(y)+abs(z));
                Complex temp= new Complex(x*r-z*ra+q*sa,x*s-z*sa-q*ra).div(new Complex(vr,vi));
                a[i][na]=temp.re();
                a[i][nn]=temp.im();
                if (abs(x) > abs(z)+abs(q)) {
                  a[i+1][na]=(-ra-w*a[i][na]+q*a[i][nn])/x;
                  a[i+1][nn]=(-sa-w*a[i][nn]-q*a[i][na])/x;
                } else {
                  temp= new Complex(-r-y*a[i][na],-s-y*a[i][nn]).div(new Complex(z,q));

View Full Code Here

Examples of com.nr.Complex.im()

                  a[i+1][na]=(-ra-w*a[i][na]+q*a[i][nn])/x;
                  a[i+1][nn]=(-sa-w*a[i][nn]-q*a[i][na])/x;
                } else {
                  temp= new Complex(-r-y*a[i][na],-s-y*a[i][nn]).div(new Complex(z,q));
                  a[i+1][na]=temp.re();
                  a[i+1][nn]=temp.im();
                }
              }
            }
            t=max(abs(a[i][na]),abs(a[i][nn]));
            if (EPS*t*t > 1)

View Full Code Here

TOP

All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.