News
Features
Documentation
Application examples
FDTD tutorial
Download
Participate

GSvit v 1.8 graphical user interface

GSvit architecture is based on use of computational core (GSvit) and a GUI (XSvit). This page provides documentation for the graphical user interface, you can also check the documentation for computational core.

Installation

Please see the download page.

Running the graphical user interface

Quick start

Check the basic program operation example to run your first computation. Other examples are linked at bottom of this page.

Documentation

XSvit is primarily a viewer of parameter files suitable for GSvit FDTD calculations. It cannot do anything except it, however it can simplify the process of parameter file creation and modification as it serves as a viewer and basic editor of parameter files. Its GUI is shown in the following image:

Technically the primary data in XSvit are the parameter file that is shown in the "Par. file" tab in left side of the screen. This represents the data directly as are used in GSvit. While loading this file or reloading it after any modification, XSvit goes through the parameters in it and constructs list of entries shown in tab "Parameters". By clicking on any of these entries, or by using the Edit parameters menu, we can adjust parameters or add some more of them. If we do so, parameter file is automatically written into a temporary location and reloaded back. This is to guarantee consistency of the parameter file with GSvit - the parameter file is still the only source of information and you can run GSvit solver commandline on the parameter file, or edit the parameter file in any text editor, leaving out XSvit viewer completely.

As an alternative to clicking on menu items or parameter list, you can also directly edit content of "Par. file" tab multiline text entry. If we do so, parameter file is dynamically reloaded as we type (and all the entries in "Parameter" tab are again refreshed).

Note that the parameter file used within all the modifications is only a temporary file, so it must be saved prior to computation or exiting the software.

Similarily to parameter file, there is also material file and list of material entries representing "vector material" if such is used in the computation. Again, material file is the primary source of information and list of entries is only to assist the user dislaying it. At this version there is no GUI to edit the material file, it can be done only directly in "Mat. file" tab, according to GSvit syntax rules described in GSvit documentation.

Which material file to show is governed by parameter file entry for "vector material". Therefore, chaning this entry also leads to loading another material file. Current temporary material file is saved automatically when we save the parameter file.

Main screen

There are few key elements of the main screen that appear when you start XSvit:

  • Parameters list and parameter file editing tool (each in separate tab). This is the place where most of the settings are adjusted, either clicking on parameter list entries or directly changing the parameter file content.
  • Material objects list and material file editing tool (each in separate tab). Here objects entered as geometric entities are input (cubes, spheres, tetraheral meshes), which is easiest way of including some objects into computational volume. This is however not the only way how to include material, you could also predefine material properties in every voxel directly.
  • 3D display window showing all the selected objects (checkboxes are in parameter and material lists). If there is an image output being saved during the calculation, it can be also visualised here.
  • Graph area to show text output files to see the computation progress.
  • GSvit output window showing the messages output from GSvit when it is started.
  • Statusbar at the bottom, showing different error messages and confirmations.

Note that at this version of XSvit, some amount of warnings and errors is still output to console, namely those related to parsing the parameter file (that is done by GSvit functions where a console output is normally used). This will change in future versions to get rid of the console (which is too unusual on Windows).

Dialogues

There are many simple dialogues that can be used for entering different computation parameters. They can be invoked by double clicking on an entry in list of parameters at left side of the screen (this works only for parameters, not for material objects in this version), or by Edit parameters menu.

To delete an entry in this version, you need to go to the parameter or materials file tab and delete it.

Most of the dialogues are shown in some of the examples (see bottom of this page) in practice, here is their list as well. Note that many more details about realisation of different parameters are in computational core documentation.


Parameter settings

Global parameters that can be set are location of GSvit and Gwyddion executable. These are used for calling GSvit and Gwyddion directly from GUI, if we only want to create parameter files and run the computation and output visualisation elsewhere, we don't need them. Please note that on Windows OS the default installation path for GSvit binary is typically C:\Program Files (x86)\GSvit\bin\gsvit.exe.


Computational domain settings

