News Features Documentation Download Participate |
## GSvit2D v 1.72D variant of GSvit is in a development phase and contains many untested algorithms. Its installed together with GSvit automatically and called the same way as 3D version, the only difference is that the parameter file must contain 2D dimensionality directive as a first item (to use this in 3D code is not mandatory, 3D dimensionality is a default option). Calculation is performed for transversally electric (TE) mode or transversally magnetic (TM) mode depending on user request (parameter TMMODE). Only subset of the algorithms provided for 3D is available now, due to fact that some methods are not working in 2D or are not yet implemented. Generally the parameter file is having very similar syntax to 3D version, however in some cases there are some small differences.
## Parameter file descriptionMain computation settings
Keywords are case sensitive. The following keywords can be used in the current version of GSvit (before issuing version 1.0 this means actual verison as obtained from SVN, but nearly all the parameters remain same for last issued binaries):
POOL
200 200 1e-6 1e-6
COMP
100
THREADS
4
- type of entity (integer): 2 - circle, 3 - triangle, 4 - rectangle
- point coordinates in pixel values (doubles): x y z values for each point (single for circle, two for rectangle and three for triangle, plus eventually radius (for circle)
- material type (0: standard material, 1: tabulated material, 2: Drude model, 3: CP3 model 4: CP model, 10: perfect electric conductor, 99: data from optical database) followed by:
- relative permittivity, relative permeability, electric and magnetic conductivity as double values for material type=0, and 1 (both representing linear material)
- epsilon, omega_p and nu for type=2 (Drude model)
- epsilon, sigma and three sets of a, phi, omega and gamma for type=3 (CP3, in development)
- epsilon, omega, gamma and two sets of a, phi, omega and Gamma for type=4 (CP)
- nothing for material type=10 (PEC)
- material name for material type=99, e.g. Al2O3
so the entry looks for example as (4 100 100 20.5 0 22.13 1 0 0) for circle with radius 20.5 pixels, position (100, 100) and relative permittivity of 22.13 (rest or material properties is same as for vacuum). Note that at present version of the 2D code the CP, CP3 model and Drude model is not implemented yet.
**none**which means no boundary treatment providing reflection similar to perfect electric conductor,**liao**which means 2nd order absorbing (Liao) providing quite good absorption.
Note that for generated source the orientation should match the
mode of calculation (transversally electric or transversally magnetic.
So for example setting
fits to transversally magnetic mode, e.g.
Similarily to do a source for transversally electric calculation, write SOURCE_POINT 50 50 1 0.633e-6 1 1.57079632679 1.57079632679
and then
Note that in present version the TSF source is in vacuum, material cannot cross it. Orientation of axes of incoming wave comes from 3D version, assuming that theta=90 degrees. It is shown in the following image, table shows some typical useful values of parameters:
Note that TSF is not yet implemented on GPU. Output
Note that in present version every output needs GPU data to be synchronized with CPU, which can take significant time. For really fast simulations try to reduce frequency of data outputs.
Note that Gwyddion shows data with top-left corner being center of coordinates, orientation of axes on what is seen in Gwyddion is show below (for 2D code we have always the z constant condition met, so axes are oriented as seen at the right image)
Graphics card use
General commands
Set text output (step by step), from full (4) to silent mode (0). ## ExamplesExamples are in preparation. Check summary of some of tested parameter files for 2D problems. (c) Petr Klapetek, 2013 |