Next Article in Journal
Improved Performance in the Detection of ACO-OFDM Modulated Signals Using Deep Learning Modules
Next Article in Special Issue
Towards Predicting Student’s Dropout in University Courses Using Different Machine Learning Techniques
Previous Article in Journal
Sleep Quality in Older Women: Effects of a Vibration Training Program
Previous Article in Special Issue
Analysis and Prediction of Engineering Student Behavior and Their Relation to Academic Performance Using Data Analytics Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Table Organization Optimization in Schools for Preserving the Social Distance during the COVID-19 Pandemic

Department of Mechanical, Computer and Aerospace Engineering, Universidad de León, 24071 León, Spain
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2020, 10(23), 8392; https://doi.org/10.3390/app10238392
Submission received: 8 November 2020 / Revised: 20 November 2020 / Accepted: 23 November 2020 / Published: 25 November 2020
(This article belongs to the Special Issue Data Analytics and Machine Learning in Education)

Abstract

:

Featured Application

The obtainment of a methodology for maximizing the social distancing by increasing the distance among the school desks in the classrooms during the coronavirus pandemic through a Genetic Algorithm optimization.

Abstract

The COVID-19 pandemic has supposed a challenge for education. The school closures during the initial coronavirus outbreak for reducing the infections have promoted negative effects on children, such as the interruption of their normal social relationships or their necessary physical activity. Thus, most of the countries worldwide have considered as a priority the reopening of schools but imposing some rules for keeping safe places for the school lessons such as social distancing, wearing facemasks, hydroalcoholic gels or reducing the capacity in the indoor rooms. In Spain, the government has fixed a minimum distance of 1.5 m among the students’ desks for preserving the social distancing and schools have followed orthogonal and triangular mesh patterns for achieving valid table dispositions that meet the requirements. However, these patterns may not attain the best results for maximizing the distances among the tables. Therefore, in this paper, we introduce for the first time in the authors’ best knowledge a Genetic Algorithm (GA) for optimizing the disposition of the tables at schools during the coronavirus pandemic. We apply this GA in two real-application scenarios in which we find table dispositions that increase the distances among the tables by 19.33% and 10%, respectively, with regards to regular government patterns in these classrooms, thus fulfilling the main objectives of the paper.

1. Introduction

COVID-19 was declared a pandemic public health menace on 11 March 2020, by the World Health Organization (WHO) [1]. This supposed the assumption of this severe acute respiratory syndrome (SARS-CoV-2) as a challenge for humanity to deal with a virus with the potential to condition our normal coexistence. In this pandemic, more than 190 countries have been affected by this outbreak and more than 1.23 M deaths have been cumulated globally according to the John Hopkins University as of 5 November 2020 [2].
People with COVID-19 commonly show fever, cough, musculoskeletal symptoms, gastrointestinal symptoms, dyspnea or anosmia/dysgeusia even causing, in the most severe infections, death [3,4,5]. These symptoms can even persist months after overcoming the infection and present novel evidence such as hair loss or cutaneous spots [6]. These symptomatic patients have been treated in hospitals or isolated at home with their close contacts in order to control the propagation of COVID-19 [7]. However, the challenge arises in the asymptomatic coronavirus patients not diagnosed, which propagate the virus without restricting their mobility [8]. Even the symptomatic patients have an initial asymptomatic phase since alveolar macrophages, which are likely the first immune cells to encounter SARS-CoV-2 during the infection, are incapable of sensing the virus in the first stages [9].
As a consequence, countries around the world have opted for severe lockdowns, frontier restrictions or contact tracing for controlling the virus propagation [10] but causing severe impacts on the economy [11] or other different health problems such as mental health or sleep disturbances [12,13].
In this context, education has been severely affected by the effects of the lockdowns that have forced new educational models in which online learning has taken special relevance [14]. These efforts have partially mitigated the school closure and home confinement. However, especially for young children, traditional evidence of studies performed during school holidays have shown that children during these periods are physically less active, modify their sleep patterns, increment their screen time or follow less favorable diets [15]. These effects have been even more pronounced during this outbreak due to the impossibility of the children of socializing with their classmates or playing outside their homes [16].
Therefore, the reopening of schools has been a priority for many governments worldwide [17]. Furthermore, schools are called to remain open during the second national lockdowns in some countries such as France, the United Kingdom, Germany and Italy [18]. However, there are some studies that argue that a strict control on precautionary measures must be done in order to control possible COVID-19 outbreaks at school [19]. Therefore, several restrictions and rules have been imposed for preserving the safe return to school such as the imposition of social distancing, wearing facemasks, reducing the students in the classrooms or hydroalcoholic gels for hand cleaning in every classroom.
In Spain, schools reopened in September facing all these measures described and fixing a social distancing of 1.5 m for reducing the exposition of the children to the virus in the school centers. This led to reorganizing the student desks for guaranteeing the social distancing inside the classrooms.
However, this is a complex problem to be addressed, which forces the reduction of the number of students in traditional classrooms. Many schools tried to find the most appropriate distribution of the tables for fulfilling the government social distancing requirements but they found problems in particularly irregular classrooms where regular patterns in the table disposition (i.e., rectangular or zig-zag grids) do not reach the best achievable results. Even, the beneficial effect of the social distancing for reducing the contagion probability in the COVID-19 pandemic [20] recommends the finding of the table organization that maximizes the distance among the student desks even in classrooms in which the attainment of a valid table disposition can be more easily found.
However, the finding of the optimal table disposition is a combinatorial problem, which is similar to the technological Node Location Problem (NLP) [21], which has been assigned as NP-Hard [22,23]. Therefore, a heuristic solution to this problem is recommended. Simulated annealing [24], the firefly algorithm [25] or the elephant herding optimization [26] have been traditionally used for the NLP even though Genetic Algorithms (GA) [27,28,29] and memetic algorithms [30] are the most recommended for this problem for their trade-off among diversification and intensification of the space of solutions, which is essential for finding optimal results.
As a consequence, in this paper, we propose for the first time, in the authors’ best knowledge a GA optimization for the Table Location Problem (TLP) for finding optimal table dispositions at school that maximizes the distance among the student desks for increasing the social distancing in the classrooms thus reducing the children exposure to COVID-19 in their daily lessons.
The remainder of the paper is organized as follows: we analyze the TLP and its similarities to the NLP together with a complexity analysis in Section 2, the problem definition and the real scenarios in which we applied our algorithm are described in Section 3, the GA for the TLP is introduced in Section 4, the results achieved are presented in Section 5, Section 6 presents the discussion for the TLP and we conclude the manuscript in Section 7.