There are two main computational domain parameters: voxel size (also called resolution in the documentation) and voxel spacing. Voxel size can be any (for any direction), however note that GSvit is 3D FDTD solver, so settings any dimension to e.g. 1 does not convert it to 2D solver. Voxel spacing should be the same for all the dimensions in this version, using different spacings is untested, even if in principle FDTD allows it. If we are in vacuum, voxel spacing should be at maximum 1/10 of the wavelength used.


Basic parameters settings

Several basic parameters influence the whole computation (all the algorithms): number of Yee algorithm steps (real time step is calculated automatically), verbosity level (0...3), use of multithreading based on OpenMP (-1...N) where -1 means use of all the available cores, and use of any of the available graphics card instead of processor. Note that verbose level should not be 0 to see anything in the text output (e.g. in log window in the GUI). On the other hand, high verbosity level might lead to some delays in GUI visualisation on MS Windows, so value of 1 seems to be optimum. If you have doubts whether graphics card is usable, you can run gsvit test 0 to perform basic test of Nvidia CUDA capable devices. Cards are counted from zero.


Boundary conditions

Computational volume typically needs to be terminated somehow and the following dialogue is used to set the boundary conditions. Different options can be chosen, starting from nothing (which reflects waves) up to Convolutional Perfectly Matched Layer conditions (for details see CPML boundary condition example). Note that even if in principle any boundary can have different boundary condition applied, it is not guaranteed that all the combinations (e.g. Liao + CPML) are working together properly.

This dialog can be used also to set some boundary conditions within the material, now only peridic boundary conditions. To do so, you must specify where from and where to the space should be treated as periodic.


Input of geometry of media

There are two ways how to input objects in this version of GSvit. These can be defined voxel-by-voxel for all the voxels in the computational domain or defined using discrete objects, here also called as vector data (spheres, cubes,... up to tetrahedral meshes). Voxel-by-voxel data are not further visualised in XSvit, vector data are loaded automatically on basis of settings in this dialogue to the bottom widget showing material entities.


Output settings

Main general output setting is the output file where all the cross-sections (images) go to. It is Gwyddion file so it can be opened only by Gwyddion.


Point source

A simple dipole-like point source realised by setting field quantities in a single cell to a prescribed time dependency. We need to set its voxel position and waveform, that can be both loaded from a text file (containing total number of steps followed by sequence of these entries: index Ex Ey Ez Hx Hy Hz). This file ca be also automatically generated providing a sine wave source or sine wave damped by gaussian envelope. Finally, two parameters theta and phi determine the source orientation.


Scattered field source

Scattered field source creates a plane wave by introducing its oposite to every face of a perfectly electric conductor. To use it you first need to have the conductor somewhere in your computational domain, and then also set the orientation of wave by the following dialogue. Other source conditions are same as for point source. For orientation details see the computational core documentation


Total/scattered field source

Plane wave can be also injected on some virtual boundary inside computational domain and removed at the other side. This approach is calle Total/scattered field source. To setup it we use the same parameters as to scattered field source, adding only definition of the virtual boundary where it is added/subtracted.


Focused total/scattered field source

This represents z direction oriented focused wave source using total/scattered field formulation with multiple plane waves decomposition according to algorithm published in I.R.Capoglu, A. Taflowe, V. Backman, Optics Express 23 (2008) 19208. In addition to definition of virtual boundary where it is added/subtracted, we need to define maximal angle of incoming light, focal distance in voxel units and wave polarisation angle. Integer parameters n and m determine number of plane waves to be integrated (more is better, values around 12, 36 should be already enough).


Layered total/scattered field source

Plane wave source using layered total/scattered field formulation as published by I.R. Capoglu, G.S. Smith: IEEE Transactions on Antennas and Propagation 02/2008, 56:1. In contrast to TSF source, here the material can be formed by N layers in z-direction (in present implementation only dielectric and non-absorbing). Incident wave can cross the sample at angle as well. Layers need to be introduced both in parameters of the source and in material description files and the values of material parameters need to match at LTSF boundary. In present GUI only two interfaces can be set (GSvit itself could handle any number of layers). Note that the source starts always in vacuum, so to add single material layer we need exactly two interfaces (one for upper boundary and one for lower boundary).


Layered and focused total/scattered field source

A mixture of TSFF and LTSF source, i.e. focused source that can travel accross layered media. GUI is a combination of settings from both the sources (see above).


Point output

