Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed optimization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

Co-simulation Framework for Community-scale Building-grid Integration [SWR-21-75]

Distributed energy resources (DERs), including rooftop solar, energy storage, and flexible loads, are gaining popularity as costs decline and as building owners and utilities realize their benefits. DERs can improve distribution system efficiency, help prevent the need for expensive grid upgrades, and increase the resilience of local communities. However, they can also cause difficulties in grid operations and can require controls to achieve their benefits. To address this challenge, NREL researchers have developed a community-scale solution that assesses the impacts of DERs and their control strategies on a distribution system. The framework has been shown to reduce solar photovoltaic (PV) curtailment to 0%, mitigate the adverse impact of solar variability on the distribution voltage, and provide up to 5-day critical load support during emergency events. Utilizing 5 different modules representing the feeder, buildings, home energy management systems, an aggregator, and a utility controller, NREL expects this simulation technology to play a critical role in the continued integration of DERs. According to the Energy Information Administration (EIA), solar curtailments accounted for 94% of the total energy curtailed in the California Independent System Operator (CAISO) in 2020. By enabling Independent System Operators (ISOs) and utility operators to bring solar curtailments to 0%, the electrical grid can become less dependent on fossil-fueled power generation sources. NREL's co-simulation framework contains five major components: Distribution Feeder Model: describes the distribution feeder topology using OpenDSS, including the locations of all DERs. Residential Building Model: simulates a large number of buildings at a high resolution using OCHRETM. The model is equipped to control equipment based on signals from an external module. The model includes major household appliances such as HVAC and a water heater, non-dispatchable load models, a distributed PV system, and a home battery system. Home Energy Management System: optimizes the controls for the devices in a home using foreseeTM. The control can adjust based on the user preferences including cost, comfort, and convenience. In hierarchical control scenarios, where the houses follow signals from an aggregator, the home energy management system provides a flexibility band with a range of power and follows the dispatch signals received from aggregator. Community-Level Aggregator: solves for optimal energy dispatch based on the flexibility bands received from each home and the grid service signal received from the utility controller. Utility-Level Controller: provides grid signals for voltage control using Distributed Energy Resources (DERs), such as solar systems, in the community.

Balamurugan, Sivasathya Pradha↗

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Seasonal Reconfiguration of Electrical Distribution Systems to Mitigate the Impact of Electric Vehicle Charging

Power grids face challenges in their infrastructure related to the integration of electric vehicles (EV). In particular, EV charging stations may induce instability in key system parameters such as substantial voltage drops, active power losses, and transformer overload due to high demand during charging periods. This article presents a seasonal reconfiguration strategy based on the differential evolution (DE) algorithm, aimed at enhancing system performance under highly variable and stochastic load profiles, particularly those driven by EV charging. The DEA algorithm is hybridized with the find-union (FU) algorithm to efficiently ensure network radiality throughout the optimization process. The proposed methodology is validated on a hybrid distribution system composed of the IEEE 33-bus network, a modified IEEE 13-bus system, and a specific 13-bus microgrid. Results have demonstrated that seasonal reconfiguration significantly reduces active power losses and mitigates transformer loading during critical demand hours, thereby quantifiably increasing the system’s performance. As an integral component of the proposed approach, an analysis of CO2 emissions associated with energy losses is included, allowing a contextualized assessment of the environmental benefits of seasonal reconfiguration in various geographical areas.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimal Control Strategy With Efficiency and Reliability Improvement for Offshore DC Microgrids

Offshore microgrids, due to their remote location and lack of external energy support, face significant challenges in wide-range load operation and maintenance. Consequently, efficiency and reliability are critical concerns for converters in offshore dc microgrids. This article presents an optimal control strategy aimed at enhancing both efficiency and reliability. A normalized nonlinear relationship between power loss and thermal stress of a paralleled converter is first established. Based on this, a dual-objective optimization function with an active weight function as well as a system overall performance index is established. The active weight function dynamically adjusts the control priority based on converter efficiency and switching device thermal stress. Then, the optimal power-sharing strategy is derived by the Lagrange multiplier method with the proposed optimal function. Additionally, to accommodate a wide load range, an optimal selection strategy for operating converter combinations is proposed, requiring only low-bandwidth communication. Experiment verification is given to validate the effectiveness of the proposed control strategy. The experiment results demonstrate that the proposed control strategy can improve the overall performance of offshore microgrids by optimizing efficiency and reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Lion Cub: Minimizing Communication Overhead in Distributed Lion

