<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Visualizing methods of multi-criteria alternatives for pairwise comparison procedure</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>A.A. Zakharova</institution>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Bryansk State Technical University</institution>
          ,
          <addr-line>Bryansk</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>Keldysh Institute of Applied Mathematics Russian Academy of Sciences</institution>
          ,
          <addr-line>Moscow</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>The article deals with the problem of choosing a preferred alternative in a pairwise comparison procedure. The difficulties of applying this procedure in a case of using alternatives with a large number of criteria are noted. It is proposed to supplement the procedure of expert pairwise comparison with visualization tools of multi-criteria alternatives. The paper considers several visualization methods for multi-criteria alternatives for pairwise comparison procedures: histograms, two-dimensional graphs, three-dimensional surfaces, probability distribution diagrams, visualization based on modifications of radar and radial diagrams, as well as combined methods. It described an experimental study of the application of the considered method for the task of determining the preferred alternative by the example of choosing one of two OpenFoam solvers (rhoCentralFoam and pisoCentralFoam), with the help of which estimates of the accuracy of calculating the inviscid flow around a cone were obtained. Еach solver is characterized by 288 criteria. It is shown that the use of some of the methods considered does not make it possible for the expert to make a choice. In this case, a good result was obtained using methods for constructing three-dimensional surfaces, probability distribution diagrams, as well as using the combined method based on modified radar diagrams. It is concluded that the rhoCentralFoam solver is more preferable if there are no additional criteria for ranking the criteria. The possibility of using the combined method in combination with the ranking procedure of criteria (or their groups) during decision-making is also noted.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        In decision theory, one of the basic tasks is ranking
alternatives [
        <xref ref-type="bibr" rid="ref1 ref2">1, 2</xref>
        ]. This allows setting their priority in
relation to the task at hand. There are various methods
that allow such ranking, some of which are expert.
Among the expert ranking methods, the method of
pairwise comparisons [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] has proven itself well, the
essence of which is to provide the expert alternately with
pairs of alternatives for comparison, during which he
prefers one of them. In particular, this procedure used in
one of the classical decision-making methods — the
hierarchy analysis method developed by T. Saati [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ], in
which it is necessary to construct matrices of pairwise
comparisons for all levels of the hierarchy.
      </p>
      <p>
        The use of pairwise comparisons is usually effective