2. Analysis and Complexity Studies on the Table Location Problem

The TLP entails the definition of the two-dimensional Cartesian coordinates for the location of each student in the classroom. This supposes the associated table distribution in the plane for ensuring the follow-up of the classroom keeping the necessary social distancing by reducing the probability for the students of being infected from the coronavirus during the pandemic.
Therefore, the TLP is a combinatorial problem in which the number of possible table distributions (P) is factorial [30]:
P ( Table Distributions ) = i = 0 n s 1 ( n P L T i )
where n P L T is the number of possible locations where a table can be located in the optimization process and n s the number of students.
This analysis shows the dimensions of the space of solutions in which the size is incremented when considering a larger number of students and when a larger number of possible table locations is defined.
Since the number of students is commonly predetermined, an adequate selection of the n P L T must be performed in order to achieve optimal results. This number should be, on the one hand, high enough for granting sufficient resolution in the table location and on the other hand, reduced for performing time-effective optimizations [29]. In this paper, we have proved different configurations selecting the final hyperparameter based on the achievement of slight modifications in the principal statistical variables when reducing the spatial resolution.
Each one of these considerations makes the TLP factorial in complexity [31] and very similar to the technological NLP, which has been assigned as NP-Hard [22,23]. As a consequence, the finding of an optimal solution for the TLP recommends a heuristic approach similarly to the NLP and can also be categorized as NP-Hard.
Many different metaheuristics have been used for the NLP, which could be applied to the TLP such as simulated annealing [24], the firefly algorithm [25], the elephant herding optimization [26], the dolphin swarm [32], the bat algorithm [33], the grey wolf optimization [34], the bacterial foraging algorithm [35] or diversified local search [36]. However, GA have been the most extended in the literature [27,28,29,37] since they provide an optimal balance between diversification and intensification of the space of solutions for the NLP [38]. In addition, their flexibility and adaptation to any problem makes the GA the best candidate to address the TLP.

3. Problem Definition and Scenario of Application

A mathematical characterization of this problem is required for correctly defining the methods followed for its solution. In this section, we define the problem and the characteristics of the real-world scenarios in which we perform the table localization optimizations.

3.1. Problem Definition

Let t i = x i , y i be the spatial coordinates of the table i considered during the optimization process, n s the number of students in the classroom and consequently the number of tables to be located, T the set including every possible combination of tables in the classroom, T i a subset containing a possible disposition of the tables in the classroom, T j the subset containing every possible combination of T except T i , f T i ( t i ) the evaluation of the quality of the location of the table i belonging to the set T i , the TLP is defined as finding the optimal T i fulfilling the following relation:
T i = t 1 , , t n s T : i = 1 n s f T i ( t i ) n s max i = 1 n s f T j ( t i ) n s
Consequently, the TLP is defined as the finding of the optimal table disposition in the classroom for maximizing a fitness function of which the main purpose is the maximization of the distances among tables but preserving the legislation of each country, which restricts the table disposition to having a minimum required distance among every table.

3.2. Scenario of Application

The TLP is a real-application problem that has been addressed in collaboration with the Marist Brothers School San José in the city of León, Spain. The government of the autonomous country of Castilla and León has promoted the disposition of the students in two kinds of grids for addressing the TLP (i.e., a rectangular and a zig-zag grid) in order to satisfy the 1.5 m separation defined by the Spanish legislation. However, this kind of table disposition is suboptimal and the regular patterns do not necessarily obtain the best results in this type of combinatorial problem [39] since better results can be achieved by using metaheuristics for solving this complex problem.
Therefore, we have analyzed two different classrooms of this school (i.e., Class A and Class B) with different characteristics, which forces different optimization goals. The first class consists of 16 tables while the second scenario needs to allocate 21 students, although the first scenario is smaller, its student density is considerably lower than the second class. Therefore, the Class A scenario is less restrictive than the second classroom in respect to the 1.5 minimum separation protocol.
In order to model these classes and the table allocation area, we have defined three different regions for each class. The first area, the bigger one, represents the class limits, both classes were approximately rectangular; thus, both scenarios are considered to be of such shape.
The second area is defined as the Table Location Environment (TLE) and it is the region within the class limits where the tables can be positioned. For this optimization, the positions of the tables are referenced from the students coordinates; thus, each point of the TLE is a possible location for a student. Therefore, a series of table measurements, shown in Table 1, have been considered in order to properly define the TLE region.
Moreover, although major obstacles of the classrooms were considered into the TLE limit definition, some prohibited areas remained inside the TLE region, thus the codification of these areas is required. The third region in the scenario is the prohibited areas or obstacles, which are defined as areas where a student cannot be allocated. In the simulations presented in this paper, we have displayed the limits of these obstacles with respect to the student’s table for visual clarification. Thus, in the following figures, neither the student nor the table can be positioned inside these obstacle regions.
As for the original student distribution, the first class, shown in Figure 1a, displays an orthogonal mesh distribution of tables, with this scenario being less restrictive in terms of student allocation. On the other hand, Class B, shown in Figure 1b, requires a greater density of students, thus utilizing a triangular mesh distribution, which generally produces more efficient results in the table optimization.
Figure 2a,b shows the classroom models used in this study.

4. Genetic Algorithm Optimization for the TLP

Genetic Algorithms have proven to be among other metaheuristic techniques an excellent trade-off between diversification and intensification of the space of solutions in the optimization procedure. GA were initially proposed by Holland [40] and refined by Goldberg [41] afterward. These algorithms are based on the theory of evolution and rely on the fitness adaptation of a population of individuals to the problem specific scenario.
These individuals contain the problem’s distinct variables, from which the solution is dependent of. Through the optimization, our population is exposed to a certain pressure selection, so that the most adapted individuals prevail and pass out their remarkable genes to the following generation of individuals. This cycle of evolution, carried out through a sufficient amount of iterations, may achieve a state where an individual or a group of individuals contains an optimal enough solution to the problem, thus concluding the optimization.
In the following section we discuss the codification and implementation of the GA proposed for the TLP.