A simple point output providing text file with specified data extracted from specified position. It does not need to be generated at every step (in some steps it can be skipped), which is useful namely when working with GPU (all the outputs lead to possibly unnecessary communication between graphics card and RAM).


Image output

Image output going to general output file. It always makes a cross-section, so its orientation and direction can be specified. Description is a string shown for the channel in Gwyddion data browser.


Summing output

Output sum of electric field intensity (components: ex, ey, ez, all) or absorption (component: abs) extracted from bounding box and material with given properties (epsilon only used for this). If epsilon is -1, no control of this parameter is done and whole box is used. Skipping controls frequency of file output only (and sync with GPU eventually), values are calculated at each time step anyway and will be output for intermediate steps as well. At present this is tested only on CPU.


Volume output

Output volume information about field components, material components (components: ex, ey, ez, hx, hy, hz, all, material) or absorption (component: abs). Skipping controls frequency of file output only (and sync with GPU eventually), ouput is performed only between start and stop step (-1 in these parameters means starting in first step and stopping in last one). Volume information can be binary or ASCII. Binary file is formed by 3D array of doubles that can be visualised e.g. by Paraview. ASCII output starts with header informing about voxel dimensions and selected component.


Force output

Outputs optical force acting on volume defined by six integers. X, y, and z component time dependence is output together with its values averaged over source period (source frequency is determined using FFT). Optical force is calculated using Maxwell stress tensor. No media should cross the boundary and scaterrer to be evaluated should be placed inside the evaluation volume.


Near-to-far field transform basic settings

Defines integration boundary for calculating far field values using time domain Near field to far field transform (NFFF). Sometimes it is useful to skip some part of boundary (or whole boundary) even if it is non-physical. In this case you can enter the area which should be skipped from the boundary (for every boundary independently).


Near-to-far field output point

Sets the position of far field integration point using NFFF.


Near-to-far field output area

Sets the generated set of far field points. Points are part of a sphere with defined voxel radius, with span given by theta and phi angles.


Periodic near-to-far field transform basic settings

Periodic NFFF is designed to be used with periodic boundary conditions and integrates the NFFF values from predefined number of repetitions of calculated far field values. It can be therefore used e.g. to calculated field from a grating while field over single pit is only calculated by Yee's algorithm.


Periodic near-to-far field output point

Sets the position of far field integration point using PNFFF.


Periodic near-to-far field output area

Sets the generated set of periodic far field points. Points are part of a sphere with defined voxel radius, with span given by theta and phi angles.

Material objects

To add material objects into computational domain you can either specify the material properties voxel-by-voxel and load it, or, which is more common, to use different individual objects (cubes, spheres, tetrahedral meshes, ...) specified in vector material file. To do this within GUI, use menu item Edit material objects.

As an example, the following dialog is used to add sphere:

Note that the position and radius is given in voxel coordinates, so the real position and dimensions is dependent on your voxel spacing. You have wide variety of choices of material properties that will be assigned to the object (sphere); here a perfect electric conductor was chosen. For a simple dielectric object, like glass sphere, the best choice is "tabulated properties" checkbox and some reasonable parameters of relative electric permittivity, unit magnetic permeability and zero conductivities (both electric and magnetic).

Another way how to add material objects within GUI is to write the vector material file directly. This means writing entities like

4 20 20 20 50 10

which means adding a perfect electric conductor sphere of radius 50 voxels to position 20,20,20. For details on materials see computational core documentation or some of the examples.

Examples

There are several examples distributed together with the application. to test different software algorithms. Please note that if you have installed the Gsvit package into the standard location (e.g. /usr/share/gsvit/tests on Linux), the software runned as ordinary user won't have permissions to write the results. Therefore, we recommend to copy the respective test directory to the home directory owned by the user and run the computation from this location.

You can also use the following step-by-step guides showing how to setup some basic computational parameters:

Basic program operation
Total/scattered field source, PEC sphere
Scattered field source, PEC sphere
Near-to-far field transformation
CPML boundary condition
Force acting on a non-spherical particle
Diffraction by a transmission grating
Diffraction by a reflection grating defined as height field
Local field enhancement close to a metallic rod







(c) Petr Klapetek, 2013