in cases where each alternative is well reflected in the
expert’s perception or is characterized by a small number
of criteria (usually no more than 10) [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. In the case when
an expert needs to compare new for him alternatives with
a large number of criteria, this can cause difficulties for
him. Therefore, in such situations, it is necessary to use
additional tools, for example, reducing the dimension, or
statistical processing of criteria values. However, even
applying these approaches, there is still a chance of not
getting the desired result. For example, in the case of
calculating statistical characteristics, we can get
conflicting data in a situation where the mathematical
expectation for an alternative is better, but the variance is
worse. Therefore, additional tools are needed that could
help the expert decide.
      </p>
      <p>
        One such methods may be visual analytics - when for
each alternative a corresponding visual image is
constructed that characterizes the set of values of its
criteria. Visualization is able to present the alternative as
a holistic image, and it will be easier for an expert to
make his choice with its help. It should be noted that the
visual comparison of alternatives is currently already
being applied and shows a good result, for example,
when comparing the site design at the stage of its design
[
        <xref ref-type="bibr" rid="ref6">6</xref>
        ].
      </p>
      <p>Thus, we have the task of visualizing data sets
characterizing alternatives. Consider and analyze several
approaches and methods that can be applied to visualize
multi-criteria alternatives.
2.</p>
    </sec>
    <sec id="sec-2">
      <title>Visualization methods</title>
      <sec id="sec-2-1">
        <title>Data preparation</title>
        <p>
          The visualization procedure begins with a step
requiring initial data preparation.
1. All input values should be given to a numeric format
- relevant methods of decision theory may be used
for these purposes [
          <xref ref-type="bibr" rid="ref7">7</xref>
          ].
2. Normalization of data per segment [0; 1] taking into
account the direction of the criterions optimization
(maximization or minimization). For these purposes,
the formulas can be used:
•  ′, =     , −, −    ,,  – in case of maximization;
•  ′, =     , −,−   ,,  – in case of minimization,
where i – alternative number (1 ≤  ≤  ), j – criteria
number (1 ≤  ≤  ), N – count of alternatives, K – count
of criteria, vmax,j, vmin,j – maximum and minimum possible
value of j-th criteria. Values vmax,j, vmin,j usually
determined on the basis of their physical meaning.
However, if there are problems with their definition, then
they can be calculated by the formulas:
   , = min   , ,
        </p>
        <p>, = max   , .</p>
        <p>After the data is prepared, you can proceed to
visualize them. For these purposes, several different
approaches and methods can be applied, in each of which
we will consider the visualization of two alternatives and
the features of their visual pairwise comparison.</p>
        <p>One of the most accessible and simple methods of
visualization of multidimensional data is diagrams, and
among them, the most accessible in the procedure of
pairwise comparison can be called histograms. In this
case, two main approaches to their application can be
distinguished.
1. Both alternatives are shown on a common
histogram. (Fig. 1). In addition, you can use the
option of overlapping columns to focus on the
deviation of values by criteria.</p>
        <p>Each alternative
histograms. (Fig. 2).</p>
        <p>is
visualized
on
separate</p>
        <p>Note that in the 2nd approach, it is important to use
identical parameters of the diagrams (color, size, etc.) in
order to reduce the subjectivity of the perception of
visual images, describing the data sets of the
corresponding alternatives, and to enable the expert to
focus on a holistic perception of visual images.</p>
        <p>When determining preferences among alternatives
based on visual images, an expert can be guided by his
perception of the degree of filling of the columns of
diagrams, including and the total area of all columns (the
second approach is more suitable for these purposes) or
the subjectively averaged deviation of the criteria
columns of the two alternatives (the first approach is
more suitable for these purposes).</p>
        <p>The use of histograms is justified for those situations
when it is necessary to see the whole alternative as a
whole, and the number of criteria is not too large (about
5-20). Moreover, the effectiveness of this visualization
method is additionally achieved if ranking approaches or
non-linear scales were used in the normalization process
since this allows you to get quite noticeable differences
in the criteria that the expert can undoubtedly notice with
a pairwise comparison.</p>
        <p>In addition, when comparing two histograms, the
expert’s perception can be significantly affected by the
order of columns (criteria), therefore, when using this
approach, it is appropriate to think over the ordering or
grouping of criteria based on their features, for example,
by a degree of importance.</p>
        <p>In addition to using this visualization method in the
pairing comparison procedure, it can also be useful in
ranking the criteria, as well as grouping them. This
makes it possible to use this method as a preparatory
stage for the construction and comparison of visual
images of alternatives.</p>
      </sec>
      <sec id="sec-2-2">
        <title>Two-dimensional graphics and threedimensional surfaces</title>
        <p>Other common data visualization methods are
twodimensional graphs and three-dimensional surfaces.
Their application can be effective when there are areas
with a noticeable difference in the values of the criteria,
and due to interpolation, these areas are more
pronounced.</p>
        <p>In the case of two-dimensional graphs, the X axis can
be interpreted as the serial number of the comparison
criterion j (1 ≤  ≤  ) and at Y axis normalized criterion
value –  ′, (1 ≤  ≤ 2) is plotted. If each criterion
corresponds to a numerical value (for example, time),
then it can also be used to determine the X coordinate
(only provided that all these values are pairwise
different) (Fig. 3).</p>
        <p>However, for the application of visualization based
on surfaces in three-dimensional space, a prerequisite is
the presence or ability to display a set of criteria on the
axis X and Y:  ( ) =  ,  ( ) =  , 1 ≤  ≤  . This
mapping can be set based on the physical meaning of the
criteria, or by using the grouping of criteria (for example,
using clustering methods) (Fig. 4).
one alternative prevails over another. Moreover, in the
case of surfaces built in three-dimensional space, a
prerequisite is the availability of tools that allow you to
rotate and zoom in on the surface so that the expert can
choose the most suitable angles for comparison. Thus, it
is important to provide an interactivity property. This
method allows you to compare visual alternatives
entirely based on the subjective perception of the
coverage area of one surface of another, and also with its
help you can determine the various relationships of the
groups of criteria that are used as the coordinates of the
measurements along the abscissa and ordinates.</p>
      </sec>
      <sec id="sec-2-3">
        <title>Probability distribution diagram</title>
        <p>(1 ≤  ≤ 

(Fig. 5).</p>
        <p>Probability distribution diagrams can also be used as
another approach to visualizing data sets describing
alternatives. For these purposes, the interval [0; 1], to
which all normalized values of the criteria ( ′, ) belong,
on M equal intervals in length (for example, by 5 or 10)
depending on the number of criteria. These intervals are
located on the abscissa axis. After that, the number Сi,m
) of hits of the normalized criteria values for
each of the intervals is determined, which then allows
you to determine the corresponding probabilities:   ,
  . , used as values on the ordinate axis when plotting
=</p>
        <p>As in previous approaches for this visualization
method, it is also possible to build diagrams on a single
diagram, or separately. The choice of preference in the
comparison of the probability distribution the expert can
give an alternative, which is characterized by the
displacement of the probability distribution to the right
(towards the interval with the highest criteria of values).
This type of chart is conveniently displayed on a single
diagram while applying transparency of the columns so
that the differences are accented (Fig. 5).</p>
        <p>The
effectiveness
of
this
approach
becomes
significantly higher when the number of criteria K is
sufficiently large (for example, hundreds and thousands).
With its help, not only visual images can be obtained, but
also</p>
        <p>quantitative probabilistic characteristics of the
dominance of one alternative over another. Also, this
method allows us to group criteria, but it does not allow
them to be ranked.
˗ radar diagram with a permutation of criteria, taking
into account their alternation (alternately clockwise
are the criteria with larger and smaller values) (Fig.
9).</p>
        <p>When using this method of visualization, the main
emphasis is on the fact that the best alternative occupies
a larger area, and also that the visual image is brighter
due to the use of gradient fills (in the center, the color is
more neutral – green, and on the periphery – more
contrast – red).</p>
        <p>Fig. 6. Visual comparison of two alternatives using a radial
diagram with sector radii proportional to the values of the
criteria
Fig. 7. Visual comparison of two alternatives using a radial
diagram with sector radii proportional roots of values of the
criteria
Fig. 8. Visual comparison of two alternatives using a radar
diagram with a permutation of criteria by grouping large
values side by side
Fig. 9. Visual comparison of two alternatives using a radar
diagram with a permutation of criteria, taking into account
their alternation</p>
        <p>If one of these visualization methods is used, the
expert, when paired, selects the alternative that seems to
him subjectively brighter and larger in area.</p>
        <p>
          The methods described in [
          <xref ref-type="bibr" rid="ref8">8</xref>
          ] represent a more
universal visualization mechanism, because they allow
one to take into account the order and grouping of
criteria, and also for them integral quantitative
characteristics (brightness, area) can be determined.
        </p>
      </sec>
      <sec id="sec-2-4">
        <title>Complex method</title>
        <p>When comparing alternatives by only one visual
image, it is not always possible to choose the preferred
one from them. This is because the comparison is usually
based on the color, shape, area or volume of the visual
objects defining the respective alternatives. At the same
time, different visualization methods have different
advantages and disadvantages, and often some of them
may not be useful in the visual comparison itself, but in
the preparatory stage, the purpose of which is to
determine the order or grouping of criteria (as histograms
and 3D surface), as well as the integral quantitative
characteristics of visual images - brightness, area of
prevalence, statistical characteristics (probability
distribution diagram), etc. And already these
characteristics allow, for example, to set a specific order
of permutation and grouping of criteria during
visualization (radial and radar diagrams).</p>
        <p>Thus, it is advisable to move from the task of
comparing a single visual object to the task of comparing
a group of visual objects that characterize an alternative
or its components. For this purpose, an integrated
approach is proposed, consisting in the sequential
presentation of a series of visual images obtained using
different methods, until the expert makes a choice.</p>
        <p>For this expert to select the first represented whole
visual images. If with their help he cannot determine the
preferred alternative, then by means of visual analysis he
tries to identify groups of general criteria, and also, if
possible, to rank and filter them. Further, the selected
groups can be visualized separately and placed in a table
grid - on the right are the images for the components of
one alternative, and on the left for the other (Fig. 10).</p>
        <p>In such a set of visual images, there is a high
probability that in a number of rows it will be possible to
choose a preference. If for one alternative there are more
such preferences than for another, then you can make a
choice in favor of this alternative.</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>Experiments</title>
      <p>
        Let us analyze the application of the considered
methods on the example of the alternative (solvers)
described in the works [
        <xref ref-type="bibr" rid="ref10 ref9">9, 10</xref>
        ]. As noted in [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], out of five
solvers, two give the best results – rhoCentralFoam and
pisoCentralFoam (rCF и pCF). Given the fact that the
number of comparison criteria for these two alternatives
is quite large, we will use visual images built on different
diagrams. In Fig. 11 is a visual comparison using
histograms.
      </p>
      <p>For most experts, this comparison will not be
unambiguous, because the images are very similar, and
at the same time on both diagrams, there are both areas
with the best values and the worst.</p>
      <p>We will get an approximately similar result when
using a two-dimensional graph, however, constructing a
surface in three-dimensional space can give a more
interesting result. This method visualization is possible
because criteria can be grouped due to the fact that they
were obtained during computational experiments by
varying two parameters – angle β (in range 10-35° with
step 5°) and Mach numbers (in range 2-7 with step 1), as
well as defined for two norms (L1, L2) four parameters
(Ux, Uy, p, ρ). Analyzing this visual image (Fig. 12), one
can notice that the blue color (rCF solver) prevails on the
surface over red (pCF solver), so the expert can choose
this alternative (rCF slover).</p>
      <p>Fig. 13 shows the results of the visualization of
alternatives using a probability distribution diagram. To
build it, we used a partition of the values of the criteria
into 10 intervals. In this diagram (Fig. 13), it can be noted
that the blue color largely prevails in the columns [0.6;
0.7), [0.8; 0.9) and [0.9; 1.0], which correspond to the
probability of falling into intervals with a higher value
(rank). This means that for the normalized values of the
rCF solver, the probability of obtaining a better solution
is higher. Therefore, the choice of an expert, in this case,
will most likely also be made in favor of this alternative
(rCF solver).
Given the possibility of decomposing criteria into
subsets, we consider the use of complex visual images
based on petal families of diagrams. The decomposition
will be carried out based on various parameters (Ux, Uy,
p, ρ) and norms (L1, L2), i.e. for each alternative, we
construct eight diagrams (Fig. 14).</p>
      <p>Analyzing these visual images, it can be noted that
for the first four rows (L1 norm) due to a more uniform
shape and subjectively somewhat larger area of images
the second alternative (pCF solver) looks preferable,
however, for the L2 norm (5-8 rows) first alternative (rCF
solver) is significantly preferable (due to a more uniform
shape and a subjectively larger area of images). If we
assume that these criteria are peer-to-peer, then it will be
difficult for an expert to determine preference.</p>
      <p>However, if these criteria can be ranked (for example,
the criteria of the L1 block are preferable to the criteria of
the L2 block, or vice versa), then it will be easier for an
expert to make a choice because for this, it will suffice to
compare either only the upper images or only the lower
ones.</p>
    </sec>
    <sec id="sec-4">
      <title>4. Conclusion</title>
      <p>The analysis of visualization methods of alternatives
for the pairwise comparison procedure showed that,
depending on the properties of the source data and their
criteria, various approaches can be both effective and not.
Therefore, it is appropriate to attempt to use several
different visualization methods and their combination in
conjunction with a decomposition of the source data. In
this case, it is possible on some methods to see that one
alternative is better than another due to the subjective
perception of the area of predominance, brightness,
smoothness of forms, etc.</p>
      <p>Fig. 14. Visual comparison of two solvers using a series of
radar diagrams</p>
      <p>The greatest effect in the pairwise comparison
procedure can be achieved by visualization of groups of
initial criteria combining with a ranking of
decomposition parameters. Research in this direction can
be quite promising. Those, we can thereby reduce the
dimension of the initial data set, which will also allow us
to apply traditional decision-making methods, and the
comparison of criteria, in this case, can be based on the
considered visualization methods or supplemented by
them.</p>
    </sec>
    <sec id="sec-5">
      <title>Acknowledgments</title>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Castañón-Puga</surname>
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sanchez</surname>
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Aguilar</surname>
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rodríguez-Díaz</surname>
            <given-names>A</given-names>
          </string-name>
          .
          <string-name>
            <surname>Applied</surname>
          </string-name>
          Decision-Making.
          <source>Applications in Computer Sciences and Engineering</source>
          . - 2019. DOI:
          <volume>10</volume>
          .1007/978-3-
          <fpage>030</fpage>
          - 17985-4.
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Alinezhad</surname>
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khalili</surname>
            <given-names>J</given-names>
          </string-name>
          .
          <article-title>New Methods and Applications in Multiple Attribute Decision Making (MADM)</article-title>
          .
          <article-title>- 2019</article-title>
          . DOI:
          <volume>10</volume>
          .1007/978-3-
          <fpage>030</fpage>
          -15009- 9.
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>David</surname>
            <given-names>H.A.</given-names>
          </string-name>
          <article-title>The Method of Paired Comparisons. (2nd edition</article-title>
          , revised). - Oxford University Press, Berlin. - 1988.
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <surname>Saaty</surname>
            ,
            <given-names>T.L.</given-names>
          </string-name>
          ,
          <year>1980</year>
          . Analytic Hierarchy Process: Planning,
          <string-name>
            <given-names>Priority</given-names>
            <surname>Setting</surname>
          </string-name>
          , and
          <string-name>
            <given-names>Resource</given-names>
            <surname>Allocation</surname>
          </string-name>
          . - New York:
          <string-name>
            <surname>McGraw-Hill</surname>
          </string-name>
          .
          <article-title>- 1980.</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Averchenkov</surname>
            <given-names>V.I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Miroshnikov</surname>
            <given-names>V.V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Podvesovskiy</surname>
            <given-names>A.G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Korostelyov</surname>
            <given-names>D.A.</given-names>
          </string-name>
          (
          <year>2014</year>
          )
          <article-title>Fuzzy and Hierarchical Models for Decision Support in Software Systems Implementations</article-title>
          .
          <source>Proceedings of the 11th Joint Conference</source>
          , JCKBSE 2014 (eds. A.G.
          <article-title>Kravets et</article-title>
          . al.),
          <source>Communications in Computer and Information Science</source>
          ,
          <volume>466</volume>
          , Springer International Publishing, pp.
          <fpage>410</fpage>
          -
          <lpage>421</lpage>
          . DOI:
          <volume>10</volume>
          .1007/978-3-
          <fpage>319</fpage>
          - 11854-3_
          <fpage>35</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          <article-title>[6] Practical application of the paired comparison method</article-title>
          . - URL: https://usabilitylab.ru/blog/boevoeprimenenie-metoda
          <string-name>
            <surname>-</surname>
          </string-name>
          parnyix-sravnenij/.
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <surname>Figuera</surname>
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Greco</surname>
            <given-names>S.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Ehrgott</surname>
            <given-names>M.</given-names>
          </string-name>
          (Eds).
          <article-title>Multiple Criteria Decision Analysis: State of the Art Surveys</article-title>
          . - New York: Springer,
          <year>2005</year>
          . - DOI: 10.1007/b100605.
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <surname>Zakharova</surname>
            <given-names>A.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Korostelyov</surname>
            <given-names>D.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fedonin</surname>
            <given-names>O.N.</given-names>
          </string-name>
          <article-title>Visualization Algorithms for Multi-criteria Alternatives Filtering (</article-title>
          <year>2019</year>
          ).
          <source>Scientific Visualization</source>
          <volume>11</volume>
          .4:
          <fpage>66</fpage>
          -
          <lpage>80</lpage>
          , DOI: 10.26583/sv.11.4.
          <fpage>06</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <surname>Bondarev</surname>
            <given-names>A.E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kuvshinnikov</surname>
            <given-names>A.E.</given-names>
          </string-name>
          <article-title>Analysis of the Accuracy of OpenFOAM Solvers for the Problem of Supersonic Flow Around a Cone /</article-title>
          / ICCS 2018, Lecture Notes in Computer Science (LNCS)
          <fpage>10862</fpage>
          . - P.
          <fpage>221</fpage>
          -
          <lpage>230</lpage>
          ,
          <year>2018</year>
          . - DOI:10.1007/978-3-
          <fpage>319</fpage>
          - 93713-7_
          <fpage>18</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <surname>Bondarev</surname>
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kuvshinnikov</surname>
            <given-names>A</given-names>
          </string-name>
          .
          <article-title>Comparative Estimation of QGDFoam Solver Accuracy for Inviscid Flow Around a Cone // IEEE The Proceedings of the 2018 Ivannikov ISPRAS Open Conference (ISPRAS-</article-title>
          <year>2018</year>
          ). - P.
          <fpage>82</fpage>
          -
          <lpage>87</lpage>
          ,
          <year>2018</year>
          . - DOI: 10.1109/ISPRAS.
          <year>2018</year>
          .
          <volume>00019</volume>
          .
          <string-name>
            <surname>Zakharova Alena</surname>
            <given-names>A.</given-names>
          </string-name>
          , Doctor of Engineering Science, Keldysh Institute of Applied Mathematics, Moscow and professor of the department “Computer Science and Software” of Bryansk State Technical University, Bryansk, Russia. Email:
          <article-title>zaa@tu-bryansk.ru. Korostelyov Dmitriy A</article-title>
          ., Candidate of Engineering Science, Keldysh Institute of Applied Mathematics,
          <article-title>Moscow and associate professor of the department “Computer Science</article-title>
          and Software” of Bryansk State Technical University, Bryansk, Russia. E-mail:
          <article-title>nigm85@mail</article-title>
          .ru.
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>