4.1. Codification of the Individuals

The population of a GA must be coded in a way so that it contains all the essential information for the optimization procedure. For this particular case, the main variable of the optimization is the spatial distribution of the tables among the specified region (i.e., the coordinates of each table). Therefore, each individual of the population contains a particular table distribution, carrying the coordinates of each table along the TLE.
Figure 3 shows the codification scheme followed in this paper, executed in the Python programming language, where we create a population list composed of n individuals. Each individual represents a distinct table distribution, carrying i lists for coordinates of each specified table. Every table array carries the two spatial coordinates required for the table positioning, both coded into binary arrays.
The binary codification is particularly useful when undergoing the genetic operators, such as crossover or mutation. Moreover, the binary structure grants a superior degree of flexibility when facing irregular scenarios, through the application of a binary scaling [29].

4.2. Evaluation of the Individuals

Once the population is coded and generated, the next step in the GA structure is to determine the value of each individual, through a specially designed fitness function.
For our particular case of study, the objective of the optimization is the distributions of the student’s tables in a way so that a minimum separation of 1.5 m is obtained, thus fulfilling the Spanish legislation. However, it is desirable in multiple ways for the distance between each pair of tables to be maximized, thus reducing contagion probability and increasing student and professor mobility through the class. Hence, the fitness function proposed must guarantee the minimum separation among every table while also seeking to enhance even beyond the 1.5 m separation distance.
Although it is possible to compute both elements of the optimization simultaneously (i.e., 1.5 minimum separation and mean distance maximization), we have concluded that it is more fruitful to implement a two-step fitness function evaluation that foremost aims to guarantee the safety distance, being this an imperative requisite for the optimization. Once the 1.5 has been obtained, the second phase of the fitness function takes place, seeking to optimize the mean distance between tables. This is similar to the process of predefining an initial population of the GA not randomized containing optimal potential solutions in particularly difficult evolutionary optimization scenarios [42].
For the first phase of the evaluation, we have considered the severity of the 1.5 m infraction into the evaluation of the individuals. In GA optimizations, it is crucial to preserve the convergence through the generations to the final solution. By introducing some progression into the evaluation, even though the individuals may not have achieved the minimum separation, we assure a rewarding mechanism for those individuals that represent a more scattered distribution.
Therefore, the following fitness function equations have been implemented for the first phase of the evaluation:
f f Step 1 = 1.5 l Pen d Pen o b s
Pen d = i = 0 L 1 j = 0 L 1 x i x j 2 + y i y j 2 ; i , j : x i x j 2 + y i y j 2 < 1.5 ; i j
where Pen d is the sum of all the distances between pairs of tables that are under 1.5 m; l is the number of distances that satisfy this statement; L is the number of tables and the number of distances previously measured that are under 1.5 m; and Pen o b s is a penalization applied to those individuals in which some of the tables have been located in a banned region, where Pen o b s values κ if a table is positioned in a prohibited region and Pen o b s values 0; otherwise, κ is a stated hyperparameter.
This fitness function evaluation may vary between [ 1.5 l i = 0 L 1 ( L 1 ) , 0 ] , being the first value obtained when every table is positioned in the same coordinates, and 0 if all the tables maintain a minimum separation of 1.5 m.
This progressing evaluation grants a steadier convergence of the population to a state where most individuals respect the 1.5 m separation.
Once this evaluation has taken place, only if a certain individual or group of individuals have achieved a value of 0, the evaluation proceeds to the following phase.
In this second phase we seek to optimize the mean distance between tables; thus, the following fitness function is proposed:
f f Step 2 = μ λ σ
μ = 1 2 L i = 0 L 1 j = 0 L 1 x i x j 2 + y i y j 2 ; i , j : x i x j 2 + y i y j 2 ρ ; i j
σ = 1 2 L i = 0 L 1 j = 0 L 1 x i x j 2 + y i y j 2 μ 2 ; i , j : x i x j 2 + y i y j 2 ρ ; i j
where μ is the mean distance between pairs of tables; σ is the standard deviation of the distribution of the distances; λ is a weight hyperparameter; and ρ is the hyperparameter that determines whether or not a pair of tables is considered to be close to each other.
In this second phase of the fitness evaluation, we measure both the mean distance and the standard deviation of the table distancing. The introduction of the standard deviation into the optimization aims to obtain a more uniform table distribution. The uniformity of the table distancing plays a vital role in the contagion probability, being undesirable the scenario where some tables are considerably closer than the mean distance. Moreover, we have encountered experimentally that this scenario (i.e., where a reduced number of tables are rather separated from the rest while the rest are substantially close to each other) to be particularly common in the different optimization stages. Therefore, the implementation of the standard deviation in our fitness evaluation addresses this phenomena by penalizing these imbalanced distributions. Figure 4 shows the two-step fitness evaluation procedure proposed for the TLP.
Furthermore, the main goal in this paper is to obtain a table distribution that minimizes the COVID-19 contagion probability, thus, maximizing the separation for every pair of tables is crucial. The implementation of the mean distance between tables as the primary value estimator results in the introduction of the selection pressure required for the GA convergence to the desired solution.

4.3. Selection and Elitism

The selection operator aims to arrange the individuals of the population in a way that enhances the optimization performance. This criterion is based on the fitness value given by the fitness function of each individual in the previous step.
However, we can encounter multiple selection methodologies throughout the literature. These methodologies differ in the pressure selection that they introduce into the optimization. GA optimization is mainly driven by two core aspects, the intensification and the diversification of the solution [43].
The diversification phenomena introduce entropy into the optimization process, this randomness slows the convergence to the final solution, allowing a greater exploration of the solution environment, enhancing the quality of the solutions obtained.
On the other hand, the intensification boosts the convergence to a solution, thus directing the genetic evolution to the optimal path in a reduced space of solutions [44].
The balance between these two factors is key in any GA optimization. An excessive approach to intensification may conclude in premature convergence of the problem into a local maximum due to the lack of exploration, while a heavy intensification focus may compromise the convergence to any solution.
Moreover, this balance depends not only on the problem studied but also may differ from different initial conditions or scenarios of applications. Therefore, we must study the performance of different genetic operators for each particular case of study [45].
Thence, we will analyze the performance of the selection techniques of Tournament 2 (T2), Tournament 3 (T3) and Roulette (R) selection, being these methodologies among the most expanded selection techniques throughout the literature [46].
Furthermore, in addition to the selection methodologies, we can introduce selection pressure into the optimization through the use of elitism. This technique aims to preserve the better adapted individuals throughout the generations, seeking to influence the optimization convergence into the optimal path.
Therefore, through elitism we preserve a certain percentage of individuals along with omitting those with a lesser fitness value [43]. The percentage used is another hyperparameter that needs to be adjusted for each particular application. An excessive value may incur in a premature convergence while an insufficient percentage may not achieve the desired results.