Communication overhead is a key challenge in distributed deep learning, especially on slower Ethernet intercon nects, and given current hardware trends, communication is likely to become a major bottleneck. While gradient compression techniques have been explored for SGD and Adam, the Lion optimizer has the distinct advantage that its update vectors are the output of a sign operation, enabling straightforward quantization. However, simply compressing updates for communication and using techniques like majority voting fails to lead to end-to-end speedups due to inefficient communication algorithms and reduced convergence. We analyze three factors critical to distributed learning with Lion: optimizing communication methods, identifying effective quantization methods, and assessing the necessity of momentum synchronization. Our findings show that quantization techniques adapted to Lion and selective momentum synchronization can significantly reduce communication costs while maintaining convergence. We combine these into Lion Cub, which enables up to 5x speedups in end-to-end training compared to Lion. This highlights Lion’s potential as a communication-efficient solution for distributed training.

97 MATHEMATICS AND COMPUTING↗

Vertical Characterization of Flow and Contaminant Concentration in New Hanford Extraction Wells

Effective groundwater remediation should focus on high-concentration, low-permeability zones within the aquifer, as these areas often present significant challenges to contaminant removal. Characterizing the vertical distribution of contamination in an aquifer is therefore critical for optimizing and monitoring pump-and-treat (P&T) systems to meet target cleanup timelines. Flow patterns driven by hydraulic gradients and variations in permeability are described to provide insight into contaminant mass removal from long-screened wells. The objective of this multi-phase work is to provide actionable recommendations for improving remedy performance monitoring and P&T operations, which ongoing evaluations will address. This report summarizes the preliminary results of vertical flow and contaminant distribution monitoring in new extraction wells at the Hanford Site.

54 ENVIRONMENTAL SCIENCES↗

Blueprinting Electrified Transit System Implementation

To achieve a more affordable and reliable transportation system, we need to smartly upgrade our power systems and install a large number of charging stations, but conventional planning methods are not up to the task. By applying advanced simulation and optimization tools, we can design a smarter, more cost-effective electric transportation network. The initial focus was on public transit systems, demonstrating how this approach can deliver broader economic, reliability, and air quality benefits nationwide.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Combined Retrievals of Boreal Forest Fire Aerosol Properties with a Polarimeter and Lidar

Absorbing aerosols play an important, but uncertain, role in the global climate. Much of this uncertainty is due to a lack of adequate aerosol measurements. While great strides have been made in observational capability in the previous years and decades, it has become increasingly apparent that this development must continue. Scanning polarimeters have been designed to help resolve this issue by making accurate, multi-spectral, multi-angle polarized observations. This work involves the use of the Research Scanning Polarimeter (RSP). The RSP was designed as the airborne prototype for the Aerosol Polarimetery Sensor (APS), which was due to be launched as part of the (ultimately failed) NASA Glory mission. Field observations with the RSP, however, have established that simultaneous retrievals of aerosol absorption and vertical distribution over bright land surfaces are quite uncertain. We test a merger of RSP and High Spectral Resolution Lidar (HSRL) data with observations of boreal forest fire smoke, collected during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS). During ARCTAS, the RSP and HSRL instruments were mounted on the same aircraft, and validation data were provided by instruments on an aircraft flying a coordinated flight pattern. We found that the lidar data did indeed improve aerosol retrievals using an optimal estimation method, although not primarily because of the constraints imposed on the aerosol vertical distribution. The more useful piece of information from the HSRL was the total column aerosol optical depth, which was used to select the initial value (optimization starting point) of the aerosol number concentration. When ground based sun photometer network climatologies of number concentration were used as an initial value, we found that roughly half of the retrievals had unrealistic sizes and imaginary indices, even though the retrieved spectral optical depths agreed within uncertainties to independent observations. The convergence to an unrealistic local minimum by the optimal estimator is related to the relatively low sensitivity to particles smaller than 0.1 ( m) at large optical thicknesses. Thus, optimization algorithms used for operational aerosol retrievals of the fine mode size distribution, when the total optical depth is large, will require initial values generated from table look-ups that exclude unrealistic size/complex index mixtures. External constraints from lidar on initial values used in the optimal estimation methods will also be valuable in reducing the likelihood of obtaining spurious retrievals.

Knobelspiesse, K.↗

The application of simulation modeling to the cost and performance ranking of solar thermal power plants

