Public Transportation Network Design and Frequency Setting: Pareto Optimality through Alternating-Objective Genetic Algorithms

The transportation network design and frequency setting problem concerns the optimization of transportation systems comprising fleets of vehicles serving a set amount of passengers on a predetermined network (e.g., public transport systems). It has been a persistent focus of the transportation planning community while, its NP-hard nature continues to present obstacles in designing efficient, all-encompassing solutions. In this paper, we present a new approach based on an alternating-objective genetic algorithm that aims to find Pareto optimality between user and operator costs. Extensive computational experiments are performed on Mandl’s benchmark test and prove that the results generated by our algorithm are 5–6% improved in comparison to previously published results for Pareto optimality objectives both in regard to user and operator costs. At the same time, the methods presented are computationally inexpensive and easily run on office equipment, thus minimizing the need for expensive server infrastructure and costs. Additionally, we identify a wide variance in the way that similar computational results are reported and, propose a novel way of reporting benchmark results that facilitates comparisons between methods and enables a taxonomy of heuristic approaches to be created. Thus, this paper aims to provide an efficient, easily applicable method for finding Pareto optimality in transportation networks while highlighting specific limitations of existing research both in regards to the methods used and the way they are communicated.


Introduction
The field of public transport optimization has attracted immense research interest for more than 50 years. The increase in the Earth's population, increasing urbanization and widespread urban development have made the need for efficient urban transportation systems ever more pressing [1,2]. To this day, and for the foreseeable future, cities in every continent are having to build or restructure their services to meet the demands of increased population sizes [3,4]. However, across the world, city center road networks are buckling from congestion issues that lead to environmental and noise pollution, increased road accident rates and ever decreasing quality of life indexes [5][6][7]. This has turned the sights of researchers, competent authorities, and key stakeholders towards improved public transport systems, as they provide a cost-effective avenue of regulating intercity mobility, stimulating sustainable growth, and increasing citizen satisfaction [8,9]. The problem is called Transportation Network Design and Frequency Setting Problem (TRNDFSP) and this paper presents a novel twist on established approaches to create efficient networks.

Evolutionary Algorithms
In computer science, evolutionary computation is a family of algorithms for global optimization inspired by biological evolution, and the subfield of artificial intelligence and soft computing studying these algorithms. Such algorithms include genetic algorithms, genetic programming, evolutionary programming, which can be adjusted to address singleobjective or multi-objective problems. In most cases, as seen previously, the TRNDFSP is handled as a multi-objective problem.
According to recent research in the field, notable optimization of multi-objective problems can be accomplished with a learning-based algorithm that includes learning automation, and later adjusting the evolutionary strategies of the algorithm to adapt to the problem characteristics, according to the feedback information during the optimizing procedure [40]. In another study, a new adaptive island evolutionary algorithm was proposed for the berth scheduling problem, aiming to minimize the total weighted service cost of vessels. The optimality gap of the algorithm is relatively low and the proposed solution algorithm exhibits statistically significant improvements in terms of the objective function values when compared to state-of-the-art metaheuristic algorithms [41]. Furthermore, a research that developed the AnD algorithm (angle-based selection strategy and a shift-based density estimation strategy), showed that it can achieve highly competitive performance when compared to other state-of-the-art many-objective evolutionary algorithms on a variety of benchmark test problems with up to 15 objectives. AnD has a simple structure, few parameters, and no complicated operators [42]. Furthermore, in the field of bioinformatics the use of machine learning-based methodologies show high performance in the distinction of fatal bacterial meningitis from viral meningitis, achieving 100% of sensitivity in detecting the former [43].
As for the transportation sector, a research was conducted to address the vehicle routing problem, by developing an optimization model and solution algorithms which utilizes a mixed-integer linear programming model and employs four metaheuristic algorithms, including the evolutionary algorithm, variable neighborhood search, tabu search, and simulated annealing to solve the model for large-scale problem instances [44].
Lastly, it is important to mention the significance of data classification in the process of optimizing a multi-objective problem. According to research, the use of swan algorithms, and specifically the Salp swam algorithm in standard classification measures over considered benchmark datasets shows supremacy in results as compared with other evolutionary algorithm based multilayer perceptron, such as genetic algorithm (GA), particle swarm optimization (PSO), differential evolution (DE), grey wolf optimization (GWO), etc [45].