4.4. Crossover and Mutation

The crossover operator seeks to create the new generation of individuals, based on the genetic characteristics of their predecessors. Being the pairing arrangement decided in the selection operator, by combining parents with different genetic properties, it is possible for their offspring to better both of their parents.
However, there coexist multiple crossover techniques, and likewise to the selection operators, their performance depends on the characteristics of the problem’s scenario. Therefore, we must study the behavior of multiple crossover methodologies in search of the most appropriate for every particular problem [47].
Therefore, we analyze the implementation of the Single-Point crossover (SP), MultiPoint crossover for 2 and 3 cross-points (MP2 and MP3, respectively) and Uniform crossover (U), being these techniques among the most expanded methodologies throughout the literature [48].
Moreover, it is possible to introduce a higher degree of entropy into the optimization through the application of mutation. In this operator, we randomly modify certain genes of a specific number of individuals. The purpose of this alteration is to perform a further exploration of the solution environment, since it is possible for a random mutation to induce a rather superior performance of the individual, influencing the optimization direction [43,49].
This operator plays a key role in the GA optimization, and although its implementation might appear less favorable, its adequate application allows the achievement of a greater solution in the whole GA optimization.
However, it is vital to balance the degree of entropy generated, being an excessive amount of mutation detrimental for the convergence; thus, we must study the appropriate amount for the desired scenario of application [43,50].

4.5. Stop Criteria

Once the new generation is created, the iteration of the GA ends, giving way to the successor generation to be evaluated in the following iteration. This cycle of evaluation, selection and crossover continues throughout multiple iterations, until the stop criteria of the GA is fulfilled.
Heuristic methodologies like GA are frequently applied in optimization problems where the optimal solution cannot be easily obtained. Furthermore, in these problems it is common that any solution obtained cannot be verified to be optimal.
For this particular problem, although we can evaluate a certain distribution by measuring its mean distance and verifying the 1.5 m separation accomplishment, we are unable to define if a certain distribution of tables is optimal. Therefore, it is necessary to define a stop criterion in the GA optimization.
In this paper, a double stop condition based on two logic parameters is proposed. Firstly, the GA optimization shall stop if a certain percentage of the population is identical, thus concluding the optimization when the algorithm has converged to a certain solution. In this state, any deviation from this solution presents a lower fitness value than the previous value; however, this phenomena is not limited to the optimal solution, being possible for this solution to be a local maximum [51].
Secondly, the GA optimization shall also stop if a certain number of iterations have been completed. The addition of this parameter is crucial in any GA optimization, since it is possible that the GA does not achieve a convergence to any solution, thus entering the GA into an infinite loop.

5. Results

In this section, we present the results of the previously detailed GA in both classroom scenarios. All algorithms were coded and executed in the Python software environment, performing every test with an Intel(R) i7 2.4 GHz of CPU and 16 GB f RAM.
Previously, we have studied the significant impact of different genetic operators into the GA global performance. Their effect on the balance between intensification and diversification of the solution is key in order to obtain an adequate solution through the GA optimization.
Therefore, in order to obtain the most appropriate genetic operators for each particular scenario, firstly, we must study the performance of each individual combination of genetic operators through multiple simulations. All simulations were executed with the following parameters, as shown in Table 2.
Hence, in search of the optimal combination of genetic operators, we have analyzed the performance of each possible individual combination among these functions, resulting in the comparison shown in Table 3.
Table 3 shows the difference in the optimization performance of the GA depending on the combination of genetic operators selected and the scenario of application.
The first classroom presents a great opportunity for table location optimization; however, the results of the optimization differs from each combination of selection and crossover techniques utilized.
The second scenario is far more restrictive than the previous one, and in some combinations, the GA did not achieve a valid solution, proving the importance of testing the performance of multiple operators.
For both scenarios, the roulette selection and the uniform crossover showed inferior performance than the other studied techniques. Generally, R. selection and U. crossover are particularly heavy focused techniques on intensification and diversification, respectively, thus their introduction to this problem may result inadequate.
On the other hand, from the other combinations, tournament-2 with three-point crossover (T2-MP3) and two-point crossover (T2-MP2) showed satisfying results for the first and second scenario analyzed, respectively. These combinations of selection and crossover methodologies exhibiting an adequate balance between diversification and intensification in the maximum and mean fitness values obtained throughout the simulations.
Furthermore, once we selected both methodologies, we evaluated the most advantageous elitism and mutation values for each scenario of application, as shown in Table 4.
The resulting configurations of the GA proposed in Table 4, executed for both scenarios obtaining the following table distribution, as shown in Figure 5. The evolution of the algorithms and convergence to the final solution is also shown in Figure 6.
These resulting distributions obtained by the GA in both scenarios thrive at optimizing the available space, achieving a mean separation, shown in Table 5, unreachable for regular meshes.
Furthermore, due to the introduction of the standard deviation into the fitness evaluation, the resulting distributions, although presenting irregular patterns, exhibit some degree of uniformity. This fact proves essential when minimizing Covid-19 contagion in our scenarios.
Results in Table 5 show an improvement up to 19.33% in the mean distance between tables from the GA optimization. This statement proves the viability of these algorithms for any scenario of application, also achieving a 10% increase in mean distance in the second, most restrictive scenario.
Moreover, the obtained distribution for the first scenario was later implemented into the original classroom studied, proving the viability of this methodology, resulting in the distribution shown in Figure 7.
Therefore, the GA proposed in this paper successfully enhances the table distribution for multiple scenarios with respect to the Covid-19 contagion probability, thus fulfilling the main objective of this research.

6. Discussion

