Package com.nr.test.test_chapter10

Source Code of com.nr.test.test_chapter10.Test_Frprmn$Funcd_Frprmn

package com.nr.test.test_chapter10;

import static com.nr.NRUtil.*;
import static java.lang.Math.*;
import static org.junit.Assert.fail;

import org.junit.After;
import org.junit.Before;
import org.junit.Test;

import com.nr.RealValueFunWithDiff;
import com.nr.min.Frprmn;
import com.nr.sf.Bessjy;

public class Test_Frprmn {

  @Before
  public void setUp() throws Exception {
  }

  @After
  public void tearDown() throws Exception {
  }

  @Test
  public void test() {
    int i,N=20;
    double theta,phi,f0,f,d0,d,pi=acos(-1.0);
    double sbeps1=1.0e-7,sbeps2=5.0e-16;
    double[] p = new double[3];
    boolean localflag=false, globalflag=false;

   

    // Test Frprmn
    System.out.println("Testing Frprmn");

    Funcd_Frprmn fd = new Funcd_Frprmn();
    Frprmn frprmn = new Frprmn(fd);
    p[0]=p[1]=p[2]=0.1;
    p=frprmn.minimize(p);
    f0=frprmn.fret;
    d0=sqrt(SQR(p[0])+SQR(p[1])+SQR(p[2]));

    for (i=0;i<N;i++) {
      theta=pi*i/N;
      phi=pi*i/N;
      p[0]=sin(phi)*cos(theta);
      p[1]=sin(phi)*sin(theta);
      p[2]=cos(phi);
      p=frprmn.minimize(p);
      f=frprmn.fret;
      d=sqrt(SQR(p[0])+SQR(p[1])+SQR(p[2]));
//      System.out.println("     " << abs(d-d0) << " " << abs(f-f0));

      localflag = localflag || abs(d-d0) > sbeps1;
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** Frprmn: First minimum of radial function reported at different radius for different starting points");
       
      }

      localflag = localflag || abs(f-f0) > sbeps2;
      globalflag = globalflag || localflag;
      if (localflag) {
        fail("*** Frprmn: Reported function value at first minimum is different for different starting points");
       
      }
    }

    if (globalflag) System.out.println("Failed\n");
    else System.out.println("Passed\n");
  }
 
  class Funcd_Frprmn implements RealValueFunWithDiff{
    Bessjy b = new Bessjy();
    double u;
    public double funk (double[] x) {
      return(b.j0(sqrt(SQR(x[0])+SQR(x[1])+SQR(x[2]))));
    }
    public void df(double[] x, double[] deriv) {
      u=sqrt(SQR(x[0])+SQR(x[1])+SQR(x[2]));
      deriv[0]= -x[0]*b.j1(u)/u;
      deriv[1]= -x[1]*b.j1(u)/u;
      deriv[2]= -x[2]*b.j1(u)/u;
    }
  }
}
TOP

Related Classes of com.nr.test.test_chapter10.Test_Frprmn$Funcd_Frprmn

TOP
Copyright © 2018 www.massapi.com. All rights reserved.
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.