View Javadoc
1   /*
2    * Licensed to the Hipparchus project under one or more
3    * contributor license agreements.  See the NOTICE file distributed with
4    * this work for additional information regarding copyright ownership.
5    * The Hipparchus project licenses this file to You under the Apache License, Version 2.0
6    * (the "License"); you may not use this file except in compliance with
7    * the License.  You may obtain a copy of the License at
8    *
9    *      https://www.apache.org/licenses/LICENSE-2.0
10   *
11   * Unless required by applicable law or agreed to in writing, software
12   * distributed under the License is distributed on an "AS IS" BASIS,
13   * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14   * See the License for the specific language governing permissions and
15   * limitations under the License.
16   */
17  package org.hipparchus.special.elliptic.jacobi;
18  
19  import org.hipparchus.complex.FieldComplex;
20  import org.hipparchus.dfp.Dfp;
21  import org.hipparchus.dfp.DfpField;
22  import org.junit.jupiter.api.Test;
23  
24  import java.io.IOException;
25  
26  import static org.junit.jupiter.api.Assertions.assertEquals;
27  
28  class JacobiEllipticFieldComplexTest {
29  
30      @Test
31      void testComplex() throws IOException {
32          final DfpField field = new DfpField(30);
33          final FieldComplex<Dfp> m = new FieldComplex<>(field.newDfp("0.3"), field.newDfp("1.0"));
34          FieldJacobiElliptic<FieldComplex<Dfp>> je = JacobiEllipticBuilder.build(m);
35          final FieldComplex<Dfp> z = new FieldComplex<>(field.newDfp("5.2"), field.newDfp("-2.5"));
36          final FieldCopolarC<FieldComplex<Dfp>> valuesC = je.valuesC(z);
37          assertEquals(-0.24609405083573348938, valuesC.sc().getRealPart().getReal(),      1.0e-15);
38          assertEquals( 0.74202229271111558523, valuesC.sc().getImaginaryPart().getReal(), 1.0e-15);
39      }
40  
41  }