This paper presents a technological solution for the TLP in the schools during the COVID-19 pandemic. The distance achieved among the students’ desks allows the reduction of the probability of contagion of this emerging infection, which collaborates with other rules such as wearing facemasks or using hydroalcoholic gels for creating safe places for the children to follow their daily lessons.
Our approach addresses a novel problem, which is firstly introduced in this paper. The location of the tables has been proved in the manuscript to be a combinatorial NP-Hard problem, which recommends a heuristic solution to find acceptable results. The huge dimension of the space of solutions dependent on the resolution of the TLE and the number of tables to be located introduces difficulties during the optimization process, which requires the definition of a two-step optimization procedure in which we first ensure a combination of the tables in a space that meets the government 1.5 m minimum separation among all the tables and we later expand this distance to minimize the contagion risks in the classrooms.
We have implemented a GA optimization for this purpose due to the flexibility of adaptation of this metaheuristic to different similar technological problems [27,28,29,37,52,53,54] and the trade-off achieved between diversification and intensification of the space of solutions. The results obtained in the manuscript have demonstrated the suitability of the application of metaheuristics to solve this kind of problem that emerged during the pandemic, improving the government solutions proposed.
In our future works, we will extend the analysis of the TLP to other different metaheuristics such as simulated annealing, diversified local search or the combination of the GA with a local search procedure to explore improvements in the optimization results. Furthermore, we will consider novel optimization scenarios such as the school canteen, restaurants, or pubs in which novel challenges for the optimization can arise. Such as multi-objective optimization for different criteria leading to algorithms such as NSGA-II, NSGA-III or MOEA.

7. Conclusions

The COVID-19 pandemic has supposed a challenge for humanity to deal with a virus that has changed the normal coexistence. As a consequence, many restrictions have been imposed globally for reducing the probability of contagion of a virus that can even cause, in the most severe cases, death. Therefore, social distancing, wearing facemasks, hydroalcoholic gels for hand cleaning or reducing the capacity allowed in indoor spaces have been some of the rules adopted for containing the virus propagation.
Education has been one of the most affected sectors. Most of the countries decided to move to online learning during the first lockdowns promoted for facing the emergency of the initial coronavirus outbreak. This has complicated the normal learning of children, their social relationships and their physical health activity. Consequently, most of the countries have considered the reopening of the schools as a priority, even keeping the schools opened during the second coronavirus outbreak that it is nowadays facing Europe. This has promoted the implementation of restrictions and rules at schools for reducing the infectious potential of the virus.
One of them has consisted of smart dispositions of the tables in the classrooms for maintaining the social distancing. In Spain, the government has fixed a minimum of 1.5 m among the students’ desks, which has led to reducing the number of children in each classroom and designing regular patterns in the disposition of the tables such as orthogonal or triangular mesh configurations.
However, the problem of the disposition of the tables is NP-Hard and a metaheuristic solution is recommended for obtaining improved results through irregular table disposition patterns that maximize the distances among the tables thus minimizing the probability of getting the coronavirus at schools.
In this paper, we introduce for the first time to the authors’ best knowledge a Genetic Algorithm optimization for the Table Location Problem for addressing the table disposition during the COVID-19 pandemic. We analyze the definition and complexity of the problem and we propose a methodology for its resolution. This methodology is applied in two different real-application scenarios (i.e., Class A and Class B) in which we prove the suitability of the optimization of the table disposition for obtaining improved results. Results show that an increase of the mean distance among tables of 19.33% in Class A and 10% in Class B can be attained following the proposed methodology thus fulfilling the main objectives of this paper.

Author Contributions

Conceptualization, R.F.-G. and J.D.-G.; formal analysis, H.P.; investigation, R.F.-G., J.D.-G. and P.V.; methodology, R.F.-G., J.D.-G. and R.Á.; project administration, H.P.; resources, P.V. and R.Á.; software, R.F.-G., J.D.-G. and R.Á.; supervision, J.D.-G. and H.P.; validation, R.Á. and H.P.; visualization, R.F.-G. and J.D.-G.; writing—original draft, R.F.-G. and J.D.-G.; writing—review and editing, P.V., R.Á. and H.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministry of Science and Innovation grant number PID2019-108277GB-C21.

Acknowledgments

The authors want to acknowledge the Marist Brothers School San José in the city of León for their collaboration in the research presented in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
GAGenetic Algorithm
MP2Two-Points MultiPoint crossover
MP3Three-Points MultiPoint crossover
NLPNode Location Problem
RRoulette selection
SARS-CoV-2Severe Acute Respiratory Syndrome
TLETable Location Environment
TLPTable Location Problem
T2Tournament-2 selection
T3Tournament-3 selection
UUniform crossover
WHOWorld Health Organization

