## OPTIMISTIC PARALLEL DISCRETE EVENT SIMULATION OF THE EVENT-BASED TRANSMISSION LINE MATRIX METHOD

### Cached

### Download Links

- [www.informs-sim.org]
- [cell-devs.sce.carleton.ca]
- [www.sce.carleton.ca]
- [www.thesimguy.com]
- DBLP

### Other Repositories/Bibliography

### BibTeX

@MISC{Henderson_optimisticparallel,

author = {S. G. Henderson and B. Biller and M. -h. Hsieh and J. Shortle and J. D. Tew and R. R. Barton},

title = {OPTIMISTIC PARALLEL DISCRETE EVENT SIMULATION OF THE EVENT-BASED TRANSMISSION LINE MATRIX METHOD},

year = {}

}

### OpenURL

### Abstract

In this paper we describe a technique for efficient parallelization of digital wave guide network (DWN) models based on an interpretation of the finite difference time domain (FDTD) method for discrete event simulation. Modeling methodologies based on FDTD approaches are typically constrained in both the spatial and time domains. This interpretation for discrete event simulation allows us to investigate the performance of DWN models in the context of optimistic parallel discrete event simulation employing reverse computation for rollback support. We present parallel performance results for a large-scale simulation of a 3D battlefield scenario, 100km 2 and at a height of 100m with a resolution of 100m in the X-, Y-planes, and 10m in the Z-plane for 754 simultaneous radio wave transmissions. 1

### Citations

888 | Virtual time - Jefferson - 1985 |

168 |
Calendar Queues: a fast O(1) priority queue implementation for the simulation event set problem
- Brown
- 1988
(Show Context)
Citation Context ...Also, each event on the wavefront surface has the same timestamp, therefore the choice of priority queue used in the simulator has a major impact on performance, as shown in Figure 3. Calendar queue (=-=Brown 1988-=-) performance is known to have a complexity of O(n) when the frequency of ties among events is high, accounting for it’s poor performance in all cases. The performance of the Heap (Deo and Prasad 1992... |

121 | ber die partiellen differenzengleichungen der mathematischen physik - Courant, Friedrichs, et al. - 1928 |

64 |
Comparative analysis of periodic state saving techniques in Time Warp simulators
- Fleischmann, Wilsey
- 1995
(Show Context)
Citation Context ...nisms include copy state saving, incremental state-saving (Steinman 1993, Gomes 1996), and a large amount of effort has been expended on infrequent state saving (Lin and Preiss 1991, Lin et al. 1993, =-=Fleischmann and Wilsey 1995-=-). More recently, the technique of reverse computation has been suggested as an approach to rollback that does not rely on storing state information (Carothers, Perumalla, and Fujimoto 1999). Here, LP... |

60 | Efficient optimistic parallel simulations using reverse computation - Carothers, Perumalla, et al. - 1999 |

59 | Selecting the Checkpoint Interval in Time Warp Simulation
- Lin, Preiss, et al.
- 1993
(Show Context)
Citation Context ...to rollback mechanisms include copy state saving, incremental state-saving (Steinman 1993, Gomes 1996), and a large amount of effort has been expended on infrequent state saving (Lin and Preiss 1991, =-=Lin et al. 1993-=-, Fleischmann and Wilsey 1995). More recently, the technique of reverse computation has been suggested as an approach to rollback that does not rely on storing state information (Carothers, Perumalla,... |

57 | On the Partial Difference Equations of Mathematical Physics - 47Courant, Friedrichs, et al. - 1967 |

37 | Optimal memory management for Time Warp parallel simulation
- Lin, Preiss
- 1991
(Show Context)
Citation Context ...The major approaches to rollback mechanisms include copy state saving, incremental state-saving (Steinman 1993, Gomes 1996), and a large amount of effort has been expended on infrequent state saving (=-=Lin and Preiss 1991-=-, Lin et al. 1993, Fleischmann and Wilsey 1995). More recently, the technique of reverse computation has been suggested as an approach to rollback that does not rely on storing state information (Caro... |

36 | A comparative study of parallel and sequential priority queue algorithms
- Rönngren, Ayani
- 1997
(Show Context)
Citation Context ... O(n) when the frequency of ties among events is high, accounting for it’s poor performance in all cases. The performance of the Heap (Deo and Prasad 1992) is a large improvement, but the Splay Tree (=-=Rönngren and Ayani 1997-=-) implementation clearly attained the best performance at almost 1.3 million events per second on a single processor core. The main reason for comparing the model performance using different priority ... |

33 | ROSS: A High-Performance, Low Memory, Modular Time Warp System - Carothers, Bauer, et al. - 2002 |

31 | Computing global virtual time in shared-memory multiprocessors
- Fujimoto, Hybinette
- 1997
(Show Context)
Citation Context ...algorithm for shared memory multiprocessors when executed on a single compute node, and the Seven O’clock GVT algorithm when executed across multiple compute nodes (Carothers, Bauer, and Pearce 2002, =-=Fujimoto and Hybinette 1997-=-, Bauer et al. 2005). The ROSS simulator is a general-purpose simulator exposing a simple API for the creation of PDES models. Kernel Processes were introduced in Carothers, Bauer, and Pearce (2002) a... |

30 |
Incremental state saving in SPEEDES using C
- Steinman
- 1993
(Show Context)
Citation Context ...tatesaving has been employed as the rollback mechanism, outlined by the numerous studies in this area. The major approaches to rollback mechanisms include copy state saving, incremental state-saving (=-=Steinman 1993-=-, Gomes 1996), and a large amount of effort has been expended on infrequent state saving (Lin and Preiss 1991, Lin et al. 1993, Fleischmann and Wilsey 1995). More recently, the technique of reverse co... |

19 |
Parallel heap: An optimal parallel priority queue
- Deo, Prasad
- 1992
(Show Context)
Citation Context ...r queue (Brown 1988) performance is known to have a complexity of O(n) when the frequency of ties among events is high, accounting for it’s poor performance in all cases. The performance of the Heap (=-=Deo and Prasad 1992-=-) is a large improvement, but the Splay Tree (Rönngren and Ayani 1997) implementation clearly attained the best performance at almost 1.3 million events per second on a single processor core. The main... |

18 | A survey of the finite-difference timedomain literature
- Shlager, Schneider
- 1998
(Show Context)
Citation Context ... the discrete event paradigm: particle-in-cell (PIC) and Event-Based Transmission Line Matrix (ETLM). The PIC method is derived from the numerical method called finite difference time domain or FDTD (=-=Shlager and Schneider 1995-=-). This approach is commonly used when modeling physical systems because it is highly accurate and provides a bounding of the error in the computed result. This error is related to the representation ... |

17 | A new asynchronous methodology for modeling of physical systems: breaking the curse of courant condition
- Karimabadi, Driscoll, et al.
(Show Context)
Citation Context ...967). In the past few years, researchers have begun to adapt these methods to the discrete event paradigm. Notably, a 2005 study applied discrete event simulation to the particlein-cell (PIC) method (=-=Karimabadi et al. 2005-=-) and achieved a two order of magnitude increase in performance. The PIC 1-4244-1306-0/07/$25.00 ©2007 IEEE David W. Bauer Jr. Ernest H. Page The MITRE Corporation 7525 Colshire Drive McLean, VA 22102... |

16 | Specification of Discrete Event Models for Fire Spreading. SIMULATION: Transactions of the Society for Modeling and Simulation
- Muzy, Innocent, et al.
- 2005
(Show Context)
Citation Context ...rocessors. The second interpretation of a numerical method for physical systems for DES is called Event-Based Transmission Line Matrix or ETLM. This method has been used to study both fire-spreading (=-=Muzy et al. 2005-=-) as well as electromagnetic wave propagation (Kuruganti and Nutaro 2006). This method is based on a discrete approximation of a continuous structure. This relaxes the constraint on the distance that ... |

15 | L.: Numerical solution of 2-dimensional scattering problems using a transmission-line matrix - Johns, Beurle - 1971 |

11 | Particle-in-Cell Charged-Particle Simulations, Plus MonteCarlo Collisions With Neutral Atoms - Birdsall - 1991 |

10 | The Transmission-line matrix method—Theory and applications - Hoefer - 1985 |

9 |
Optimizing incremental state-saving and restoration
- Gomes
- 1996
(Show Context)
Citation Context ...been employed as the rollback mechanism, outlined by the numerous studies in this area. The major approaches to rollback mechanisms include copy state saving, incremental state-saving (Steinman 1993, =-=Gomes 1996-=-), and a large amount of effort has been expended on infrequent state saving (Lin and Preiss 1991, Lin et al. 1993, Fleischmann and Wilsey 1995). More recently, the technique of reverse computation ha... |

8 | Seven-o’clock: A new distributed gvt algorithm using network atomic operations
- Bauer, Yaun, et al.
- 2005
(Show Context)
Citation Context ...ultiprocessors when executed on a single compute node, and the Seven O’clock GVT algorithm when executed across multiple compute nodes (Carothers, Bauer, and Pearce 2002, Fujimoto and Hybinette 1997, =-=Bauer et al. 2005-=-). The ROSS simulator is a general-purpose simulator exposing a simple API for the creation of PDES models. Kernel Processes were introduced in Carothers, Bauer, and Pearce (2002) and aggregate the LP... |

6 |
The solution of Maxwell’s equations in three space dimensions and time by the TLM method of numerical analysis
- Akhtarzad, Johns
- 1975
(Show Context)
Citation Context ...oposed a numerical modeling technique based on Huygen’s model for solving 2-D scattering problems. This work was followed up with extensions to three dimensions and the concept of dielectric loading (=-=Akhtarzed and Johns 1975-=-, Johns 1974c, Johns 1974b, Akhtarzad and Johns 1974). Numerous researchers since Bauer and Page 677 have added features and extended the method to include concepts such as variable mesh size, simplif... |

6 | Event reconstruction in time warp, in - Li, Tropper - 2004 |

5 |
Solution of 6-components electromagnetic fields in three dimensions and time by the T.L.M
- Akhtarzad, Johns
- 1974
(Show Context)
Citation Context ...n’s model for solving 2-D scattering problems. This work was followed up with extensions to three dimensions and the concept of dielectric loading (Akhtarzed and Johns 1975, Johns 1974c, Johns 1974b, =-=Akhtarzad and Johns 1974-=-). Numerous researchers since Bauer and Page 677 have added features and extended the method to include concepts such as variable mesh size, simplified nodes, error correction techniques, and anisotro... |

4 | A case study in understanding OSPF and BGP interactions using efficient experiment design
- Bauer, Yuksel, et al.
- 2006
(Show Context)
Citation Context ...ed event processed. Several reverse computation models have been designed for a variety of models, including communications network modeling (Yaun et al. 2003, Yaun, Carothers, and Kalyanaraman 2003, =-=Bauer et al. 2006-=-), and physical systems (Tang et al. 2006). Additional optimizations for reducing memory consumption related to the storing of processed events was proposed in Bauer and Page (2007) for reverse comput... |

4 | Omelchenko, Optimistic simulations of physical systems using reverse computation, Simulation 82
- Tang, Perumalla, et al.
- 2006
(Show Context)
Citation Context ...ation models have been designed for a variety of models, including communications network modeling (Yaun et al. 2003, Yaun, Carothers, and Kalyanaraman 2003, Bauer et al. 2006), and physical systems (=-=Tang et al. 2006-=-). Additional optimizations for reducing memory consumption related to the storing of processed events was proposed in Bauer and Page (2007) for reverse computation and Li and Tropper (2004) for state... |

3 | The Solution of Inhomogeneous Waveguide Problems Using a Transmission-Line Matrix,” Microwave Theory and Techniques - Johns - 1974 |

3 | A discrete event method for wave simulation - Nutaro - 2006 |

2 | English translation of the 1928 German original - Journal |

2 | A new mathematical model to describe the physics of propagation, Radio Electron - Johns - 1974 |

1 | Application of the transmission-line matrix method to homogeneous waveguides of arbitrary crosssection - Johns |

1 |
Validation radio wave propagation model. Oak Ridge National Laboratory
- Kuruganti, Nutaro
- 2006
(Show Context)
Citation Context ...physical systems for DES is called Event-Based Transmission Line Matrix or ETLM. This method has been used to study both fire-spreading (Muzy et al. 2005) as well as electromagnetic wave propagation (=-=Kuruganti and Nutaro 2006-=-). This method is based on a discrete approximation of a continuous structure. This relaxes the constraint on the distance that can be used between points in the grid used to model the spatial domain.... |

1 |
Large-scale Bauer and Page network simulation techniques: examples of tcp and ospf models
- Yaun, Bauer, et al.
- 2003
(Show Context)
Citation Context ...omputing the inverse operations for each improperly executed event processed. Several reverse computation models have been designed for a variety of models, including communications network modeling (=-=Yaun et al. 2003-=-, Yaun, Carothers, and Kalyanaraman 2003, Bauer et al. 2006), and physical systems (Tang et al. 2006). Additional optimizations for reducing memory consumption related to the storing of processed even... |

1 | is a Senior Simulation Systems Engineer at the MITRE Corporation. He earned a Ph.D.,M.S., and B.S. from Rensselaer Polytechnic Institute in 2005, 2004, and 2000, respectively. Prior to joining MITRE, he was a research scientist at AT&T and GE. His researc - JR |