A computer simulation code was employed to evaluate several generic types of solar power systems (up to 10 MWe). Details of the simulation methodology, and the solar plant concepts are given along with cost and performance results. The Solar Energy Simulation computer code (SESII) was used, which optimizes the size of the collector field and energy storage subsystem for given engine-generator and energy-transport characteristics. Nine plant types were examined which employed combinations of different technology options, such as: distributed or central receivers with one- or two-axis tracking or no tracking; point- or line-focusing concentrator; central or distributed power conversion; Rankin, Brayton, or Stirling thermodynamic cycles; and thermal or electrical storage. Optimal cost curves were plotted as a function of levelized busbar energy cost and annualized plant capacity. Point-focusing distributed receiver systems were found to be most efficient (17-26 percent).

Rosenberg, L. S.↗

Energy characteristics of multi-chiller load distribution algorithms in a large office building

Here, this study evaluates the energy efficiency of multi-chiller systems in large office buildings, focusing on their optimization across various climate zones as defined by ASHRAE. Using EnergyPlus for simulations, the research examines five different load distribution algorithms in multi-chiller systems that range from one to ten chillers, aiming to understand their effectiveness in 15 distinct climate zones. The primary objectives of the study include identifying the energy efficiency of multi-chiller systems in each climate zone, determining the appropriate number of chillers for each zone, and evaluating the performance of the load distribution algorithms. Based on the U.S. Department of Energy’s commercial building model, the results suggest that multi-chiller systems can significantly reduce cooling energy consumption in various climates. Among the algorithms evaluated, the Sequential Uniform Part Load Ratio (SUPLR) algorithm demonstrates notable efficiency, especially in the 4A climate zone (Baltimore), where it achieves substantial energy savings. Applying the SUPLR algorithm in a multi-chiller setup with four chillers in this zone leads to an estimated 24.5 % reduction in energy usage, equivalent to 183 MW annually. The research indicates that a range of 3 to 5 chillers is typically optimal for most climate zones. In-depth analysis in the 4A climate zone highlights the importance of minimizing operation hours at low Part Load Ratios (PLR) to ensure that chillers operate at a high Coefficient of Performance (COP). This strategy underscores the potential of well-designed multi-chiller systems to reduce cooling energy demand, particularly in climates with transitional seasons. This study provides an overview of the energy-saving potential of multi-chiller systems, applicable across a variety of climatic scenarios.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Design and analysis of a wake steering controller with wind direction variability

Wind farm control strategies are being developed to mitigate wake losses in wind farms, increasing energy production. Wake steering is a type of wind farm control in which a wind turbine's yaw position is misaligned from the wind direction, causing its wake to deflect away from downstream turbines. Current modeling tools used to optimize and estimate energy gains from wake steering are designed to represent wakes for fixed wind directions. However, wake steering controllers must operate in dynamic wind conditions and a turbine's yaw position cannot perfectly track changing wind directions. Research has been conducted on robust wake steering control optimized for variable wind directions. In this paper, the design and analysis of a wake steering controller with wind direction variability is presented for a two-turbine array using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented wake model. First, the authors propose a method for modeling the turbulent and low-frequency components of the wind direction, where the slowly varying wind direction serves as the relevant input to the wake model. Next, we explain a procedure for finding optimal yaw offsets for dynamic wind conditions considering both wind direction and yaw position uncertainty. We then performed simulations with the optimal yaw offsets applied using a realistic yaw offset controller in conjunction with a baseline yaw controller, showing good agreement with the predicted energy gain using the probabilistic model. Using the Gaussian wake model in FLORIS as an example, we compared the performance of yaw offset controllers optimized for static and dynamic wind conditions for different turbine spacings and turbulence intensity values, assuming uniformly distributed wind directions. For a spacing of five rotor diameters and a turbulence intensity of 10 %, robust yaw offsets optimized for variable wind directions yielded an energy gain equivalent to 3.24 % of wake losses recovered, compared to 1.42 % of wake losses recovered with yaw offsets optimized for static wind directions. In general, accounting for wind direction variability in the yaw offset optimization process was found to improve energy production more as the separation distance increased, whereas the relative improvement remained roughly the same for the range of turbulence intensity values considered.

17 WIND ENERGY↗

Simultaneous Retrieval of Aerosol and Cloud Properties During the MILAGRO Field Campaign

Estimation of Direct Climate Forcing (DCF) due to aerosols in cloudy areas has historically been a difficult task, mainly because of a lack of appropriate measurements. Recently, passive remote sensing instruments have been developed that have the potential to retrieve both cloud and aerosol properties using polarimetric, multiple view angle, and multi spectral observations, and therefore determine DCF from aerosols above clouds. One such instrument is the Research Scanning Polarimeter (RSP), an airborne prototype of a sensor on the NASA Glory satellite, which unfortunately failed to reach orbit during its launch in March of 2011. In the spring of 2006, the RSP was deployed on an aircraft based in Veracruz, Mexico, as part of the Megacity Initiative: Local and Global Research Observations (MILAGRO) field campaign. On 13 March, the RSP over flew an aerosol layer lofted above a low altitude marine stratocumulus cloud close to shore in the Gulf of Mexico. We investigate the feasibility of retrieving aerosol properties over clouds using these data. Our approach is to first determine cloud droplet size distribution using the angular location of the cloud bow and other features in the polarized reflectance. The selected cloud was then used in a multiple scattering radiative transfer model optimization to determine the aerosol optical properties and fine tune the cloud size distribution. In this scene, we were able to retrieve aerosol optical depth, the fine mode aerosol size distribution parameters and the cloud droplet size distribution parameters to a degree of accuracy required for climate modeling. This required assumptions about the aerosol vertical distribution and the optical properties of the coarse aerosol size mode. A sensitivity study was also performed to place this study in the context of future systematic scanning polarimeter observations, which found that the aerosol complex refractive index can also be observed accurately if the aerosol optical depth is larger than roughly 0.8 at a wavelength of (0.555 m).

Knobelspiesse, K.↗

Protein intake is more stable than carbohydrate or fat intake across various US demographic groups and international populations

The optimal macronutrient composition of the diet is controversial and many adults attempt to regulate the intake of specific macronutrients for various health-related reasons. The objective was to compare stability and ranges of intakes of different macronutrients across diverse adult populations in the USA and globally. US dietary intake data from NHANES 2009–2014 were used to determine macronutrient intake as a percentage of total energy intake. Variability in macronutrient intake was estimated by calculating the difference between 75th and 25th percentile (Q3–Q1) IQRs of macronutrient intake distributions. In addition, intake data from 13 other countries with per capita gross domestic product (GDP) over $10,000 US dollars (USD) were used to assess variability of intake internationally since there are large differences in types of foods consumed in different countries. Protein, carbohydrate, and fat intake (NHANES 2009–2014) was 15.7 ± 0.1, 48.1 ± 0.1, and 32.9 ± 0.1% kcal, respectively, in US adults. The IQR of protein intake distribution (3.73 ± 0.11% kcal) was 41% of carbohydrate intake distribution (9.18 ± 0.20% kcal) and 58% of fat intake distribution (6.40 ± 0.14% kcal). The IQRs of carbohydrate and fat intake distributions were significantly (P <0.01) influenced by age and race; however, the IQR of protein intake was not associated with demographic and lifestyle factors including sex, race, income, physical activity, and body weight. International mean protein intake was 16.3 ± 0.2% kcal, similar to US intake, and there was less variation in protein than carbohydrate or fat intake. Protein intake of the US population and multiple international populations, regardless of demographic and lifestyle factors, was consistently ~16% of total energy, suggesting biological control mechanism(s) tightly regulate protein intake and, consequently, influence intake of other macronutrients and food constituents. Substantial differences in intake of the other macronutrients observed in US and international populations had little influence on protein intake. This trial was registered at the ISRCTN registry as ISRCTN46157745

dietary preferences↗

Multi-variance replica exchange SGMCMC for inverse and forward problems via Bayesian PINN

Physics-informed neural network (PINN) has been successfully applied in solving a variety of nonlinear non-convex forward and inverse problems. However, the training is challenging because of the non-convex loss functions and the multiple optima in the Bayesian inverse problem. In this work, we propose a multi-variance replica exchange stochastic gradient Langevin dynamics method to tackle the challenge of the multiple local optima in the optimization and the challenge of the multiple modal posterior distribution in the inverse problem. Replica exchange methods are capable of escaping from the local traps and accelerating the convergence; two chains with different temperatures are designed where the low temperature chain aims for the local convergence, and the target of the high temperature chain is to travel globally and explore the whole loss function entropy landscape. However, it may not be efficient to solve mathematical inversion problems by using the vanilla replica method directly since the method doubles the computational cost in evaluating the forward solvers (likelihood functions) in the two chains. To address this issue, we propose to make different assumptions on the energy function estimation and this facilities one to use solvers of different fidelities in the likelihood function evaluation. More precisely, one can use a solver with low fidelity in the high temperature chain while using a solver with high fidelity in the low temperature chain. Our proposed method significantly lowers the computational cost in the high temperature chain, meanwhile preserving the accuracy and converging very fast. Here we give an unbiased estimate of the swapping rate and give an estimation of the discretization error of the scheme. To verify our idea, we design and solve four inverse problems which have multiple modes. The proposed method is also employed to train the Bayesian PINN to solve the forward and inverse problems; faster and more accurate convergence has been observed when compared to the stochastic gradient Langevin dynamics (SGLD) method and vanilla replica exchange methods.

97 MATHEMATICS AND COMPUTING↗

Enhanced charge carrier extraction and transport with interface modification for efficient tin-based perovskite solar cells

Interface modification improves charge carrier extraction in tin-based perovskite solar cells. Tin-based perovskites have become the most promising non-lead perovskites due to their ideal band gap and low toxicity. Although the open circuit voltage of tin-based perovskite solar cells (TPSCs) continues to approach the theoretical value, the short-circuit current is still far from the theoretical value. Here, we describe an interface modification method by regulating the property of hole transport layer, PEDOT:PSS, which improves the surface molecular morphology and the energy level alignment of PEDOT:PSS/perovskite interface. Advanced GIWAXS and IR s-SNOM characterization are conducted to achieve multi-dimensional characterization of nanoscale surface morphology and chemical distribution of PEDOT:PSS. With the multi-attribute optimization, charge carrier extraction and non-radiative recombination are also improved. The resultant TPSCs exhibit a higher power conversion efficiency of 13.32% in compared with the control device of 10.50%, accompanied with an increase in the short-circuit current from 18.10 to 20.50 mA cm −2 and FF from 68.23% to 76.43%. This work demonstrates a reliable strategy for improving charge carrier extraction and device performance for lead-free TPSCs.

Zhao, Zhenzhu↗

Fast jet tagging with MLP-Mixers on FPGAs

We explore the innovative use of MLP-Mixer models for real-time jet tagging and establish their feasibility on resource-constrained hardware like FPGAs. MLP-Mixers excel in processing sequences of jet constituents, achieving state-of-the-art performance on datasets mimicking Large Hadron Collider conditions. By using advanced optimization techniques such as High-Granularity Quantization and Distributed Arithmetic, we achieve unprecedented efficiency. These models match or surpass the accuracy of previous architectures, reduce hardware resource usage by up to 97%, double the throughput, and half the latency. Additionally, non-permutation-invariant architectures enable smart feature prioritization and efficient FPGA deployment, setting a new benchmark for machine learning in real-time data processing at particle colliders.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Medium-Duty Delivery Truck Integrated Bidirectional Wireless Power Transfer System With Grid and Stationary Energy Storage System Connectivity

Electric vehicles (EVs) can provide power to the grid or buildings similar to distributed energy resources (DER) for energy balancing applications or optimizing the operation of the microgrids in harmony with the other DER assets. This article presents the operating modes of a bidirectional wireless power transfer (WPT) system designed for a medium-duty package delivery vehicle. The WPT system designed for this study can transfer 20 kW of power across 11 in of air gap using custom-designed double-D (DD) couplers with LCC–LCC tuning networks. The proposed system utilizes a 480-V three-phase grid connection, a plug-in hybrid delivery truck with bidirectional WPT, and a stationary energy storage system (SESS) that can be connected to the primary-side dc link. Due to the differences in primary and secondary dc bus voltages, and considering the voltage of the SESS, asymmetric voltage gains were used in the system. Sensitivity analyses of this system with respect to these voltage levels are presented. Five different operating modes of the grid, SESS, and the EV battery are investigated with experimental results. Control algorithms are described for grid-to-vehicle (G2V) and vehicle-to-grid (V2G) operating modes. A bidirectional WPT system is operated with a power factor of 0.99 on the grid side in every operating mode. The EV battery was charged with 20.3 kW with an overall efficiency of 93.02% in the G2V operating mode. Finally, in the V2G operating mode, the WPT system provided 12.82 kW of power back to the grid with an overall efficiency of 89.08%.

25 ENERGY STORAGE↗

Bug fixes / enhancements to Parallel tempering Markov chain Monte Carlo sampler (PTMCMCSampler) v3

Markov chain Monte Carlo sampler useful for a variety of problems for optimization, numerical integration, and generating draws from a probability distribution. They are often used in physical and mathematical problems, especially when it is difficult or impossible to use other approaches. This disclosure is for new bug fixes and enhancements to an existing open source library that is already in the public domain under the MIT license. https://github.com/jellis18/PTMCMCSampler

Forrer, Mark↗