References

  1. McKay, B.; Calfas, J.; Ansari, T. Coronavirus declared pandemic by World Health Organization. Wall Street J. 2020, 11, 59. [Google Scholar]
  2. Johns Hopkins University & Medicine. Covid Global Map—COVID-19 Dashboard by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU). Available online: https://coronavirus.jhu.edu/map.html (accessed on 5 November 2020).
  3. Carfì, A.; Bernabei, R.; Landi, F.; The Gemelli against COVID-19 Post-Acute Care Study Group. Persistent symptoms in patients after acute COVID-19. Jama 2020, 324, 603–605. [Google Scholar]
  4. Klopfenstein, T.; Toko, L.; Royer, P.Y.; Lepiller, Q.; Gendrin, V.; Zayet, S. Features of anosmia in COVID-19. Med. Mal. Infect. 2020, 50, 436–439. [Google Scholar] [CrossRef] [PubMed]
  5. Ruan, Q.; Yang, K.; Wang, W.; Jiang, L.; Song, J. Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China. Intensive Care Med. 2020, 46, 846–848. [Google Scholar] [CrossRef] [Green Version]
  6. Suchonwanit, P.; Leerunyakul, K.; Kositkuljorn, C. Cutaneous manifestations in COVID-19: Lessons learned from current evidence. J. Am. Acad. Dermatol. 2020. [Google Scholar] [CrossRef]
  7. Kretzschmar, M.E.; Rozhnova, G.; Bootsma, M.C.; van Boven, M.; van de Wijgert, J.H.; Bonten, M.J. Impact of delays on effectiveness of contact tracing strategies for COVID-19: A modelling study. Lancet Public Health 2020, 5, e452–e459. [Google Scholar] [CrossRef]
  8. Gandhi, M.; Yokoe, D.S.; Havlir, D.V. Asymptomatic transmission, the Achilles’ heel of current strategies to control COVID-19. N. Engl. J. Med. 2020, 382, 2158–2160. [Google Scholar] [CrossRef]
  9. Dalskov, L.; Møhlenberg, M.; Thyrsted, J.; Blay-Cadanet, J.; Poulsen, E.T.; Folkersen, B.H.; Skaarup, S.H.; Olagnier, D.; Reinert, L.; Enghild, J.J.; et al. SARS-CoV-2 evades immune detection in alveolar macrophages. EMBO Rep. 2020, e51252. [Google Scholar] [CrossRef]
  10. Hellewell, J.; Abbott, S.; Gimma, A.; Bosse, N.I.; Jarvis, C.I.; Russell, T.W.; Munday, J.D.; Kucharski, A.J.; Edmunds, W.J.; Sun, F.; et al. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. Lancet Glob. Health 2020, 8, e488–e496. [Google Scholar] [CrossRef] [Green Version]
  11. Nicola, M.; Alsafi, Z.; Sohrabi, C.; Kerwan, A.; Al-Jabir, A.; Iosifidis, C.; Agha, M.; Agha, R. The socio-economic implications of the coronavirus pandemic (COVID-19): A review. Int. J. Surg. 2020, 78, 185. [Google Scholar] [CrossRef]
  12. Gualano, M.R.; Lo Moro, G.; Voglino, G.; Bert, F.; Siliquini, R. Effects of Covid-19 lockdown on mental health and sleep disturbances in Italy. Int. J. Environ. Res. Public Health 2020, 17, 4779. [Google Scholar] [CrossRef] [PubMed]
  13. Marelli, S.; Castelnuovo, A.; Somma, A.; Castronovo, V.; Mombelli, S.; Bottoni, D.; Leitner, C.; Fossati, A.; Ferini-Strambi, L. Impact of COVID-19 lockdown on sleep quality in university students and administration staff. J. Neurol. 2020. [Google Scholar] [CrossRef] [PubMed]
  14. Viner, R.M.; Russell, S.J.; Croker, H.; Packer, J.; Ward, J.; Stansfield, C.; Mytton, O.; Bonell, C.; Booy, R. School closure and management practices during coronavirus outbreaks including COVID-19: A rapid systematic review. Lancet Child Adolesc. Health 2020, 4, 397–404. [Google Scholar] [CrossRef]
  15. Brazendale, K.; Beets, M.W.; Weaver, R.G.; Pate, R.R.; Turner-McGrievy, G.M.; Kaczynski, A.T.; Chandler, J.L.; Bohnert, A.; von Hippel, P.T. Understanding differences between summer vs. school obesogenic behaviors of children: The structured days hypothesis. Int. J. Behav. Nutr. Phys. Act. 2017, 14, 100. [Google Scholar] [CrossRef]
  16. Wang, G.; Zhang, Y.; Zhao, J.; Zhang, J.; Jiang, F. Mitigate the effects of home confinement on children during the COVID-19 outbreak. Lancet 2020, 395, 945–947. [Google Scholar] [CrossRef]
  17. Sheikh, A.; Sheikh, A.; Sheikh, Z.; Dhami, S. Reopening schools after the COVID-19 lockdown. J. Glob. Health 2020, 10. [Google Scholar] [CrossRef]
  18. New York Times. Why Is Europe Keeping Its Schools Open, Despite New Lockdowns? Available online: https://www.nytimes.com/2020/10/29/world/europe/schools-coronavirus-europe-lockdowns.html (accessed on 4 November 2020).
  19. Stein-Zamir, C.; Abramson, N.; Shoob, H.; Libal, E.; Bitan, M.; Cardash, T.; Cayam, R.; Miskin, I. A large COVID-19 outbreak in a high school 10 days after schools’ reopening, Israel, May 2020. Eurosurveillance 2020, 25, 2001352. [Google Scholar] [CrossRef]
  20. Wilder-Smith, A.; Freedman, D.O. Isolation, quarantine, social distancing and community containment: Pivotal role for old-style public health measures in the novel coronavirus (2019-nCoV) outbreak. J. Travel Med. 2020, 27. [Google Scholar] [CrossRef]
  21. Landa-Torres, I.; Gil-Lopez, S.; Salcedo-Sanz, S.; Del Ser, J.; Portilla-Figueras, J.A. A novel grouping harmony search algorithm for the multiple-type access node location problem. Expert Syst. Appl. 2012, 39, 5262–5270. [Google Scholar] [CrossRef]
  22. Tekdas, O.; Isler, V. Sensor placement for triangulation-based localization. IEEE Trans. Autom. Sci. Eng. 2010, 7, 681–685. [Google Scholar] [CrossRef]
  23. Nguyen, N.T.; Liu, B.H. The mobile sensor deployment problem and the target coverage problem in mobile wireless sensor networks are NP-hard. IEEE Syst. J. 2018, 13, 1312–1315. [Google Scholar] [CrossRef]
  24. Wang, X.; Ma, J.J.; Wang, S.; Bi, D.W. Distributed particle swarm optimization and simulated annealing for energy-efficient coverage in wireless sensor networks. Sensors 2007, 7, 628–648. [Google Scholar] [CrossRef] [Green Version]
  25. Tuba, E.; Tuba, M.; Beko, M. Two stage wireless sensor node localization using firefly algorithm. In Smart Trends in Systems, Security and Sustainability; Springer: Berlin/Heidelberg, Germany, 2018; pp. 113–120. [Google Scholar]
  26. Correia, S.D.; Beko, M.; da Silva Cruz, L.A.; Tomic, S. Elephant herding optimization for energy-based localization. Sensors 2018, 18, 2849. [Google Scholar]
  27. Yoon, Y.; Kim, Y.H. An efficient genetic algorithm for maximum coverage deployment in wireless sensor etworks. IEEE Trans. Cybern. 2013, 43, 1473–1483. [Google Scholar] [CrossRef] [PubMed]
  28. Domingo-Perez, F.; Lazaro-Galilea, J.L.; Wieser, A.; Martin-Gorostiza, E.; Salido-Monzu, D.; de la Llana, A. Sensor placement determination for range-difference positioning using evolutionary multi-objective optimization. Expert Syst. Appl. 2016, 47, 95–105. [Google Scholar] [CrossRef]
  29. Díez-González, J.; Álvarez, R.; González-Bárcena, D.; Sánchez-González, L.; Castejón-Limas, M.; Perez, H. Genetic algorithm approach to the 3D node localization in TDOA systems. Sensors 2019, 19, 3880. [Google Scholar] [CrossRef] [Green Version]
  30. Díez-González, J.; Verde, P.; Ferrero-Guillén, R.; Álvarez, R.; Pérez, H. Hybrid Memetic Algorithm for the Node Location Problem in Local Positioning Systems. Sensors 2020, 20, 5475. [Google Scholar] [CrossRef]
  31. Díez-Gonzalez, J.; Álvarez, R.; Perez, H. Optimized Cost-Effective Node Deployments in Asynchronous Time Local Positioning Systems. IEEE Access 2020, 8, 154671–154682. [Google Scholar] [CrossRef]
  32. Kannadasan, K.; Edla, D.R.; Kongara, M.C.; Kuppili, V. M-curves path planning model for mobile anchor node and localization of sensor nodes using dolphin swarm algorithm. Wirel. Netw. 2019, 26, 2769–2783. [Google Scholar] [CrossRef]
  33. Mihoubi, M.; Rahmoun, A.; Lorenz, P. An Enhanced Bat Algorithm for Parallel Localization Based on a Mobile Beacon Sensor in Wireless Sensor Networks. In International Conference on Advanced Communication and Networking; Springer: Berlin/Heidelberg, Germany, 2019; pp. 43–61. [Google Scholar]
  34. Rajakumar, R.; Amudhavel, J.; Dhavachelvan, P.; Vengattaraman, T. GWO-LPWSN: Grey wolf optimization algorithm for node localization problem in wireless sensor networks. J. Comput. Netw. Commun. 2017, 2017, 7348141. [Google Scholar] [CrossRef] [Green Version]
  35. Kulkarni, R.V.; Venayagamoorthy, G.K.; Cheng, M.X. Bio-inspired node localization in wireless sensor networks. In Proceedings of the 2009 IEEE International Conference on Systems, Man and Cybernetics, San Antonio, TX, USA, 11–14 October 2009; pp. 205–210. [Google Scholar]
  36. Laguna, M.; Roa, J.O.; Jiménez, A.R.; Seco, F. Diversified local search for the optimal layout of beacons in an indoor positioning system. Iie Trans. 2009, 41, 247–259. [Google Scholar] [CrossRef]
  37. Kratica, J.; Tošic, D.; Filipović, V.; Ljubić, I. Solving the simple plant location problem by genetic algorithm. RAIRO Oper. Res. Rech. Opér. 2001, 35, 127–142. [Google Scholar] [CrossRef] [Green Version]
  38. Zhang, H.; Liu, Y.; Zhou, J. Balanced-evolution genetic algorithm for combinatorial optimization problems: The general outline and implementation of balanced-evolution strategy based on linear diversity index. Nat. Comput. 2018, 17, 611–639. [Google Scholar] [CrossRef]
  39. He, T.; Huang, C.; Blum, B.M.; Stankovic, J.A.; Abdelzaher, T. Range-free localization schemes for large scale sensor networks. In Proceedings of the 9th Annual International Conference on Mobile Computing and Networking, San Diego, CA, USA, 14–19 September 2003; pp. 81–95. [Google Scholar]
  40. Holland, J.H. Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  41. Goldberg, D.E. Genetic Algorithms; Pearson Education India: New Delhi, India, 2006. [Google Scholar]
  42. Álvarez, R.; Díez-González, J.; Strisciuglio, N.; Perez, H. Multi-Objective Optimization for Asynchronous Positioning Systems Based on a Complete Characterization of Ranging Errors in 3D Complex Environments. IEEE Access 2020, 8, 43046–43056. [Google Scholar] [CrossRef]
  43. Rudolph, G. Convergence analysis of canonical genetic algorithms. IEEE Trans. Neural Netw. 1994, 5, 96–101. [Google Scholar] [CrossRef] [Green Version]
  44. Ahn, C.W.; Ramakrishna, R.S. Elitism-based compact genetic algorithms. IEEE Trans. Evol. Comput. 2003, 7, 367–385. [Google Scholar]
  45. Contreras-Bolton, C.; Parada, V. Automatic combination of operators in a genetic algorithm to solve the traveling salesman problem. PLoS ONE 2015, 10, e0137724. [Google Scholar] [CrossRef] [Green Version]
  46. Shukla, A.; Pandey, H.M.; Mehrotra, D. Comparative review of selection techniques in genetic algorithm. In Proceedings of the 2015 International Conference on Futuristic Trends on Computational Analysis and Knowledge Management (ABLAZE), Noida, India, 25–27 February 2015; pp. 515–519. [Google Scholar]
  47. Ferrero-Guillén, R.; Díez-González, J.; Álvarez, R.; Pérez, H. Analysis of the Genetic Algorithm Operators for the Node Location Problem in Local Positioning Systems. In Hybrid Artificial Intelligent Systems; de la Cal, E.A., Villar Flecha, J.R., Quintián, H., Corchado, E., Eds.; Springer: Cham, Switzerland, 2020; pp. 273–283. [Google Scholar]
  48. Umbarkar, A.J.; Sheth, P.D. Crossover operators in genetic algorithms: A review. ICTACT J. Soft Comput. 2015, 6. [Google Scholar] [CrossRef]
  49. Bäck, T. Optimal mutation rates in genetic search. In Proceedings of the Fifth International Conference on Genetic Algorithms, SanMateo, CA, USA, 17–21 July 1993. [Google Scholar]
  50. Mirjalili, S. Genetic algorithm. In Evolutionary Algorithms and Neural Networks; Springer: Berlin/Heidelberg, Germany, 2019; pp. 43–55. [Google Scholar]
  51. Greenhalgh, D.; Marshall, S. Convergence criteria for genetic algorithms. SIAM J. Comput. 2000, 30, 269–282. [Google Scholar] [CrossRef]
  52. Díez-González, J.; Álvarez, R.; Prieto-Fernández, N.; Perez, H. Local Wireless Sensor Networks Positioning Reliability Under Sensor Failure. Sensors 2020, 20, 1426. [Google Scholar] [CrossRef] [Green Version]
  53. Marinakis, Y.; Marinaki, M. A bilevel genetic algorithm for a real life location routing problem. Int. J. Logist. Res. Appl. 2008, 11, 49–65. [Google Scholar] [CrossRef]
  54. Muniyandi, R.C.; Qamar, F.; Jasim, A.N. Genetic Optimized Location Aided Routing Protocol for VANET Based on Rectangular Estimation of Position. Appl. Sci. 2020, 10, 5759. [Google Scholar] [CrossRef]