Problem Description
The bus network design is formulated as a non-convex, NP-hard optimization problem dealing with the minimization of all resources and costs related to a public transport system with fixed demand [29]. Common problem inputs are the distribution of trip generation within an area, its topology along with a set of objectives and constraints [46]. The objectives that appear more frequently in relevant literature can be summarized as follows [21,47]: (1) User benefit maximization (2) Operator cost minimization (3) Maximization of total system welfare (4) Maximization of network capacity (5) Energy conservation (6) Parameter-focused optimization From all the above, recent literature has aimed at improving overall system welfare, as it encompasses both user and operator costs [26,48].
Regarding the problem formulation, a plethora of authors [18,22,23,25,27,31,49] have described the problem in, each time, slightly different language. The formulation that is presented in the following paragraphs aims to unify the problem description under a common, mathematical formulation. It draws from graph theory and problem descriptions of previ-ous papers to present a complete formulation of the transportation network design and frequency setting problem.
The road network that is studied can be modelled as an undirected planar graph G (flow network) consisting of a set of vertices, V, and a set of edges, E, along with their respective flows ( f ). This way, we can mathematically describe the OD matrix, containing the passenger flows from each vertex (i.e., bus stop, centroid) of the matrix to any other vertex along with the way that these are linked; with null values between vertices that are not directly linked [50,51]. Additionally, the necessary time matrix for the design of the transportation network is included in the mathematical description: where E is a set of edges e i : where e i , ith edge of the edge set; u i s , start vertex of the i-th edge; u i s , end vertex of the i-th edge; w i , weight of the edge, that represented the travel distance; l i , the length of the edge; p i , edge passengers travelling from u i s to u i e ; id i , identifies the edge among other edges; u, vertex notation.
Additionally, we consider the family of sets of routes between the vertices that fulfill a set of constraints described in a following section. These are denoted as follows: where R n , is the n-th set of routes; V n , vector of vertices that belong in the bus route; Q n , vector of route passenger load for every vertex of the route; f n , route fleet; h n , headway between two vehicles of the route fleet; R n , the n-th family of sets of routes; S, the set of all possible routes between the graph vertices.
Given the above we consider two objective functions that represent the interests of the passengers as well as the interests of the operator. The passenger cost is represented through the generalized user cost equation, while the operator cost is represented through the fleet that is necessary for the network to operate. The objective equations and the overall goal are as follows: where T, is the total passenger cost; ϑ w , is the weight associated with the time users spend waiting at a bus stop; T w , the total time users spend waiting at a bus stop; T IVT , is the total time users spend in a vehicle; P 1t , P 1t are the total time penalties associated with users needing to transfer routes once or twice respectively; P ud , is the time penalty for unsatisfied demand by users; O, is the total operator cost, measured in vehicles; f R i n , is the total vehicle fleet required for the selected route set R i . Therefore, our aim is to find a set of Pareto-optimal network solutions (i.e., route sets) between these two conflicting objectives. Constraints t min , t max , minimum and maximum time-length of a route; w R n i , the time-length of the i-th route in the n-th route set; Q R n i , passenger load observed in the i-th route of the n-th route set; f R n i , fleet size the i-th route of the n-th route set; CAP R n i , vehicle passenger capacity for the i-th route of the n-th route set; LF, maximum load factor set by the operator; h min , h max , minimum and maximum headway of a route; h R i n , the headway of the i-th route in the n-th route set; r min , r max , maximum number of routes allowed in a route set:

Solution Approach-Methodology
The Transportation Network Design and Frequency Setting problem is well known to be non-convex and computationally NP-hard. In the scope of this study, a heuristic and meta-heuristic approach leading to a novel, alternating-objective, genetic algorithm has been chosen. Its goal is to produce a wide array of feasible network solutions that satisfy the demands of both the transportation network user as well as its operator while attempting to find near-optimal solutions. These solutions are known as Pareto Optimal solutions, as they describe the edge situation where each party's interests cannot be further improved without compromising the interests of the other party.
The proposed solution can be broken down into two distinct parts: (1) Route database & population generation (2) Alternating-objective genetic algorithm with elitism The first part utilizes heuristics-mainly Yen's k-shortest path [52] and Dijkstra's algorithm [53]-to create a vast database of routes connecting any two nodes of the OD matrix. Consequently, the population algorithm attempts to create feasible (i.e., fulfilling the constraints) network solutions by selecting routes generated in the preceding step. Each new feasible solution is called a "chromosome" while the set of solutions generated is called a "generation". In the second part, a modified genetic algorithm cycles between finding the best (i.e., lowest cost) solutions in the generation for the users and the best for the operator. The best solutions are kept, while new ones are created for the next generation. By the end of the algorithm, a Pareto frontier has been created, outlining a set of near-optimal solutions for both the users and the operator. A flowchart of the algorithm is presented in Figure 1.

Input Data
The GA requires a set of data matching the concepts developed under the problem formulation section. Therefore, the following data are required for the algorithm to run: (1) A symmetrical Origin-Destination matrix ( It should be noted that, while, the list of parameters and variables is extensive, the ones that define the robustness of the results are the OD and time matrices. The rest are merely tools that allow for the fine-tuning of the results and assist in drastically reducing the time required for such a resource-intensive algorithm to run. It is possible for all the parameters (i.e., 3a to 4g) to be defined as variables and produce a wider variety of results in exponentially more time.

Route Database & Population Generation
This first stage of the algorithm aims to create a wide range of routes connecting OD nodes, from which feasible solutions can be created. Naturally, the greater the size and variety within the pool, the greater the algorithm's ability to approach optimality. At the same time, a large pool size can significantly lower computational performance.
The initial routes are generated through Yen's k-shortest path algorithm and are subject to the constraints that have been described. The benefit of the k-shortest approach is that it allows for "k-1" shortest routes between two nodes to be created. An increase in the value of "k" is reflected by an increase in the influence of lines generated as, by varying from the shortest path, they cover more nodes.
Consequently, the "population" algorithm looks to bundle together routes in such a way that the transportation network constraints are met. This process involves randomly selecting routes from the route database so that a pre-set percentage of nodes are covered. Next, given the network created, the passenger loads outlined in the OD matrix are assigned to the routes (including 1-and 2-transfer routes should the OD pair not be covered directly) that service their destinations. These routes start off with pre-set headways and frequencies. Should more than one route service the same OD set, the following equation is used to calculate the possibility of passengers selecting each one: where P route , is the possibility of passengers selecting the ith solution of the route set. Following passenger assignment, new headways are calculated that, more accurately, service passenger demand. Given these new headways, user costs are recalculated (mainly waiting cost) and so are the route choice probabilities. This loop can continue perpetually, as each route choice is further optimized to match demand. In Figure 2, the average difference that each subsequent loop makes for each chromosome are shown, as tested on 1000 different chromosomes. Future Transp. 2021, 1, FOR PEER REVIEW 8 perpetually, as each route choice is further optimized to match demand. In Figure 2, the average difference that each subsequent loop makes for each chromosome are shown, as tested on 1000 different chromosomes. Lastly, in case any of the constrains and checks are violated, the route set is discarded and a new one is created. This process continues until the amount of feasible route sets created reaches the specified "chromosome number per generation" that the user has provided as input.

Alternating Objective Genetic Algorithm with Elitism
The alternating-objective genetic algorithm with elitism (AOGAE) is built to find optimal Pareto-efficient solutions within the pool generated during the previous step. This is achieved through iterative comparisons spanning different sets ("generations") of solutions and scripts that allow the best solutions to move onto later generations (elitism), thus gradually augmenting the median performance of each solution set. The process of the best solutions "surviving" while the rest are discarded has been inspired by Darwin's evolutionary theory and modern genetic processes; it is the reason that this algorithmic approach is called a "genetic algorithm".
In the scope of this paper, the algorithm was modified so that it alternates between minimizing user cost and minimizing operator cost. This allowed the methodology to retain only the solutions that minimized both costs at once and, finally, produce efficient networks. In the following paragraphs, the way the algorithm works will be outlined and points that allowed us to produce efficient results will be outlined.
The first step of the AOGAE algorithm is chromosome crossover. It aims to improve chromosome solutions by changing a route from one chromosome with a route from another from the same generation. Initially, each generation's chromosomes are sorted based on the current optimization objective. Each chromosome is assigned a probability of being the participating in a crossover instance. This serves to incentivize more radical changes in chromosomes with subpar performance. Lastly, in case any of the constrains and checks are violated, the route set is discarded and a new one is created. This process continues until the amount of feasible route sets created reaches the specified "chromosome number per generation" that the user has provided as input.

Alternating Objective Genetic Algorithm with Elitism
The alternating-objective genetic algorithm with elitism (AOGAE) is built to find optimal Pareto-efficient solutions within the pool generated during the previous step. This is achieved through iterative comparisons spanning different sets ("generations") of solutions and scripts that allow the best solutions to move onto later generations (elitism), thus gradually augmenting the median performance of each solution set. The process of the best solutions "surviving" while the rest are discarded has been inspired by Darwin's evolutionary theory and modern genetic processes; it is the reason that this algorithmic approach is called a "genetic algorithm".
In the scope of this paper, the algorithm was modified so that it alternates between minimizing user cost and minimizing operator cost. This allowed the methodology to retain only the solutions that minimized both costs at once and, finally, produce efficient networks. In the following paragraphs, the way the algorithm works will be outlined and points that allowed us to produce efficient results will be outlined.
The first step of the AOGAE algorithm is chromosome crossover. It aims to improve chromosome solutions by changing a route from one chromosome with a route from another from the same generation. Initially, each generation's chromosomes are sorted based on the current optimization objective. Each chromosome is assigned a probability of being the participating in a crossover instance. This serves to incentivize more radical changes in chromosomes with subpar performance.
Should a pair of chromosomes be selected, the segments connecting a random OD pair that is present in both chromosomes are interchanged in what is called a two-point crossover. To guarantee that the new chromosome remains within the boundaries set by the initial constraints, an edge recombination operator is used.
Another method used by the AOGAE to avoid local minima, is called "mutation". This process imitates the way that parts of genetic code randomly mutate to form the basis of evolutionary theory. This process is controlled by a global mutation probability for all chromosomes in each generation and a mutation probability related to the fitness of each chromosome. The worse the chromosome the more likely it is to undergo large scale mutations. Chromosomes selected for the mutation process exchange one of their routes with another from the route pool. Consequently, all the initial constraints are checked and, in the case that one is violated, the exchange is discarded and the process restarts.
Lastly, the AOGAE method implements a final way forcing the algorithm out of local minima by constantly monitoring the similarity of the results between consecutive generations. Should the similarity limit be exceeded, the algorithm discards the entirety of the generation and looks to create new chromosomes from the route pool. It should be noted that the maximum amount of consequent generations exhibiting increased similarity, heavily impacts algorithm results. This occurs, as discarding generations too early prevents them from reaching optimality, while discarding them too late may lead to results focused on local minima.
We suggest, relating max similarity occurrences with the number of generations that each instance of the algorithm will go through. All the processes described in this section alternate between minimizing user cost and minimizing operator cost. At the end of the algorithm, a wide array of solutions is presented that service different user/operator needs at performance levels that appear more robust than current literature.

Results
The AOGAE algorithm was tested on Mandl's benchmark network (Figure 3), consisting of 15 vertices and 21 bidirectional links with no capacity constraint. The network was based on the locations and roads between 15 swiss towns. Mandl's network was selected as it is the most used network for testing and benchmarking new approaches to the existing state of the art. By running the AOGAE algorithm on this network the authors benchmarked their results as presented in Table 1 The algorithm was tested on an Intel i5 7300HQ processor, 10 times with each run taking 10 min to complete. This reflects the potential of the algorithm to be integrated into widely adopted personal computer processors and therefore speaks to its potential ease of adoption from colleagues around the world. The results from the best run are presented in the following tables and figures. Future Transp. 2021, 1, FOR PEER REVIEW 10    The variables in Table 2 were used as inputs for the algorithm. The results of the AOGAE algorithm are compared with similar methodologies [12,22,27,31,34,35,49,[54][55][56] in Table 1. It should be noted that, while a host of research teams have benchmarked algorithms in the Mandl network, results are, more often than not, presented in different formats. An option to circumvent this divide would be to try and reproduce research results based on the route sets provided by the research teams. However, even in that case, considerable assumptions would have to be made about other characteristics of their networks (e.g., headways, maximum route detour factors etc.). Therefore, we chose to only present the results that were reported by our colleagues.
The comparison is performed using the following measures: (1) Fleet size-commonly the most significant measure of operator cost The results presented in Table 3, show improved solutions in regards to user cost, direct percentage and fleet size. The differences in average and maximum headway in comparison with preceding literature show that more optimal solutions can be found when the majority of network capacity is utilized. Additionally, the data visualization software package Tableau Public (Tableau Software, Seattle, WA, USA, 2017) was utilized as a method of enabling a more concise, appealing and interactive presentation of results. All such Figures can be found at the link provided at the Supplementary Materials Section. Figure 4 presents the observed improvement in direct network service that can be achieved through the AOGAE method while reducing fleet size. Naturally, very few chromosomes have second transfer percentages, and no chromosomes leave demand unsatisfied as the network is small in size; for larger networks it is considerably more likely to observe unsatisfied demand and increased percentages of 2-transfer trips. Other factors contributing to the distribution of trip service are the initial constraints regarding the minimum node coverage of each chromosome, the limits set on the number and time of routes and the penalties relating to transfers that each trip requires. achieved through the AOGAE method while reducing fleet size. Naturally, very few chromosomes have second transfer percentages, and no chromosomes leave demand unsatisfied as the network is small in size; for larger networks it is considerably more likely to observe unsatisfied demand and increased percentages of 2-transfer trips. Other factors contributing to the distribution of trip service are the initial constraints regarding the minimum node coverage of each chromosome, the limits set on the number and time of routes and the penalties relating to transfers that each trip requires.     Figure 6 presents the set of Pareto-optimal solutions that were generated via the AOGAE method. As the non-convex nature of the problem prohibits single global optimal solutions, here Pareto-efficient solutions for different fleet sizes are demonstrated. Larger (>84 vehicles) and smaller (<67 vehicles) fleet sizes have few or no solutions as they do not fit the greater Pareto optimality objective of the algorithm. A trend emerging from the results which should be highlighted is the diminishing, user cost, returns that an increase in fleet size provides. This can be attributed to the nearly complete (approx. 95%) direct coverage of the network. It also explains why results converge in a small part of potential fleet sizes, as shown in Figure 7 Figure 6 presents the set of Pareto-optimal solutions that were generated via the AOGAE method. As the non-convex nature of the problem prohibits single global optimal solutions, here Pareto-efficient solutions for different fleet sizes are demonstrated. Larger (>84 vehicles) and smaller (<67 vehicles) fleet sizes have few or no solutions as they do not fit the greater Pareto optimality objective of the algorithm. A trend emerging from the results which should be highlighted is the diminishing, user cost, returns that an increase in fleet size provides. This can be attributed to the nearly complete (approx. 95%) direct coverage of the network. It also explains why results converge in a small part of potential fleet sizes, as shown in Figure 7.  Figure 6 presents the set of Pareto-optimal solutions that were generated via the AOGAE method. As the non-convex nature of the problem prohibits single global optimal solutions, here Pareto-efficient solutions for different fleet sizes are demonstrated. Larger (>84 vehicles) and smaller (<67 vehicles) fleet sizes have few or no solutions as they do not fit the greater Pareto optimality objective of the algorithm. A trend emerging from the results which should be highlighted is the diminishing, user cost, returns that an increase in fleet size provides. This can be attributed to the nearly complete (approx. 95%) direct coverage of the network. It also explains why results converge in a small part of potential fleet sizes, as shown in Figure 7 Figure 7 also highlights the fact that the algorithm converges to a few different fleet sizes. Figure 8 presents the average passenger load per vehicle per generation, with color marking the percentage of trips covered by direct routes. It is shown that the algorithm favors higher-density, direct routes as a method of optimizing the network. This result is   Figure 8 presents the average passenger load per vehicle per generation, with color marking the percentage of trips covered by direct routes. It is shown that the algorithm favors higher-density, direct routes as a method of optimizing the network. This result is based on the assumption that passenger satisfaction does not vary with changes in the vehicle's level of service. Figure 7. Histogram of number of chromosomes per fleet size (color notes user cost, with blue and red marking low and high user costs respectively). Figure 7 also highlights the fact that the algorithm converges to a few different fleet sizes. Figure 8 presents the average passenger load per vehicle per generation, with color marking the percentage of trips covered by direct routes. It is shown that the algorithm favors higher-density, direct routes as a method of optimizing the network. This result is based on the assumption that passenger satisfaction does not vary with changes in the vehicle's level of service.    The most common fleet size being 71 shows that a fleet size of 71 is more likely to lead to Pareto-optimal solutions. This is also demonstrated in the user costs relative to fleet size that are represented by the color of each bar, with higher cost being marked by red. The reason that solutions with high fleet sizes may also have high user costs is that  The most common fleet size being 71 shows that a fleet size of 71 is more likely to lead to Pareto-optimal solutions. This is also demonstrated in the user costs relative to fleet size that are represented by the color of each bar, with higher cost being marked by red. The reason that solutions with high fleet sizes may also have high user costs is that they are mostly generated in the first generations of the algorithm and then quickly discarded as they do not reflect Pareto optimality on the operator side. Figure 10 is illustrating how more Pareto-optimal solutions come hand-in-hand with vehicles being more crowded; this is reflected by the Level of Service (LoS) index in Figure 10, which reflects the amount of crowding in each bus with a maximum of 1.25 times capacity as described in Table 2. The fluctuation in the line graphs is attributed to the algorithm cycling optimization objectives with each new generation.

Policy Implications
Despite its very promising results, implementing the algorithm to real-world scenarios remains a challenge as many of the assumptions that have been made would need to be re-evaluated. The main ones being the following: 1) Same operational speed in all network links 2) All links are bidirectional 3) Symmetrical OD-matrix 4) Bus is the single mode of public transport 5) No walking time to/from network nodes Additionally, any real-world application would need to be implemented in collaboration with a plethora of authorities responsible for various parts of the transportation network, including the city's administration, its transportation authority, and national/state representatives. They would provide, inform, and curate the algorithm's input data in order to improve the quality and applicability of results generated.
Nevertheless, transportation networks that have seen few changes over the decades, stand to gain greatly from applying state-of-the-art methods to their planning, both in terms of quality-of-service and in terms of minimizing stakeholder costs. As our comparisons show, networks designed through analytic approaches (e.g., Mandl's) are significantly less efficient in comparisons to more modern methods. Similar studies [7,29] that

Policy Implications
Despite its very promising results, implementing the algorithm to real-world scenarios remains a challenge as many of the assumptions that have been made would need to be re-evaluated. The main ones being the following: (1) Same operational speed in all network links (2) All links are bidirectional (3) Symmetrical OD-matrix (4) Bus is the single mode of public transport (5) No walking time to/from network nodes Additionally, any real-world application would need to be implemented in collaboration with a plethora of authorities responsible for various parts of the transportation network, including the city's administration, its transportation authority, and national/state representatives. They would provide, inform, and curate the algorithm's input data in order to improve the quality and applicability of results generated.
Nevertheless, transportation networks that have seen few changes over the decades, stand to gain greatly from applying state-of-the-art methods to their planning, both in terms of quality-of-service and in terms of minimizing stakeholder costs. As our comparisons show, networks designed through analytic approaches (e.g., Mandl's) are significantly less efficient in comparisons to more modern methods. Similar studies [7,29] that applied their methods to real-world scenarios have shown improvements ranging up to 30% in network performance metrics.

Discussion and Conclusions
In this paper, we propose an alternating-objective genetic algorithm that employs heuristics to calculate Pareto optimal solutions for the transportation network design and frequency setting problem. The first part of the algorithm deals with generating a wide pool of feasible solutions while the second part deals with finding the Pareto frontier within this pool. We showed that the algorithm provides high quality results within reasonable computational times and, consequently, compared them to recent literature. Lastly, the route to applying the algorithm to real-world applications and its benefits was analyzed and the parts in need of re-evaluation were pinpointed.
Despite its promising results, the research presented in this paper has a number of limitations that extend to its immediate applicability to real cities. The transportation network that is considered is simplified in its structure (i.e., same operational speed on all network links, all links are bidirectional, symmetrical OD-matrix) and properties (i.e., only bus as a mode of public transportation, no walking time to and from nodes) while multiple intrinsic characteristics of both passengers and operators and not considered in this study (e.g., walking time, operational costs other than fleet). Nevertheless, the results presented in this paper provide an easily applicable and computationally inexpensive tool to integrate in future research.
Future research should focus on including multiple modes of public transport and their interactions, deeper modelling of passenger behaviour and integrating transport econometrics into the operator cost functions. Additionally, while some existing approaches can be run on office computers, thus reducing the necessity of costly servers, additional research should focus on aligning existing computational transportation research with the state-of the-art of Information Technology practices, such as edge computing and containers.
The field of public transportation planning is lacking in both multimodal approaches and real-life applications. As inputs increase, a host of unpredicted variables could significantly skew results. Therefore, in the future, we aim to modify the methodology that has been described to effectively plan for n-mode choice networks (e.g., bus-metro-tramfunicular) and apply it to larger, more complex networks. Additionally, we intend to include vehicle level-of-service in the generalized user cost objective function, to better reflect changes in user perception as buses get increasingly crowded.
Lastly, a next step of this study would be to perform a case study for the city of Thessaloniki, Greece to test the algorithm and its hypotheses. Thessaloniki is a prime candidate for the AOGAE algorithm as it is in the midst of building new metro lines throughout the city while the restoration of tram transport is an ever-persistent discussion in city planning groups. project administration, T.V., K.G. and I.P.; funding acquisition, T.V., K.G., I.P. and G.G. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.