|
| 1 | +import java.lang.*; |
| 2 | + |
| 3 | +import imasjava.*; |
| 4 | +import imasjava.wrapper.LowLevel; |
| 5 | + |
| 6 | +public class TestGyrokineticsValidate { |
| 7 | + |
| 8 | + public static void main(String[] args) { |
| 9 | + imas.gyrokinetics_local ids = new imas.gyrokinetics_local(); |
| 10 | + |
| 11 | + System.out.println("### Validating Gyrokinetics validate simplest case...\n"); |
| 12 | + try { |
| 13 | + ids.validate(); |
| 14 | + } catch (Exception e) { |
| 15 | + } |
| 16 | + |
| 17 | + ids.ids_properties.homogeneous_time = LowLevel.IDS_TIME_MODE_HOMOGENEOUS; |
| 18 | + System.out.println("### Validating HOMOGENEOUS time array...\n"); |
| 19 | + try { |
| 20 | + ids.validate(); |
| 21 | + } catch (Exception e) { |
| 22 | + System.out.println("Error, HOMOGENEOUS time array exception " + e.getMessage()); |
| 23 | + } |
| 24 | + |
| 25 | + System.out.println("### Validating species_all...\n"); |
| 26 | + ids.species_all = new imas.gyrokinetics_local.species_allClass(); |
| 27 | + ids.species_all.angle_pol_equilibrium = new Vect1DDouble(1); |
| 28 | + ids.species_all.angle_pol_equilibrium.setElementAt(0, 10.0); |
| 29 | + try { |
| 30 | + ids.validate(); |
| 31 | + } catch (Exception e) { |
| 32 | + System.out.println("Error, species_all validate exception " + e.getMessage()); |
| 33 | + } |
| 34 | + |
| 35 | + System.out.println("### Validating species...\n"); |
| 36 | + ids.species = new imas.gyrokinetics_local.speciesClass[2]; |
| 37 | + for (int i = 0; i < 2; i++) { |
| 38 | + ids.species[i] = new imas.gyrokinetics_local.speciesClass(); |
| 39 | + ids.species[i].potential_energy_norm = new Vect1DDouble(1); |
| 40 | + ids.species[i].potential_energy_norm.setElementAt(0, 10.0); |
| 41 | + } |
| 42 | + try { |
| 43 | + ids.validate(); |
| 44 | + } catch (Exception e) { |
| 45 | + System.out.println("Error, species validate exception " + e.getMessage()); |
| 46 | + } |
| 47 | + |
| 48 | + System.out.println("### Validating non_linear...\n"); |
| 49 | + ids.non_linear = new imas.gyrokinetics_local.non_linearClass(); |
| 50 | + ids.non_linear.fluxes_3d.particles_phi_potential = new Vect3DDouble(0,0,0); |
| 51 | + |
| 52 | + try { |
| 53 | + ids.validate(); |
| 54 | + } catch (Exception e) { |
| 55 | + System.out.println("Error, non_linear validate exception " + e.getMessage()); |
| 56 | + } |
| 57 | + |
| 58 | + System.out.println("### Validating linear...\n"); |
| 59 | + ids.linear = new imas.gyrokinetics_local.linearClass(); |
| 60 | + |
| 61 | + ids.linear.wavevector = new imas.gyrokinetics_local.linearClass.wavevectorClass[1]; |
| 62 | + for(int i = 0; i < 1; i++) { |
| 63 | + ids.linear.wavevector[i] = new imas.gyrokinetics_local.linearClass.wavevectorClass(); |
| 64 | + ids.linear.wavevector[i].eigenmode = new imas.gyrokinetics_local.linearClass.wavevectorClass.eigenmodeClass[1]; |
| 65 | + for (int j = 0; j < 1; j++) { |
| 66 | + ids.linear.wavevector[i].eigenmode[j] = new imas.gyrokinetics_local.linearClass.wavevectorClass.eigenmodeClass(); |
| 67 | + ids.linear.wavevector[i].eigenmode[j].poloidal_turns = 5; |
| 68 | + ids.linear.wavevector[i].eigenmode[j].angle_pol = new Vect1DDouble(2); |
| 69 | + ids.linear.wavevector[i].eigenmode[j].time_norm = new Vect1DDouble(2); |
| 70 | + |
| 71 | + ids.linear.wavevector[i].eigenmode[j].moments_norm_gyrocenter_bessel_1 = new imas.gyrokinetics_local.linearClass.wavevectorClass.eigenmodeClass.moments_norm_gyrocenter_bessel_1Class(); |
| 72 | + ids.linear.wavevector[i].eigenmode[j].moments_norm_gyrocenter_bessel_1.pressure_perpendicular = new Vect3DComplex(2,2,0); |
| 73 | + ids.linear.wavevector[i].eigenmode[j].moments_norm_gyrocenter_bessel_1.density = new Vect3DComplex(2,2,0); |
| 74 | + } |
| 75 | + } |
| 76 | + try { |
| 77 | + ids.validate(); |
| 78 | + } catch (Exception e) { |
| 79 | + System.out.println("Error, linear validate exception " + e.getMessage()); |
| 80 | + } |
| 81 | + |
| 82 | + System.out.println("### Gyrokinetics validate test completed\n"); |
| 83 | + } |
| 84 | +} |
0 commit comments