Figure 1. Picture of both classrooms from the Marist Brothers School San José, used for the scenario modeling and optimization in this paper.
Figure 1. Picture of both classrooms from the Marist Brothers School San José, used for the scenario modeling and optimization in this paper.
Applsci 10 08392 g001
Figure 2. Models for the two classrooms under study, Class A (a) deployed an orthogonal mesh distribution while Class B (b) displayed a triangular mesh table allocation.
Figure 2. Models for the two classrooms under study, Class A (a) deployed an orthogonal mesh distribution while Class B (b) displayed a triangular mesh table allocation.
Applsci 10 08392 g002
Figure 3. Binary codification of the population in the Genetic Algorithm (GA) proposed for the Table Location Problem (TLP).
Figure 3. Binary codification of the population in the Genetic Algorithm (GA) proposed for the Table Location Problem (TLP).
Applsci 10 08392 g003
Figure 4. Flux diagram of the GA fitness evaluation proposed for the TLP.
Figure 4. Flux diagram of the GA fitness evaluation proposed for the TLP.
Applsci 10 08392 g004
Figure 5. Table distributions obtained for each scenario. The application of each optimization is based on the configuration proposed in Table 4.
Figure 5. Table distributions obtained for each scenario. The application of each optimization is based on the configuration proposed in Table 4.
Applsci 10 08392 g005
Figure 6. Evolution of the GA optimization for both scenarios and convergence to the final solution. Class B requires additional generations to obtain an adequate solution, due to the restrictive limitations of the scenario. The rapid increase in fitness value for both functions when surpassing the 0 axis is due to the change in the fitness evaluation technique.
Figure 6. Evolution of the GA optimization for both scenarios and convergence to the final solution. Class B requires additional generations to obtain an adequate solution, due to the restrictive limitations of the scenario. The rapid increase in fitness value for both functions when surpassing the 0 axis is due to the change in the fitness evaluation technique.
Applsci 10 08392 g006
Figure 7. Implementation of the optimized distribution into the original classroom.
Figure 7. Implementation of the optimized distribution into the original classroom.
Applsci 10 08392 g007
Table 1. List of parameters measured from the studied classrooms, implemented into the scenario modeling and table distribution optimization.
Table 1. List of parameters measured from the studied classrooms, implemented into the scenario modeling and table distribution optimization.
ParametersValues
Table Length69.5 cm
Table Width49.5 cm
Student Space Radius10 cm
Distance from Student to Table20 cm
Table 2. GA hyperparameters used for all simulations, the hyperparameter values were optimized and adjusted experimentally. Elitism and Mutation values are provisional and require further adjustment.
Table 2. GA hyperparameters used for all simulations, the hyperparameter values were optimized and adjusted experimentally. Elitism and Mutation values are provisional and require further adjustment.
GA HyperparametersClass AClass B
Stop criteria800 Generations
30% Population equal
Number of Individuals200
Elitism10%
Mutation10%
ρ 2 m
κ 500 m
Number of Tables1621
TLP Points409616,384
Number of Possible Combinations6.09 × 10 57 3.14 × 10 88
Table 3. Comparison of the fitness value obtained by the multiple selection and crossover genetic operators previously proposed.
Table 3. Comparison of the fitness value obtained by the multiple selection and crossover genetic operators previously proposed.
Crossover OperatorsTournament 2Tournament 3Roulette
MaxMeanMaxMeanMaxMean
Single point—Class A1.6481.5531.6021.2641.5151.169
Two-point—Class A1.6341.2291.5721.2371.5631.158
Three-point—Class A1.6521.571.6461.5791.6381.545
Uniform—Class A1.6221.0891.6690.7391.5060.11
Single point—Class B−0.262−4.761.681−2.312−1.296−.611
Two-point—Class B1.676−0.628−0.876−2.0171.675−1.264
Three-point—Class B1.677−1.046−0.482−3.1021.645−2.04
Uniform—Class B−0.657−2.961−0.866−4.204−1.64−4.59
Table 4. Adjustment of the most rewarding elitism and mutation values for each particular scenario, for the methodologies selected. Values were obtained experimentally.
Table 4. Adjustment of the most rewarding elitism and mutation values for each particular scenario, for the methodologies selected. Values were obtained experimentally.
ScenarioMethodologyElitismMutation
Class AT2-MP315%12.5%
Class BT2-MP27.5%7.5%
Table 5. Comparison of the mean table separation achieved between the GA optimized distribution and the initial mesh models.
Table 5. Comparison of the mean table separation achieved between the GA optimized distribution and the initial mesh models.
ScenarioGA OptimizationOriginal DistributionImprovement
Class A1.79 m1.5 m19.33%
Class B1.65 m1.5 m10%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ferrero-Guillén, R.; Díez-González, J.; Verde, P.; Álvarez, R.; Perez, H. Table Organization Optimization in Schools for Preserving the Social Distance during the COVID-19 Pandemic. Appl. Sci. 2020, 10, 8392. https://doi.org/10.3390/app10238392

AMA Style

Ferrero-Guillén R, Díez-González J, Verde P, Álvarez R, Perez H. Table Organization Optimization in Schools for Preserving the Social Distance during the COVID-19 Pandemic. Applied Sciences. 2020; 10(23):8392. https://doi.org/10.3390/app10238392

Chicago/Turabian Style

Ferrero-Guillén, Rubén, Javier Díez-González, Paula Verde, Rubén Álvarez, and Hilde Perez. 2020. "Table Organization Optimization in Schools for Preserving the Social Distance during the COVID-19 Pandemic" Applied Sciences 10, no. 23: 8392. https://doi.org/10.3390/app10238392

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop