Splines and quasi-interpolation operators are important both in approximation theory and applications. In this paper, we construct a family of quasi-interpolation operators for the bivariate quintic spline spaces S₅³ (∆_mn⁽²⁾). Moreover, the properties of the proposed quasi-interpolation operators are studied, as well as its applications for solving two-dimensional Burgers’ equation and image reconstruction. Some numerical examples show that these methods, which are easy to implement, provide accurate results.

Because Canadian Nitrogen (N) fertilizer use has more than doubled since 1990 (1.2 to 2.9 MtN by 2021), a better understanding of this trend is needed. A comprehensive set of Recommended N Rates (RNR) that agreed with fertilizer sales data from 1996 and 2001 was compared with the Fertilizer Use Survey (FUS). The FUS was conducted from 2014 to 2021, with 2017 being the most representative year for these data. Using non-parametric statistics, confidence intervals were derived from the histograms used to present the FUS data. N application rates from the RNR for Canola, Spring and Duram wheat, and Oats in the west were all below their respective FUS confidence intervals, whereas N application rates for Grain corn showed almost no difference in N use between the RNR and FUS. Crop-specific N application rates interpolated from the RNR and FUS were integrated over their respective crop areas and plotted against national fertilizer sales records from 1990 to 2021. The rapid increase in N use between 2001 and 2017 (0.89 MtN), 90% of it (0.80 MtN) in Western Canada, was primarily due to increased application rates per crop, rather than crop area changes. The FUS-RNR interpolations were a good approximation of N sales records and could improve farm GHG emissions modelling. The economically important crops in Western Canada should be the main focus for N-related GHG reduction measures, but production losses need to be avoided.

This paper aims to demonstrate how the detection of communities for epidemiological networks can be used to reduce the dimensionality of optimal control problems. A range of planted partition stochastic blockmodels are generated, where the underlying communities in the model vary in detectability. A high-performance computing workflow is then used to simulate and estimate reduced model dynamics for the networks, which in turn is used to find optimal control solutions.

To reduce solar power production invariance, it is critical to study varying patterns of power production in the concerned region. Analyzing the patterns of past power production trends can help simulate power production scenarios for future. The current study area is around Amsterdam, located in Netherlands. PVoutput.org website is used to mine 6 months of solar power production data for 120 stations around Amsterdam city. FME Workbench software is used to actively fetch the data from the mentioned website and manage in a MySQL database. Solar attenuation maps created using ArcGIS, helped to graphically visualize the variations in solar power production at different times and locations. Further, spatial autocorrelation is checked between the stations using semi-variograms in geostatistical tool of ArcMap. This feature allows to check whether the stations located close to each other are more correlated to each other rather than stations which are far apart. The statistical data analysis of power production can aid solar power production companies to better interpolate and predict solar power in advance for the concerned study region.

Nucleation is a fundamental and general process at the initial stage of first order phase transition. Although various models based on the classical nucleation theory (CNT) have been proposed to explain the energetics and kinetics of nucleation, detailed understanding on nanoscales are still required. Here we focus on homogeneous bubble nucleation, in which evaluation of the size dependence of bubble free energy is the key issue. We propose application of a formula in stochastic thermodynamics, the Jarzynski’s equality, for data analysis of molecular dynamics (MD) simulation to evaluate the free energy of bubble nucleation. As a test case, we performed a series of MD simulations with Lennard-Jones (LJ) fluid system. By applying an external spherical force field to equilibrated LJ liquid, we evaluated the free energy change during bubble growth as the Jarzynski’s ensemble average of required works. A fairly smooth free energy curve was obtained as a function of bubble radius in metastable liquid of mildly negative pressure conditions.

We developed a novel approach in the field of spatiotemporal modelling, based on the spatialisation of time: the Timescape algorithm. It is especially aimed at sparsely distributed datasets in ecological research, whose spatial and temporal variability is strongly entangled. The algorithm is based on the definition of a spatiotemporal distance that incorporates a causality constraint and that is capable of accommodating the seasonal behaviour of the modelled variable as well. The actual modelling is conducted exploiting any established spatial interpolation technique, substituting the ordinary spatial distance with our Timescape distance, thus sorting, from the same input set of observations, those causally related to each estimated value at a given site and time. The notion of causality is expressed topologically and it has to be tuned for each particular case. The Timescape algorithm originates from the field of stable isotopes spatial modelling (isoscapes), but in principle it can be used to model any real scalar random field distribution.

Weak signal spectral detection is a serious problem in many applications, especially the target possesses certain frequency band and under low signal-to-noise ratio (SNR). A kind of novel technique based on pre-whitening scale-transformation symmetric bistable parallel stochastic resonance (SR) is proposed to solve the problem. Firstly, the pre-whitening can ensure the even distribution of the receiving signal to fit the requirement of SR processing. Secondly, the scale- transformation can help to utilize the weak signal improvement property especially under low frequency band more effectively. Thirdly, the symmetric bistable parallel SR can try to reduce the variance of the clutter and additive noise fluently. Fourthly, by subtracting the steady state re-sponse of the selected symmetric bistable SR system from the parallel SR output, the weak signal spectral detection objective can be realized successfully. Simulation results based on the real sea clutter radar data guarantee the effectiveness and applicability of the proposed symmetric bistable PSR processing approach.

This paper investigates the potential of utilising multilevel quasi-interpolation techniques on Chebyshev sparse grids for complex numerical computations. The paper starts by laying down the motivations for choosing Chebyshev sparse grids and quasi-interpolation methods with Gaussian kernels. It delves into the practical aspects of implementing these techniques. Various numerical experiments are performed to evaluate the efficiency and limitations of the multilevel quasi-sparse interpolation methods with dimensions two dimension and three dimension. The work ultimately aims to provide a comprehensive understanding of the computational efficiency and accuracy achievable through this approach, comparing its performance with traditional methods.

Subsurface temperature is the key parameter in geothermal exploration. An accurate estimation of the reservoir temperature is of high importance and usually done either by interpolation of borehole temperature measurement data or numerical modeling. However, temperature measurements at depths which are of interest for deep geothermal applications (usually deeper than 2 km) are generally sparse. A pure interpolation of such sparse data always involves large uncertainties and usually neglects knowledge of the 3D reservoir geometry or the rock and reservoir properties governing the heat transport. Classical numerical modeling approaches at regional scale usually only include conductive heat transport and do not reflect thermal anomalies along faults created by convective transport. These thermal anomalies however are usually the target of geothermal exploitation. Kriging with trend does allow including secondary data to improve the interpolation of the primary one. Using this approach temperature measurements of depths larger than 1,000 m of the federal state of Hessen/Germany have been interpolated in 3D. A 3D numerical conductive temperature model was used as secondary information. This way the interpolation result reflects thermal anomalies detected by direct temperature measurements as well as the geological structure. This results in a considerable quality increase of the subsurface temperature estimation.

Quasi-interpolation is a powerful tool for approximating functions using radial basis functions (RBFs) such as Gaussian kernels. This avoids solving large systems of equations as in RBF interpolation. However, quasi-interpolation with Gaussian kernels on compact intervals can have significant errors near the boundaries. This paper proposes a quasi-interpolation method with Gaussian kernels using Chebyshev points and boundary corrections to improve the approximation near the boundaries. The boundary corrections use a linear approximation of the function beyond the interval to estimate the truncation error and add correction terms. Numerical studies on test functions show that the proposed method reduces errors significantly near boundaries compared to quasi-interpolation without corrections, for both equally spaced and Chebyshev points. The convergence and accuracy with the boundary corrections are generally better with Chebyshev points compared to equally spaced points. The proposed method provides an efficient way to perform quasi-interpolation on compact intervals while controlling the boundary errors.

Searching over a sorted list is a classical problem in computer science. Binary Search takes at most log2n+1 tries to find an item in a sorted list of size n. Interpolation Search achieves an average time complexity of Ο(loglogn) for uniformly distributed data. Hybrids of Binary Search and Interpolation Search are also available to handle data with unknown distributions. This paper analyzes the computation cost of these methods and shows that interpolation can significantly affect their performance — accordingly, a new method, Interpolation Once Binary Search (IOBS), is proposed. The experimental results show that IOBS outperforms the hybrids of Binary Search and Interpolation Search for nonuniformly distributed data.

Gain-scheduled autopilots have emerged as a dominant strategy to achieve adaptive control of coupled, non-linear engineering complexities, owing to an ability to adapt to changing operational conditions and uncertainties. This study focuses on utilizing bilinear interpolation of gain-scheduled autopilots, emphasizing enhanced system performance and robustness. Through a comprehensive investigation and comparative analysis using three disparate cases, advantages over conventional methods are revealed. Strengths and weaknesses of both simple and specialized variants (such as linear, and real-time gain-scheduling) are introduced. Three missile guidance case–studies utilize simulation time and miss distance figures of merit. Potential to achieve precise control across various mission scenarios, while ensuring reduced computational complexity is revealed by nearly two–hundred percent improved missile miss distances with comparable distances traveled and slightly improved computational burden.

In this study, a cholera infection model with the bilinear infection rate is developed by considering the perturbation of the infection rate by the mean-reverting process. First of all, we give the existence of a globally unique positive solution for a stochastic system at arbitrary initial value. On this basis, the sufficient condition for the model to have an ergodic stationary distribution is given by constructing proper Lyapunov functions and tight sets. This indicates in a biological sense the long-term persistence of cholera infection. Furthermore, after transforming the stochastic model to an relevant linearised system, an accurate expression for the probability density function of the stochastic model around a quasi-endemic equilibrium is derived. Subsequently, the sufficient condition to make the disease extinct is also derived. Eventually, the theoretical findings are shown by numerical simulations. Numerical simulations show that the impact of regression speed and fluctuation intensity on stochastic systems.

The paper presents the problem of damage detection in steel girders. Static displacements at the selected point of the structure play the role of measured variables. Structural response signal decomposition is performed according to the Mallat pyramid algorithm, which is used to perform the discrete wavelet transform (DWT). This procedure allows us to quite well determine the location of structural damage. The geometry and the placement of any defective part of the structure may have a random character. It can be assumed that the random processes occurring in the broadly understood structure mechanics are Gaussian in nature. The first four probabilistic moments are estimated using three approaches independent: semi-analytical (SAM), perturbative (SPT), and Monte-Carlo simulations (MCS). The semi-analytical random approach seems to be the most optimal due to the necessary computation time. The incorporation of the mathematical stochastic apparatus into the classical (deterministic) analysis of the statics of the structure makes it possible to estimate the reliability measures of the analyzed girder.

High Efficiency Video Coding (HEVC) uses an 8-point filter and a 7-point filter, which are based on the discrete cosine transform (DCT), for the 1/2-pixel and 1/4-pixel interpolations, respectively. In this paper, discrete sine transform (DST)-based interpolation filters (IF) are proposed. The first proposed DST-based IFs (DST-IFs) use 8-point and 7-point filters for the 1/2-pixel and 1/4-pixel interpolations, respectively. The final proposed DST-IFs use 12-point and 11-point filters for the 1/2-pixel and 1/4-pixel interpolations, respectively. These DST-IF methods are proposed to improve the motion-compensated prediction in HEVC. The 8-point and 7-point DST-IF methods showed average BD-rate reductions of 0.7% and 0.3% in the random access (RA) and low delay B (LDB) configurations, respectively. The 12-point and 11-point DST-IF methods showed average BD-rate reductions of 1.4% and 1.2% in the RA and LDB configurations for the Luma component, respectively.

The article is devoted to the new results of the author, which add his previously published ones, of studying hidden rules and symmetries in structures of long single-stranded DNA sequences in eukaryotic and prokaryotic genomes. The author uses the existence of different alphabets of n-plets in DNA: the alphabet of 4 nucleotides, the alphabet of 16 douplets, the alphabet of 64 triplets, etc. Each of such DNA alphabets of n-plets can serve for constructing a text as a chain of these n-plets. Using this possibility, the author represents any long DNA nucleotide sequence as a bunch of many so-called n-texts, each of which is written on the basis of one of these alphabets of n-plets. Each of such n-texts has its individual percents of different n-plets in its genomic DNA. But it turns out that in such multi-alphabetical or multilayer presentation of each of many genomic DNA, analyzed by the author, universal rules of probabilities and symmetry exist in interrelations of its different n-texts regarding their percents of n-plets. In this study, the tensor product of matrices and vectors is used as an effective analytical tool borrowed from the arsenal of quantum mechanics. Some additions to the topic of algebra-holographic principles in genetics are also presented. Taking into account the described genomic rules of probability, the author puts also forward a concept of the important role of stochastic resonances in genetic informatics.

Here we are concerning with two problems of a coupled system of random and stochastic nonlinear differential equations with two coupled systems of nonlinear nonlocal random and stochastic integral conditions. The existence of solutions will be studied. The sufficient condition for the uniqueness of the solution will be given. The continuous dependence of the unique solution on the nonlocal conditions will be proved.

Some physicists surmise that gravity lies outside of quantum mechanics. Thus theories like the standard semiclassical theory of quantum to gravity coupling (that of Rosenfeld and Møller) are possible real models of interaction, rather than a mere approximation of a theory of quantum gravity. Unfortunately, semiclassical gravity creates inconsistencies such as superluminal communication. Alternatives by authors such as Diósi, Martin, Penrose, and Wang often use the term 'stochastic' to set themselves apart from the standard semiclassical theory. These theories couple to fluctuations caused by for instance continuous spontaneous localization, hence the term 'stochastic'. This paper looks at stochastic gravity in the framework of a class of emergent or ontological quantum theories, such as those by Bohm, Cetto, and de Broglie. It is found that much or all of the trouble in connecting gravity with a microscopic system falls away, as Einstein's general relativity is free to react directly with the microscopic beables. The resulting continuous gravitational wave radiation by atomic and nuclear systems does not, in contrast to Einstein's speculation, cause catastrophic problems. The small amount of energy exchanged by gravitational waves may have measurable experimental consequences. A very recent experiment by Vinante et al. performed on a small cantilever at mK temperatures shows a surprising non-thermal noise component, the magnitude of which is consistent with the stochastic gravity coupling explored here.

We suppose that a Cooper pair (CP) will experience a damping force exerted by the condensed matter. A Langevin equation of a CP in two dimensional condensed matter is established. Following a method similar to Nelson's stochastic mechanics, generalized Schr\"{o}dinger equation of a CP in condensed matter is derived. If the CPs move with a constant velocity, then the corresponding direct current (DC) electrical conductivity can be calculated. Therefore, a Drude like formula of resistivity of CPs is derived. We suppose that the damping coefficient of CPs in two dimensional cuprate superconductors is a linear function of temperature. Then the resistivity and scattering rate of CPs turn out to be also linear-in-temperature. The origin of linear-in-temperature resistivity and Planckian dissipation in cuprate superconductors may be the linear temperature dependence of the damping coefficient of CPs.

The Vacuum Energy Parameter (VEP) of standard cosmology denotes the fraction of the critical density attributed to the accelerated expansion of the Universe. Astrophysical evidence sets the numerical range of VEP at 0.692 +, - 0.012, yet the root cause of these values is currently unknown. Drawing from the stochastic interpretation of early-Universe cosmology, we develop here a derivation of the VEP based on classical diffusion theory and the Langevin equation. Predictions are shown to be in reasonable agreement with observations.

The relentless challenges posed by the COVID-19 pandemic have placed immense pressure on densely populated regions, such as Hong Kong, particularly in the context of an aging population. The burgeoning demand for healthcare services has necessitated the optimization of resource al-location within healthcare systems. This research presents an innovative analytical queueing model specifically tailored to address the intricate gridlock dynamics prevalent in healthcare systems, with a focused application in the challenging landscape of COVID-19 management. Hong Kong, like numerous other regions, has experienced a demo-graphic transition leading to heightened healthcare demands. Traditional methodologies for scrutinizing gridlock issues in healthcare settings, including simulation models, often prove suboptimal within the context of sophisticated optimization frameworks. This study, however, delves into analytical queueing models, offering a more efficient and effective approach. These models have been underexplored, primarily due to the intricacies involved in modeling gridlock propagation within restricted capacity, a critical concern magnified by the COVID-19 pandemic. The novel Queueing Model introduced in this research incorporates the maintenance of queue capacities and topological aspects as exogenous parameters, rendering it exceptionally adaptable for diverse gridlock-related applications. This adaptability is especially pertinent to healthcare management during the COVID-19 crisis. The model meticulously integrates the concept of gridlock, furnishing a comprehensive comprehension of its origins and repercussions. Validation of this model involves rigorous comparison with established methodologies, exact results, and simulation data, consistently high-lighting its versatility and precision. By applying the model to healthcare management, we scrutinize the dynamics of patient flow within a Hong Kong hospital, placing specific emphasis on the critical issue of bed blockage, further exacerbated by the COVID-19 pandemic. The model's gridlock decomposition unearths the nuanced impact of bed blockage on distinct hospital units, addressing the imperative requirement for quantifying in-patient bed blockage within healthcare systems during a pandemic.

The intercomparison among different atmospheric monitoring systems is key for instrument calibration and validation. Common cases involve satellites, radiosonde, radio occultation (e.g. GNSS-RO) and atmospheric model outputs (e.g. ERA5 reanalysis). Since instruments and/or measures are not perfectly collocated, the miss-collocation uncertainty must be considered in the related intercomparison uncertainty budgets.

Studies have shown that in many practical applications, data interpolation by splines leads to better approximation and higher computational efficiency as compared to data interpolation by a single polynomial. Data interpolation by splines can be significantly improved if knots are allowed to be free rather than at a priori fixed locations such as data points. In practical applications, the smallest possible curvature is often desired. Therefore, optimal splines are determined by minimizing a derivative of continuously differentiable functions comprising the spline of the required order. The problem of obtaining an optimal spline is tantamount to minimizing derivatives of a nonlinear differentiable function over a Banach space on a compact set. While the problem of data interpolation by quadratic splines has been accomplished analytically, interpolation by splines of higher orders or in higher dimensions is challenging. In this paper, to overcome difficulties associated with the complexity of the interpolation problem, the interval over which data points are defined, is discretized and continuous derivatives are substituted by their discrete counterparts. It is shown that as the mesh of the discretization approaches zero, a resulting near-optimal spline approaches an optimal spline. Splines with the desired accuracy can be obtained by choosing an appropriate mesh of the discretization. By using cubic splines as an example, numerical results demonstrate that the linear programming (LP) formulation, resulting from the discretization of the interpolation problem, can be solved by linear solvers with high computational efficiency and resulting splines provide a good approximation to the optimal splines.

In this work a new hydroelectric basin modelling approach is described and applied to the Pontecosi basin, Italy. Several types of data sources were used to learn the model: a number of weather stations, satellite observations, Reanalysis dataset and basin data. With the goal of predicting the water level of the basin, the model was composed by three cascade modules. Firstly, different spatial interpolation methods, such as Kriging, Radial Basis Function and Natural Neighbours, were compared and applied to interpolate the weather stations data nearby the basin area to infer the main environmental variables (air temperature, air humidity, precipitation and wind speed) in the basin area. Then, using these variables as inputs, a neural network was trained to predict the mean soil moisture concentration over the area, also to improve the low availability due to satellite orbits. Finally, a non-linear auto regressive exogenous input (NARX) model was trained to simulate the basin level with different prediction horizons, using the data from the previous modules and past basin data (water level, discharge flow rate, turbine flow rate). Accurate predictions of the basin water level were achieved within 1 to 6 hours ahead, with mean absolute errors (MAE) between 2cm and 10cm respectively.

Prior studies have focused on the overall behavior of randomly moving particle swarms. However, the characteristics of the stochastic-constrained particles that form ubiquitously within these swarms remain oblivious. This study demonstrates a generalized diffusion equation for stochastic-constrained particles that considers the velocity and location aggregation effects observed from their parent particle swarm (i.e., a completely random particle swarm). This equation can be approximated as the form of Schr\"odinger equation in the microcosmic case (low relative density) and describe the dynamics of the total mass distribution in the macrocosmic case (high relative density). The predicted density distribution of the particle swarm in the stable aggregation state is consistent with the total mass distribution of massive, relaxed galaxy clusters (at least in the range of r < r_s), preventing cuspy problems in the empirical Navarro-Frenk-White (NFW) profile. This study opens a window to observe the dynamics of stochastic-constrained particles from a third perspective, from which the aggregation effect of particles without gravitation can be saw.

We apply the stochastic thermodynamics formalism to describe the dynamics of systems of complex Langevin and Fokker-Planck equations. We provide in particular a simple and general recipe to calculate thermodynamical currents, dissipated and propagating heat for networks of nonlinear oscillators. By using the Hodge decomposition of thermodynamical forces and fluxes, we derive a formula for entropy production that generalises the notion of non-potential forces and makes transparent the breaking of detailed balance and of time reversal symmetry for states arbitrarily far from equilibrium. Our formalism is then applied to describe the off-equilibrium thermodynamics of a few examples, notably a continuum ferromagnet, a network of classical spin-oscillators and the Frenkel-Kontorova model of nano friction.

Capsule endoscopy, which uses a wireless camera to take images of the digestive tract, is emerging as an alternative to traditional colonoscopy. The diagnostic values of these images depend on the quality of revealed underlying tissue surfaces. In this paper, we consider the problem of enhancing the visibility of detail and shadowed tissue surfaces for capsule endoscopy images. Using concentric circles at each pixel for random walks combined with stochastic sampling, the proposed method enhances the details of vessel and tissue surfaces. The framework decomposes the image into two detail layers that contain shadowed tissue surfaces and detail features. The target pixel value is recalculated for the smooth layer using similarity of the target pixel to neighboring pixels by weighting against the total gradient variation and intensity differences. In order to evaluate the diagnostic image quality of the proposed method, we used clinical subjective evaluation with a rank order on selected KID image database and compared to state of the art enhancement methods. The result showed that the proposed method provides a better result in terms of diagnostic image quality and objective quality contrast metrics and structural similarity index.

Several methods have been tried to estimate air temperature using satellite imagery. In this paper, the results of two machine learning algorithms, Support Vector Machine and Random Forest, are compared with Multivariate Linear Regression, TVX and Ordinary kriging. Several geographic, remote sensing and time variables are used as predictors. The validation is carried out using four different statistics on a daily basis allowing the use of ANOVA to compare the results. The main conclusion is that Random Forest with residual kriging produces the best results (R$^2$=0.612 $\pm$ 0.019, NSE=0.578 $\pm$ 0.025, RMSE=1.068 $\pm$ 0.027, PBIAS=-0.172 $\pm$ 0.046), whereas TVX produces the least accurate results. The environmental conditions in the study area are not really suited to TVX, moreover this method only takes into account satellite data. On the other hand, regression methods (Support Vector Machine, Random Forest and Multivariate Linear Regression) use several parameters that are easily calculated from a Digital Elevation Model, adding very little difficulty to the use of satellite data alone. The most important variables in the Random Forest Model were satellite temperature, potential irradiation and cdayt, a cosine transformation of the julian day.

Grid refinement is used to save computing costs while preserving high precision of fluid simulation. In the Lattice Boltzmann Method, grid refinement often uses interpolated values. We found a method, where interpolation in space and time is not required. For this purpose, we used the moment matching condition and rescaling the non-equilibrium part of populations, and developed the recalibration procedure that allows to transfer information between different LBM stencils in the simulation domain. Then, we build non-uniform lattice which uses stencils with different shapes on the transition. The resulting procedure is verified on the 2D Poisselle flow and the advected vortex benchmark. In prospect, grids with adaptive geometry can be built with the use of the proposed method.

The main sources of intrinsic noise in excitable cells at the microcircuit and network levels are the stochastic characteristics of ion channel gating and activation of the synaptic conductance. Studies using in vivo, in vitro, and in silico methods to examine the effects of synaptic background activity were not adequately investigated in non-neuronal excitatory cells, where neurotransmitter-based innervation also occurs. We created a mathematical model to replicate the background synaptic noise dynamics in a non-neuronal cell. We utilized the stochastic Ornstein-Uhlenbeck process to represent excitatory synaptic conductance, which was incorporated into a whole-cell model to produce spontaneous and evoked cellular electrical activities. The single-cell model includes many biophysically detailed ion channels represented by a set of ordinary differential equations in Hodgkin-Huxley and Markov formalisms. This paradigm effectively induced irregular spontaneous depolarizations (SDs) and spontaneous action potential (sAP) resembling in vitro-like electrical activity in the cells. The input resistance decreased by multiple factors, and the spontaneous action potential firing rate elevated. The potential to reach the action potential threshold is altered. Background synaptic activity can alter the input/output characteristics of non-neuronal excitatory cells. Suppressing these baseline activities would facilitate the discovery of new pharmaceutical targets for different clinical diseases.

The exact, non perturbative, response theory developed within the field of nonequilibrium molecular dynamics, also known as TTCF (Transient Time Correlation Function) applies to quite general dynamical systems. Its key element is called dissipation function, because it represents the power dissipated by external fields acting on the particle system of interest, whose coupling with the environment is given by deterministic thermostats. This theory has been initially developed for time independent external perturbations, and then it has been extended to time dependent perturbations. It has also been applied to dynamical systems of different nature, and to oscillator models undergoing phase transitions, which cannot be treated withe.g. linear response theory. The present work includes time dependent stochastic perturbations, in the theory, using the Karhunen-Loeve theorem. This leads to three different investigations of a given process. One in which a single realization of the stochastic coefficients is fixed, and averages are taken only over the initial conditions, as in a deterministic process. In the second, the initial condition is fixed, and averages are taken with respect to the distribution of stochastic coefficients. In the last investigations, one averages over both initial conditions and stochastic coefficients. We conclude illustrating the applicability of the resulting exact response theory with simple examples.

This study investigates the stochastic dynamics of hepatitis B virus (HBV) infection using a newly proposed stochastic model. Contrary to deterministic models that fail to encapsulate the inherent randomness and fluctuations in biological processes, our stochastic model provides a more realistic representation of HBV infection dynamics. It incorporates random variability, thereby acknowledging the changes in viral and cellular populations and uncertainties in parameters such as infection rates and immune responses. We examine the solution's existence, uniqueness, and positivity for the proposed model, followed by a comprehensive stability analysis. We provide the necessary and sufficient conditions for local and global stability, offering deep insight into the infection dynamics. Furthermore, we utilize numerical simulations to corroborate our theoretical results. This research provides a robust tool for understanding the complex behavior of HBV dynamics, contributing significantly to the ongoing quest for more effective HBV control and prevention strategies.

An exploding demand for processing capabilities related to the emergence of the IoT, AI and big data, has led to the quest for increasingly efficient ways to expeditiously process the rapidly increasing amount of data. These ways include different approaches like improved devices capable of going further in the more Moore path, but also new devices and architectures capable of going beyond Moore and getting more than Moore. Among the solutions being proposed, Stochastic Computing has positioned itself as a very reasonable alternative for low-power, low-area, low-speed, and adjustable precision calculations; four key-points beneficial to edge computing. On the other hand, chaotic circuits and systems appear to be an attractive solution for (low-power, green) secure data transmission in the frame of edge computing and IoT in general. Classical implementations of this class of circuits require intensive and precise calculations. This paper discusses the use of the SC framework for the implementation of nonlinear systems, showing that it can provide results comparable to those of classical integration, with much simpler hardware, paving the way for relevant applications.

The paper presents the author’ original modification of the residual method in order to construct an efficient tool for assessing the condition of photovoltaic cells, in particular – panels made of amorphous cells which are characterized by significant random operation parameter changes. The solution constructed has been verified through experiments, which show the usefulness of the proposed methodology as it allows to minimize the risk of obtaining false diagnostic information in the process of assessing the condition of photovoltaic panels.

Fuzzy Rule Interpolation (FRI) offers an innovative methodology for inferring outcomes from data points that match no established rules in a sparse fuzzy rule base. Generally, a traditional fuzzy rule-based system suffers from the inability to perform inference while facing an unmatched observation, given a sparse fuzzy rule base that is derived from insufficient data or inadequate human experience and that fails to cover parts of a specific problem domain. Fortunately, starting with a sparse rule base, Fuzzy Rule Interpolation (FRI) can help to formulate interpolated rules specifically in instances where observations fail to activate any existing rule. An interpolated rule gives a result called consequent, which is related to the unmatched observation so that the inference system will not fail to generate an outcome corresponding to the observation as the output. Conventionally, these interpolated rules will be discarded after obtaining outcomes. However, valuable information may be potentially embedded within those discarded rules addressing the limitation regarding the coverage of the original sparse knowledge space. Moreover, without exploiting the potential value of the interpolated rules, the fuzzy inference system based on a sparse rule base will never improve its efficiency and robustness in the long run. Therefore, this underscores the necessity for Dynamic Fuzzy Rule Interpolation (D-FRI), which aims to collate interpolated rules that cover particular knowledge spaces absent in the sparse rule base. Intuitively, those frequently appeared interpolated rules are considered valuable because future observations are more likely to hit the respective problem areas. Promoting selected interpolated rules to the sparse rule base can enhance the overall coverage of the knowledge space and increase the inference efficiency over time. This paper introduces a framework for Dynamic Fuzzy Rule Interpolation that integrates the widely used Transformation-based Fuzzy Rule Interpolation (T-FRI) with an effective form of clustering of interpolated results. Experimental findings validate the effectiveness of this approach.

The paper's main goal is to introduce a novel type of interpolative convex contraction and build up some fresh findings for it using the interpolative convex contractions' progressive approach. We have established certain conclusions using orbitally S-complete and Suzuki type contractions in F-metric spaces. My research aims to examine the fixed point method's efficacy in solving fractional differential equation with boundary conditions.

The fitting and modeling of skewed, complex, symmetric, and asymmetric datasets is an exciting research topic in many fields of applied sciences, notably lifetime, medical, and financial sciences. The Heavy tailed Nadarajah Haghighi model is introduced by compounding the heavy tailed family and Nadarajah Haghighi distribution in this paper. The so obtained model has three parameters accounting for the scale and shape of the distribution. The proposed distribution's fundamental characteristics such as probability density, cumulative distribution, hazard rate, and survival functions are provided, and several key statistical properties are established, and several entropy information measures are proposed. The estimation of model parameter is performed via the ML procedure. Further, different simulation experiments are conducted to demonstrate the performance of proposed estimator’s using some measure, like the average estimate, the average bias, and associated mean square error. Finally, we are applied our proposed model to analyze three different real datasets. In our illustration, we are compared the practicality of the recommended model with several well-known competing models.

To solve the problem that sampling-based Rapidly-exploring Random Tree (RRT) method is difficult to guarantee optimality. This paper proposed the Post Triangular Processing of Midpoint Interpolation method minimized the planning time and shorter path length of the sampling-based algorithm. The proposed Post Triangular Processing of Midpoint Interpolation method makes a closer to the optimal path and somewhat solves the sharp path problem through the interpolation process. The experiments were conducted to verify the performance of the proposed method. Applying the method proposed in this paper to the RRT algorithm increases the efficiency of optimization compared to the planning time.

Real assets in the energy market are subject to ecological uncertainty due to the penetration of renewables. We illustrate this point by analysing electrolysers, a class of assets which recently became the subject of a large interest, as they lead to the production of the desirable green hydrogen and green ammonia. The latter has the advantage of being easily stored and has huge potential in decarbonising both the fertilizer and shipping industries. We consider the optimization of green ammonia production with different types of electricity procurement in the context of stochastic power and ammonia markets, a necessary assumption to translate the features of renewable, hence intermittent, electricity. We emphasise the importance of using stochastic prices to effectively model the volatile nature of the price dynamics, illustrating the project risks that can be mitigated by hedging activities. This study shows the pivotal role of flexibility when dealing with fluctuating renewable production and volatile electricity prices to maximise profits and better manage risks.

Dealing with randomness is a crucial aspect that cost-effectiveness analysis (CEA) tools need to address, but existing stochastic CEA tools have rarely examined risk and return from the perspective of global benefits. This paper proposes a stochastic CEA tool that supports medical decision-making from the perspective of global benefits of risk and return, the risk-adjusted incremental cost-effectiveness ratio (ICER). The tool has a traditional form of ICER but uses the risk-adjusted expected cost. Theoretically, we prove that the tool can provide marginal medical decisions that promote the risk-return level on global benefits within any intervention structure and can also serve as a discriminating condition for the optimal intervention structure. Numerical simulations within a framework of mean-variance support the conclusions in this paper.

We outline and test a new methodology for genuine simulation of stochastic processes with any dependence and any marginal distribution. We reproduce time dependence with a generalized, time symmetric or asymmetric, moving-average scheme. This implements linear filtering of non-Gaussian white noise, with the weights of the filter determined by analytical equations in terms of the autocovariance of the process. We approximate the marginal distribution of the process, irrespective of its type, using a number of its cumulants, which in turn determine the cumulants of white noise in a manner that can readily support the generation of random numbers from that approximation, so that it be applicable for stochastic simulation. The simulation method is genuine as it uses the process of interest directly without any transformation (e.g. normalization). We illustrate the method in a number of synthetic and real-world applications with either persistence or antipersistence, and with non-Gaussian marginal distributions that are bounded, thus making the problem more demanding. These include distributions bounded from both sides, such as uniform, and bounded form below, such as exponential and Pareto, possibly having a discontinuity at the origin (intermittence). All examples studied show the satisfactory performance of the method.

In line with Pomeau’s conjecture about the relevance of directed percolation (DP) to turbulence onset/decay in wall-bounded flows, we propose a minimal stochastic model dedicated to the interpretation of the spatially intermittent regimes observed in channel flow before its return to laminar flow. Numerical simulations show that a regime with bands obliquely drifting in two stream-wise symmetrical directions bifurcates into an asymmetrical regime, before ultimately decaying to laminar flow. The model is expressed in terms of a probabilistic cellular automaton evolving von Neumann neighbourhoods with probabilities educed from a close examination of simulation results. It implements band propagation and the two main local processes: longitudinal splitting involving bands with the same orientation, and transversal splitting giving birth to a daughter band with orientation opposite to that of its mother. The ultimate decay stage observed to display one-dimensional DP properties in a two-dimensional geometry is interpreted as resulting from the irrelevance of lateral spreading in the single-orientation regime. The model also reproduces the bifurcation restoring the symmetry upon variation of the probability attached to transversal splitting, which opens the way to a study of the critical properties of that bifurcation, in analogy with thermodynamic phase transitions.

Partial differential representation of self-similarity feature has been derived from notion of the homogenous function in general sense. This representation allows consideration of stochastic self-similar systems. As well as the partial differential representation allows consideration of Stochastic Partial Differential Equation.

Global pose refinement is a significant challenge within Simultaneous Localization and Mapping (SLAM) frameworks. For LIDAR-based SLAM systems, pose refinement is integral to correcting drift caused by the successive registration of 3D point clouds collected by the sensor. A divergence between the actual and calculated platform path characterizes this error. In response to this challenge, we propose a linear, parameter-free model, which uses a closed circuit for global trajectory corrections. Our model maps rotations to quaternions and uses Spherical Linear Interpolation (SLERP) for transitions between them. The intervals are established from the constraint set by the Least Squares (LS) method on rotation closure and are proportional to the circuit's size. Translations are globally adjusted in a distinct linear phase. Additionally, we suggest a coarse-to-fine pairwise registration method, integrating Fast Global Registration and Generalized ICP with multiscale sampling and filtering. The proposed approach is tested on three varied datasets of point clouds, including Mobile Laser Scanners and Terrestrial Laser Scanners. These diverse datasets affirm the model's effectiveness in 3D pose estimation with substantial pose differences and efficient pose optimization in larger circuits.

Obtaining classification confidence at the pixel level is a challenging task for accuracy assessment in remote sensing image classification. Among the various methods for estimating classification confidence at the pixel level, interpolation-based methods have drawn special attention in the literature. Even though they have been widely recognized in the literature, their usefulness has not been rigorously evaluated. This paper conducts a comprehensive evaluation of three interpolation-based methods: local error matrix method, bootstrap method, and geostatistical method. We applied each of the three methods to three representative datasets with different spatial resolutions, spectral bands, and the number of classes. We then derive the estimated classification confidence and true classification confidence and compared the results with each other using both exploratory data analysis (bi-histogram) and statistical analysis (Willmott's d and Binned classification quality). The results indicate that the three interpolation methods provide some interesting insights on various aspects of estimating per-pixel classification confidence. Unfortunately, the interpolation assumes that classification confidence is smooth across the space, which is usually not true in practice. In other words, interpolation-based methods have limited practical use.

As a model for economic and ecological systems, replicator dynamics represents a basic form of agent competition for finite resources. Here, we investigate the effects of stochastic resetting in this kind of processes. Random reset events abruptly lead individual resources to a small value from which dynamics must start anew. Numerical results show that resource distribution over the population of competing agents develops highly nonuniform profiles, exhibiting clustering and fluctuations with anomalous dependence on the population size. This non-standard statistical behavior jeopardizes an analytical treatment based on mean-field assumptions. We propose alternative simplified analytical approaches which provide a stylized description of entropy evolution for the clustered distribution of resources and explain the unusually slow decrease of fluctuations.

Intratumoral heterogeneity can exist along multiple axes: Cancer Stem Cells (CSCs)/non-CSCs, drug-sensitive/drug-tolerant states and a spectrum of epithelial-hybrid-mesenchymal phenotypes. Further, these diverse cell-states can switch reversibly among one another, thereby posing a major challenge to therapeutic efficacy. Therefore, understanding the origins of phenotypic plasticity and heterogeneity remains an active area of investigation. While genomic components (mutations, chromosomal instability) driving heterogeneity have been well-studied, recent reports highlight the role of non-genetic mechanisms in enabling both phenotypic plasticity and heterogeneity. Here, we discuss various processes underlying phenotypic plasticity such as stochastic gene expression, chromatin reprogramming, asymmetric cell division and the presence of multiple “attractors”. These processes can facilitate a dynamically evolving cell population such that a subpopulation of (drug-tolerant) cells can survive lethal drug exposure and recapitulate population heterogeneity on drug withdrawal, leading to relapse. These drug-tolerant cells can be both pre-existing and also induced by the drug itself through cell-state reprogramming. The dynamics of cell-state transitions both in absence and presence of the drug can be quantified through mathematical models. Such a dynamical systems approach to elucidating patterns of intratumoral heterogeneity by integrating longitudinal experimental data with mathematical models can help design effective combinatorial and/or sequential therapies for better clinical outcomes.

In this paper, we compare the performance between systems of ordinary and (Caputo) fractional differential equations depicting the susceptible-exposed-infectious-recovered (SEIR) models of diseases. In order to understand the origins of both approaches as mean-field approximations of integer and fractional stochastic processes, we introduce the fractional differential equations as approximations of some type of fractional nonlinear birth--death processes. Then, we examine validity of the two approaches against empirical courses of epidemics; we fit both of them to case counts of three measles epidemics that occurred during the pre-vaccination era in three different locations. While FDEs appear more flexible in fitting empirical data, our ODEs offered better fits to two out of three data sets. Important differences in transient dynamics between these modeling approaches are discussed.

Healthcare leaders are faced with many financial challenges in the contemporary environment, leading to many bankruptcies in recent years. What is not well understood are the specific conditions that lead to organizational economic failure. Ideally, hospital leaders could examine the literature for guidance to predict bankruptcy before it occurs. However, existing regression methods have proven inadequate, especially when the finance variables follow a nonnormal frequency pattern. Furthermore, the regression approach encounters difficulties due to multicollinearity and nonnormality issues. Therefore, an alternate stochastic approach for predicting the probability of hospital bankruptcy is warranted. This new method includes computing and interpreting the correlations between the hospital’s revenue and expenses for bivariate lognormal data, calculating the risk of bankruptcy due to the mismatch between revenues and expenses, determining the likelihood of a hospital’s expense exceeding the state’s median expense level, and evaluating a hospital’s financial memory level. Ultimately in our research, we demonstrate a novel approach to anticipate hospital bankruptcy, which can be particularly useful for health policy analysts and hospital leaders in making informed financial decisions and proactively managing risks to ensure the sustainability and stability of their institutions.

The present paper reports an alternative solution of the measurement problem in quantum theory. The measurement problem can be characterized as that the U-procedure verses the R-procedure. The R-procedure is tested by the Yes/No measurement originally proposed by von Neuman and discussed in details by Penrose. We propose a novel stochastic sampling method to tackle the measurement paradox. Each testing sample produces a pair of Yes-number c and No-number d, which in turn generates a sample phase with respect to the exponential form of (c+id). All the sample phases forms a group G_R. Taking G_R as the sampling potential, combined with the dynamic phase group G_U，we establish a unified model of U-procedure and R-procedure for the evolution of wavefunction. Based on the Born rule, the probability of a sample is given by the squared magnitude of （c+id）. Some metaproperties of the new U-R model, such as natural transformation, consistency, and completeness, are presented. The present work provides a new picture of quantum mechanics and alike. The measurement problem widely exists in dynamic analysis, higher cognition, and large data modeling. Thus, the work reported here has a wide range of applications. Keywords: The measurement problem; U-R procedures; the Yes/No type measurement; stochastic sampling; sample phase; Born probability.

It is known that the source of a tectonic earthquake in the framework of the theory of elasticity and viscoelasticity is considered as a displacement along a certain fault surface. Usually, when describing the source, the geometry of the fault surface is simplified to a flat rectangular area. The displacement vector is assumed to be constant. In this paper, we propose a model of an earthquake source in the form of a displacement with a constant vector along a stochastic surface. A number of standard assumptions were made when modeling. We take into account only the elastic properties of the medium. We consider the Earth’s crust as a half-space and assume that the medium is homogeneous and isotropic. For the mathematical description of the earthquake source, we use the classical force equivalent of displacement along the fault. This is the distribution of double pairs of forces. The field of displacements under the action of body forces is found through a combination of Mindlin nuclei of strain. The paper presents solutions for strike-slip fault. To obtain a stochastic fault surface, we propose to introduce a random deformation of a rectangular flat surface. The paper presents the results of a computational experiment comparing the levels and regions of relative deformations of the Earth’s crust in the case of displacement along a flat rupture surface and along a stochastic one. In the case of a stochastic surface, the regions of relative deformations become asymmetric. The developed model will help to indirectly take into account various mechanical properties of the Earth’s crust when modeling deformations caused by earthquakes.

In this paper, we first establish upper stochastic bounds on the lifetime of a used cold standby system with arbitrary age, using the likelihood ratio order and the usual stochastic order. Then, stochastic comparisons are made between the lifetime of a used cold standby system with age t and the lifetime of a cold standby system consisting of used components with age t using the likelihood ratio order and the usual stochastic order. We use illustrative examples to explore the results presented.

A stochastic dominance (SD) relation can be defined by two different perspectives: One from the view of distributions, and the other one from the view of expected utilities. In early days, Fishburn investigated SD from the view of distributions and we refer this perspective as Fishburn’s SD. One of his many results was the development of fractional-order SD for continuous distributions. However, discrete fractional-order SD may not be generalized directly since some properties of fractional calculus do not have a discrete counterpart. In this paper, we develop a discrete analogue of fractional-order SD from the view of distributions. We generalize the order of SD by Lizama’s fractional delta operator, show the preservation of SD hierarchy, and formulate the utility classes that are congruent with our SD relations. This work brings a message that some results of discrete SD cannot be generalized directly from continuous SD. We characterize the difference between discrete and continuous fractional-order SD, as well as the way to handle them.

Due to the increased use of memristors, and its many applications, the use of emulators has grown in parallel to avoid some of the difficulties presented by real devices such as variability and reliability. In this paper, we present a memristive emulator designed using a Switched Capacitor (SC), this is, an analog component/ block and a control part or block implemented using stochastic computing (SCo) and therefore fully digital. Our design is thus a mixed signal circuit. Memristor equations are implemented using stochastic computing to generate the control signals necessary to work with the controllable resistor implemented as switched capacitor.

Cellular Nonlinear Networks (CNN) are a concept introduced in 1988 by Leon Chua and Lin Yang as a bio-inspired architecture, capable of massively parallel computation. Later on, CNN have been enhanced by incorporating designs that incorporate memristors to profit from their processing and memory capabilities. In addition, Stochastic Computing (SC) can be used to optimize the quantity of required processing elements; thus it provides a lightweight approximate computing framework, quite accurate and effective, though. In this work, we propose utilization of SC in designing and implementing a memristor-based CNN. As a proof of the proposed concept, an example of application is presented. This application combines Matlab and a FPGA in order to create the CNN. The implemented CNN has then been used to perform three different real-time applications on a 512x512 gray-scale and a 768x512 color image: storage of the image, edge detection, and image sharpening. It has to be pointed out that the same CNN has been used for the three different tasks, with the sole change of some programmable parameters. Results show an excellent capability with significant accompanying advantages, like the low number of needed elements further allowing for a low cost FPGA-based system implementation, something confirming the system’s ability for real time operation.

We present a class of efficient parametric closure models for 1D stochastic Burgers equations. Casting it as statistical learning of the flow map, we derive the parametric form by representing the unresolved high wavenumber Fourier modes as functionals of the resolved variables’ trajectory. The reduced models are nonlinear autoregression (NAR) time series models, with coefficients estimated from data by least squares. The NAR models can accurately reproduce the energy spectrum, the invariant densities, and the autocorrelations. Taking advantage of the simplicity of the NAR models, we investigate maximal and optimal space-time reduction. Reduction in space dimension is unlimited, and NAR models with two Fourier modes can perform well. The NAR model’s stability limits time reduction, with a maximal time step smaller than that of the K-mode Galerkin system. We report a potential criterion for optimal space-time reduction: the NAR models achieve minimal relative error in the energy spectrum at the time step where the K-mode Galerkin system’s mean CFL number agrees with the full model’s.

In recent years, the field of deep learning achieved considerable success in pattern recognition, image segmentation and may other classification fields. There are a lot of studies and practical applications of deep learning on images, video or text classification. In this study, we suggest a method for changing the architecture of the most performing CNN models with the aim of designing new models to be used as stand-alone networks or as a component of an ensemble. We propose to replace each activation layer of a CNN (usually a ReLu layer) by a different activation function stochastically drawn from a set of activation functions: in this way the resulting CNN has a different set of activation function layers.

Regression models form the core of the discipline of econometrics. One of the basic assumptions of classical linear regression model is that the values of the explanatory variables are fixed in repeated sampling. However, in most of the real life cases, particularly in economics the assumption of fixed regressors is not always tenable. Under a non-experimental or uncontrolled environment, the dependent variable is often under the influence of explanatory variables that are stochastic in nature. There is a huge literature related to stochastic regressors in various aspects. In this paper, a historical perspective on some of the works related to stochastic regressor is being tried to pen down based on literature search.

The objective of this study is to evaluate the effects of the 3D point cloud density derived from unmanned aerial vehicle (UAV) photogrammetry and structure from motion (SfM) and multi-view stereopsis (MVS) techniques, the interpolation method for generating a digital terrain model (DTM), and the resolution (grid size) of the derived DTM on the accuracy of estimated heights in small areas, where a very accurate high spatial resolution is required. A UAV-photogrammetry project was carried out on a bare soil of 13 × 13 m with a rotatory wing UAV at 10 m flight altitude (equivalent ground sample distance = 0.4 cm). The 3D point cloud was derived, and five sample replications representing 1, 2, 3, 4, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80 and 90% of the original cloud were extracted to analyze the effect of cloud density on DTM accuracy. For each of these samples, DTMs were derived using four different interpolation methods (Inverse Distance Weighted (IDW), Multiquadric Radial Basis Function (MRBF), Kriging (KR), and Triangulation with Linear Interpolation (TLI)) and 15 DTM grid size (GS) values (20, 15, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0.67, 0.5, and 0.4 cm). Then, 675 DTMs were analyzed. The results showed, for each interpolation method and each density, an optimal GS value (most of the cases equal to 1 cm) for which the Root Mean Square Error (RMSE) is minimum. IDW was the interpolator which yielded best accuracies for all combination of densities and GS. Its RMSE, considering the raw cloud, was 1.054 cm. The RMSE increased 3% when a point cloud with 80% extracted from the raw cloud was used to generate the DTM. When the point cloud included the 40% of the raw cloud, RMSE increased 5%. For densities lower than 15%, RMSE increased exponentially (45% for 1% of raw cloud). The grid size minimizing RMSE for densities of 20% or higher was 1 cm, which represents 2.5 times the ground sample distance of the pictures used for developing the photogrammetry project.

Autism spectrum disorder (ASD) is characterized by diverse patterns of social interaction, communication, and behavior. This study explores the application of clustering and stochastic migration models to visualize and understand the social and neural connectivity patterns associated with autism. The clustering model, with its strong internal connections and limited external connections, mirrors the tight-knit social groups and communication barriers often observed in individuals with autism. In contrast, the stochastic model, characterized by more dispersed connections and greater adaptability, represents the variability and flexibility seen in social interactions across typical individuals. By comparing these models, we can gain deeper insights into the unique challenges and strengths of individuals with autism, highlighting the importance of personalized interventions. Additionally, these models provide a visual framework for understanding the differences in brain connectivity patterns observed in autism, with implications for both behavioral and neural correlates. This study underscores the value of computational models in advancing our understanding of autism and guiding effective support strategies.

Here are analysed data taken in two hydrocarbon fields ("A" and "B"), located in the western part of Sava Depression (North Croatia). They are in the secondary phase of production. The selected reservoirs "L" (in the “A” Field) and "K" (“B”) are of the Lower Pontian (Upper Miocene) age and belong to Kloštar-Ivanić Formation. Due to strong tectonics, there are numerous tectonic block, relatively rarely sampled with well and laboratory tests. Here are selected two variables for interpolation - reservoirs permeabilities and the injected volumes of field water. The following interpolation methods are described, compared and applied: Nearest Neighbourhood, Natural Neighbour (the first time in the Sava Depression) and Inverse Distance Weighting. The last one has been proven as the most appropriate for datasets with size lower than 20 points.

Analyzing the records of Advanced LIGO and Virgo gravitational observatories, we found a specific time-domain asymmetry inherent only to the signals of their gravitational detectors. Experiments with different periodic signals, Gaussian and non-Gaussian noises, made it possible to conclude that the noise of gravitational detectors is an unusual mixture of signals. We also developed a specialized Pearson correlation analyzer to recognize the gravitational-wave events. It turned out that the LIGO detectors’ output signals include a significant (– 6 dB) component, which has the properties of records of reliably recognized gravitational waves. It allows us to argue that the gravitational background noise is largely due to the processes of merging astronomical objects. Since the specific signal is registered by the detectors continuously, we can consider the sub-kilohertz band gravitational oscillations field as detected. Our analysis method also allows us to estimate the contribution of the gravitational background component to the total signal energy. With its help, it will be possible not only to provide the radio-frequency estimation of the magnitude of gravitational disturbances but also, possibly, to construct a map of the gravitational noise of the sky.

In this paper we present a successful implementation of the Hindmarsh–Rose model within a SC environment. The merits of the proposed approach are design simplicity, due to stochastic computing, and the ease of implementation. Simulation results showed that the approximation achieved is equivalent to introducing a noise source into the original model. A study for the level of noise introduced, according to the number of bits in the stochastic sequence, has been performed. Additionally, we demonstrate that such an approach, even though it is noisy, it reproduces the behaviour of biological systems, which are intrinsically noisy. It is also demonstrated that a speedup of x2 compared to biological systems is easily achievable, with a very small number of gates, thus paving the road for the in silico implementation of large neuron networks.

Semantic segmentation is a very popular topic in modern computer vision and it has applications to many fields. Researchers proposed a variety of architectures over time, but the most common ones exploit an encoder-decoder structure that aims to capture the semantics of the image and it low level features. The encoder uses convolutional layers, in general with a stride larger than one, to extract the features, while the decoder recreates the image by upsampling an using skip connections with the first layers. In this work, we use DeepLab as architecture to test the effectiveness of creating an ensemble of networks by randomly changing the activation functions inside the network multiple times. We also use different backbone networks in our DeepLab to validate our findings. We manage to reach a dice coefficient of 0.888, and a mean Intersection over Union (mIoU) of 0.825, in the competitive Kvasir-SEG dataset. Results in skin detection also confirm the performance of the proposed ensemble, which is ranked first with respect to other state-of-the-art approaches (including HardNet) in a large set of testing datasets. The developed code will be available at https://github.com/LorisNanni.

Water is essential to all life-forms including various ecological, geological, hydrological, and climatic processes/activities. With changing climate, associated El Nino/Southern Oscillation (ENSO) events appear to stimulate highly uncertain patterns of precipitation (P) and evapotranspiration (EV) processes across the globe. Changes in P and EV patterns are highly sensitive to temperature variation and thus also affecting natural streamflow processes. This paper presents a novel suite of stochastic modelling approaches for associating streamflow sequences with climatic trends. The present work is built upon a stochastic modelling framework HMM_GP that integrates a Hidden Markov Model with a Generalised Pareto distribution for simulating synthetic flow sequences. The GP distribution within HMM_GP model is aimed to improve the model's efficiency in effectively simulating extreme events. This paper further investigated the potentials of Generalised Extreme Value Distribution (EVD) coupled with an HMM model within a regression-based scheme for associating impacts of precipitation and evapotranspiration processes on streamflow. The statistical characteristic of the pioneering modelling schematic has been thoroughly assessed for their suitability to generate/predict synthetic river flows sequences for a set of future climatic projections. The new modelling schematic can be adapted for a range of applications in the area of hydrology, agriculture and climate change.

In this paper, we introduce the evolve-then-filter (EF) regularization method for reduced order modeling of convection-dominated stochastic systems. The standard Galerkin projection reduced order model (G-ROM) yield numerical oscillations in a convection-dominated regime. The evolve-then-filter reduced order model (EF-ROM) aims at the numerical stabilization of the standard G-ROM, which uses explicit ROM spatial filter to regularize various terms in the reduced order model (ROM). Our numerical results based on a stochastic Burgers equation with linear multiplicative noise. It shows that the EF-ROM is significantly better results than G-ROM.

This study purposes the use of plug-in electric vehicles for demand side management (DSM) considering uncertainties in demand as well as uncertainties due to mobility of PEV to mitigate peak demand. The solution also seeks to reduce electric cost in addition to reducing the effects of greenhouse gases. In recent years DSM using distributed storage system such as battery energy management system (BESS) and plugged-in electric vehicles (PEV) have become very prevalent with most implementations resorting to deterministic load forecast. These methods do not consider the potential growth in demand making their solutions less robust. In this study we propose a real-time density demand forecast and stochastic optimization for robust operation of PEV for a building. This method accounts for demand uncertainties in addition to uncertainties in mobile energy storage as found in PEV, making the resulting solution robust as compared to the deterministic case. A case study on a real site in South Korea is used for verification and testing. The proposed study is verified and tested against existing algorithms. The result verifies the effectiveness of the proposed approach

A Demand-side management technique are deployed along with battery energy-storage systems (BESSs) to lower the electricity cost by mitigating the peak load of a building. Most of the existing methods rely on manual operation of the BESS, or even an elaborate building energy-management system resorting to a deterministic method that is susceptible to unforeseen growth in demand. In this study we propose a real-time optimal operating strategy for BESS based on density demand forecast and stochastic optimization. This method takes into consideration uncertainties in demand when accounting for an optimal BESS schedule, making it robust compared to the deterministic case. The proposed method is verified and tested against existing algorithms. Data obtained from a real site in South Korea is used for verification and testing. The results show that the proposed method is effective, even for the cases where the forecasted demand deviates from the observed demand

The classical Lotka-Volterra predator-prey model is globally stable and uniformly persistent. However, in real-life biosystems extinction of species due to stochastic effects is possible. In this paper, we consider the classical Lotka-Volterra predator-prey model under stochastic perturbations. For this model, using an analytical technique based on the direct Lyapunov method and a development of the ideas of R.Z. Khasminskii, we find the precise sufficient conditions for the stochastic extinction of one and both species and the precise necessary conditions for the stochastic system persistence. The stochastic extinction occurs via a process known as the stabilization by noise of the Khasminskii type, and, in order to establish the sufficient conditions for extinction, we found the conditions for the stabilization. The analytical results are illustrated by numerical simulations. {\bf Keywords:} stochastic perturbations, white noise, Ito's stochastic differential equation, the Lyapunov functions method, stability in probability, stabilization by noise, stochastic extinction, persistence.

This study explores the potential trajectories of human civilization influenced by the development and competition of advanced artificial intelligence (AI) systems. Using a system of stochastic partial differential equations (SPDEs), we model the probability density function representing the state of civilization in a phase space defined by prosperity and knowledge. The model incorporates diffusion, growth, saturation, and drift terms, alongside stochastic noise to reflect the uncertainties and random fluctuations in AI impacts.The simulation results are visualized in a probability density graph, revealing the likelihood of various outcomes ranging from extinction and regression, to prosperity and technological advancement. The analysis highlights the balanced prospects of human future, with significant probabilities in both positive and negative directions. Positive scenarios suggest potential for increased prosperity and knowledge, emphasizing the importance of effective AI management and international cooperation. Conversely, the notable risks of regression and extinction underline the need for strategic interventions to mitigate adverse impacts. Our findings stress the stochastic nature of future developments and the critical role of adaptive and flexible policies in steering human civilization towards favorable outcomes. This study provides a simple yet comprehensive framework for understanding the complex dynamics at play and underscores the importance of proactive strategies in the age of AI.

The container relocation problem (CRP) is one of the important topics in container terminals, which aims to find the optimal sequence of movements to minimize the total number of relocations needed to retrieve all the containers within a bay according to a given specific layout and retrieval order. However, determining such an order may be difficult due to uncertainty. Various CRP solution methods proposed so far show a common pitfall concerning the complexity of the solving process, whose scale grows exponentially with the number of containers. Considering the retrieval operation process of the inbound containers in a bay consists of several operation rounds, this paper focuses on the real-time operation scenario of each round to minimize the actual number of relocations during the pick-up process, i.e., starting from the initial layout of the bay, we minimize the sum of the incremental number of blocking containers due to the layout change and the number of reallocated containers unrelated to the target container. We present a heuristic algorithm (called SPFH) based on the three rules (LL, MSS, and FSS), considering the adjustment of the pick-up order of containers. Numerical experiments on a set of instances generated from the literature are performed. The proposed methodology excels in both execution efficiency and solution quality and is capable of solving instances that were impossible to solve before. This also makes our methods appealing for being applied in practice. The source codes are available at https://github.com/zhou-sf/rt-scrp.

Iterative Stochastic decoding is an alternative to standard fixed-point decoding of Low-density parity check codes (LDPC) to reduce inter-node routing. In this paper, we propose a Flexible Field-Programmable Gate Array (FPGA)-Based Stochastic decoding (SD) hardware architecture for LDPC codes in the Fifth-Generation (5G) New Radio (NR) Standard that supports decoding of set of various code rates. This decoder’s runtime flexibility is desirable to switch a better performing code rate automatically based on the channel conditions without the extra time needed for reprogramming of FPGA. An offline design method is implemented to generate the hardware description language (HDL) code description of the decoder for the required code-rate set, which is further synthesized and integrated into the Xilinx Kintex-7 series FPGA board to determine the hardware resource utilisation and processing throughput. The Synopsys design tools were employed during both the simulation and synthesis stages, in combination with TSMC 65-nm CMOS standard cell technology, to facilitate comparative analysis. Compared with state-of-the-art designs, the proposed architecture reduces hardware utilization by up to 26% and achieved energy efficiency by 52%.

The need to find the global minimum in multi-dimensional functions is a critical problem in many fields of science and technology. Effectively solving this problem requires the creation of initial solution estimates, which are subsequently used by the optimization algorithm to search for the best solution in the solution space. In the context of this article, a novel approach to generating the initial solution distribution is presented which is applied to a genetic optimization algorithm. Using the k-means clustering algorithm, a distribution based on data similarity is created. This helps in generating initial estimates that may be more tailored to the problem. Additionally, the proposed method employs a rejection sampling algorithm to discard samples that do not yield better solution estimates in the optimization process. This allows the algorithm to focus on potentially optimal solutions, thus improving its performance. Finally, the article presents experimental results from the application of this approach to various optimization problems, providing the scientific community with a new method for addressing this significant problem.

Machinery parts gradually wear out over time due to regular usage. To improve machinery health and prevent critical issues, a reliable prognosis framework can be implemented by monitoring the behaviour of machinery parts and issuing warnings before they reach a critical state. To achieve this, vibration data from roller bearings experiencing various fault conditions have been collected. Different techniques from the literature were combined to analyze the distinct configurations in the vibration data sets and identify the main defects in roller bearings. The significant features extracted from this analysis were then used to create optimized stochastic model equations, separately regressing inner and outer race fault features to healthy bearing features under random conditions. These models can help engineers design more dependable systems, optimize their performance, and minimize the risk of failures and downtime.

In this paper we consider a failed cold standby systems and obtain some stochastic bounds on the idle time of such systems. It is assumed that, at a time t the failure of a cold standby system is observed, but the actual time of the failure is a time before t. The idle time of the inactive cold standby system is measured as the time between the exact time of failure, which is a random variable, and time t. We use the idle times of the components of the system to obtain lower and upper bounds on the lifetime of the standby system to which they are assembled. We bring some examples to show that the results and the conditions imposed to get them are achievable.

A thermal load dispatch problem minimizes the number of objectives viz operating cost and emission of gaseous pollutants together while allocating the power demand among the committed generating units subject to physical and technological system constraints. A stochastic thermal load dispatch problem is undertaken while taking into consideration, the uncertainties, errors in data measurements and nature of load demand which is random. Owing to uncertain load demand, variance due to mismatch of power demand termed as risk, is considered as another conflicting objective to be minimized. Generally multiobjective problems generate a set of non-inferior solutions are generated and supplied to a decision maker to select the best solution from the set of non-inferior solutions. This paper proposes opposition-based greedy heuristic search (OGHS) method to generate a set of non-inferior solutions. Opposition-based learning is applied to generate initial population to select good candidates. Migration to maintain diversity in the set of feasible solutions is also based on opposition-based learning. Mutation strategy is implemented by perturbing the genes heuristically in parallel and better one solution is sought for each member. Feasible solutions are achieved heuristically by modifying the generation-schedules in such a manner that violation of operating generation limits are avoided. The OGHS method is simple to implement and provides global solutions derived from the randomness of the population generated without tuning of parameters. Decision maker exploits fuzzy membership functions to decide the final decision. Validity of the method has been demonstrated by analysing systems in different scenarios consisting of six generators and forty generators.

It is quite easy to stochastically distort an original count variable to obtain a new count variable with relatively more variability than in the original variable. Many popular overdispersion models (variance greater than mean) can indeed be obtained by mixtures, compounding or randomlystopped sums. There is no analogous stochastic mechanism for the construction of underdispersed count variables (variance less than mean), starting from an original count distribution of interest. This work proposes a generic method to stochastically distort an original count variable to obtain a new count variable with relatively less variability than in the original variable. The proposed mechanism, termed condensation, attracts probability masses from the quantiles in the tails of the original distribution and redirect them toward quantiles around the expected value. If the original distribution can be simulated, then the simulation of variates from a condensed distribution is straightforward. Moreover, condensed distributions have a simple mean-parametrization, a characteristic useful in a count regression context. An application to the negative binomial distribution resulted in a distribution allowing under, equi and overdispersion. In addition to graphical insights, fields of applications of special cases of condensed Poisson and condensed negative binomial distributions were pointed out as an indication of the potential of condensation for a flexible analysis of count data

Dams design and their operation cause strong environmental alteration and therefore a long-term debate is ongoing for the scale of these projects. At the same time, the concept of Environmental Flow Assessment (EFA) is a crucial element of modified ecosystems featuring large infrastructure such as dams and reservoirs for mitigating potential environmental degradation while they operate. Nowadays, integrated scientific frameworks are required to quantify the risks caused by large infrastructure. Through the use of stochastic analysis, it is possible to quantify these uncertainties, and present a solution that incorporates long-term persistence and environmental sustainability into a balanced reservoir simulation model. In this work, an attempt is made to determine a benchmark reservoir size incorporating hydrological and ecological criteria though stochastic analysis. The primary goal is to ensure the best possible conditions for the ecosystem, and then secondarily to allow a steady supply of water for other uses. Using a synthetic timeseries based on historical inputs, it is possible to determine and preserve essential statistical characteristics of a river’s streamflow, and use these to detect the optimal reservoir capacity that maximizes environmental and local water demand reliability.

Interpolation is procedure that depends on spatial and/or statistical properties of analysed variable(s). It is special challenging task for data that included low number of samples, like dataset with less than 20 data. This problem is especially emphasized in the subsurface geological mapping, i.e. in the cases where data are taken solely from wells. Successful solutions of such mapping problems ask for knowledge about interpolation methods designed primarily for small datasets and dataset itself. Here are compared two methods, namely Inverse Distance Weighting and Modified Shepard’s Method, applied for three variables (porosity, permeability, thickness) measured in the Neogene sandstone hydrocarbon reservoirs (Northern Croatia). The results showed that pure cross-validation is not enough condition for appropriate map selection, but also geometrical features need to be considered, for datasets with less than 20 points.

Seasonal snow cover is a fundamental component of both the global energy budget and the water cycle. Its properties such as fractional snow cover or albedo are particularly affected by climate change. Several methods based on satellite data products are available to estimate these properties, each one with its pros and cons. This work presents a novel methodology that integrates three indexes applied to MODIS satellite data (Spectral Mixture Analysis (SMA), Normalized Difference Snow Index (NDSI) and Melt Area Detection Index (MADI)), to perform a spatio-temporal reconstruction of the fractional snow cover and albedo at 250 m spatial resolution in the Brunswick Peninsula, southwest Patagonia during the cold season (April-October) for the period 2000-2020. Three main steps are included: (1) the increase of the spatial resolution of MODIS (MOD09) data to 250 m using a spectral fusion technique; (2) the snow-cloud discrimination; (3) the daily spatio-temporal reconstruction of snow extent and its albedo signature with subpixel detection using the endmembers extraction and spectral mixture analysis. The results show a 98% agreement between MODIS snow detection and ground-based snow measurement at the automatic weather station Tres Morros (53.3174 °S, 71.2790 °W), with fractional snow cover values between 20% and 50%, showing a close relationship between snow and vegetation type. The number of snow days compiled from the MODIS data indicates a good performance (Pearson correlation of 0.9) compared with the number of skiing days at Cerro Mirador ski centre near Punta Arenas. Although the number of seasonal snow days show a significant increase trend of 0.54 days/year in Brunswick Peninsula during the 2000-2020 period a significant decreasing trend of -4.64 days/year was detected during 2010-2020, and also below the 400 m a.s.l. elevation, which is the most affected area. A reconstruction of the monthly mean temperature using the ERA5 Land reanalysis product shows a significant warming trend in May (0.068 ºC/year) and October (0.098 ºC/year) in the 2000-2020 period. Under the future emission scenario RCP8.5, the regional climate model RegCM4 predicts further warming during 2020-2050 of 0.059 ºC/year in July, 0.088 ºC/year in August, and 0.019 ºC/year in October, which will further reduce snow cover.

We study epidemic spreading in complex networks by a multiple random walker approach. Each walker performs an independent simple Markovian random walk on a complex undirected (ergodic) random graph where we focus on Barabási-Albert (BA), Erdös-Rényi (ER) and Watts-Strogatz (WS) types. Both, walkers and nodes can be either susceptible (S) or infected and infectious (I) representing their states of health. Susceptible nodes may be infected by visits of infected walkers, and susceptible walkers may be infected by visiting infected nodes. No direct transmission of the disease among walkers (or among nodes) is possible. This model mimics a large class of diseases such as Dengue and Malaria with transmission of the disease via vectors (mosquitos). Infected walkers may die during the time span of their infection introducing an additional compartment D of dead walkers. Infected nodes never die and always recover from their infection after a random finite time. This assumption is based on the observation that infectious vectors (mosquitos) are not ill and do not die from the infection. The infectious time spans of nodes and walkers, and the survival times of infected walkers, are represented by independent random variables. We derive stochastic evolution equations for the mean-field compartmental populations with mortality of walkers and delayed transitions among the compartments. From linear stability analysis, we derive the basic reproduction numbers R M , R 0 with and without mortality, respectively, and prove that R M < R 0 . For R M , R 0 > 1 the healthy state is unstable whereas for zero mortality a stable endemic equilibrium exists (independent of the initial conditions) which we obtained explicitly. We observe that the solutions of the random walk simulations in the considered networks agree well with the mean-field solutions for strongly connected graph topologies, whereas less well for weakly connected structures and for diseases with high mortality. Our model has applications beyond epidemic dynamics, for instance in the kinetics of chemical reactions, the propagation of contaminants, wood fires, among many others.

This article deals with a finite-source stock-dependent stochastic inventory system with heterogeneous servers and retrial facility. The system can store up to S number of items, and the lifetime of each item follows an exponential distribution. The primary customers arrive in a finite waiting hall from a finite source and get service from the multi-channel (c-heterogeneous servers). The customer’s arrival rate depends on the present stock level. Those customers who arrive while the waiting hall is full enter into the finite orbit. Also, customers in the waiting hall can become impatient and enter the orbit. The (s,Q) ordering policy is adopted for the replenishment of stocks. In the steady state, the joint probability distribution of a number of items in the inventory, busy servers, customers in the waiting hall, and orbit is computed. The comparative numerical analysis of heterogeneous and homogeneous servers on the expected total cost, the average number of customers in the waiting hall and orbit, the fraction of successful retrial rate, and the average impatient customer rate are also carried out.

Continuous greenhouse gas emissions are causing global warming and impacting the habitats of many animals. Researchers in the field of electric power are making efforts to mitigate this situation. Operating and maintaining the power grid in an economic, low-carbon, and stable is challenging. To address the issue, we propose a grid dispatching technique that combines prediction technology, reinforcement learning, and optimization technology. Prediction technology can forecast future power demand and solar power generation, while reinforcement learning and optimization technology can make charging and discharging decisions for energy storage devices based on current and future grid conditions. In the power system, the aggregation of distributed energy resources increases uncertainty, particularly due to the fluctuating generation of renewable energy. This requires the use of advanced predictive control techniques to ensure long-term economic and decarbonization goals. In this paper, we present a real-time dispatching framework that integrates deep learning-based prediction, reinforcement learning-based decision-making, and stochastic optimization techniques. The framework can rapidly adapt to target uncertainty caused by various factors in real-time data distribution and control processes. The proposed framework achieved global Champion in the NeurIPS Challenge 2022 competition and demonstrated its effectiveness in practical scenarios of intelligent building energy management.

A plant’s health and productivity is influenced by its associated microbes. Although the common microbiome is often thought to be the most influential, significant numbers of rare or uncommon microbes (eg. specialized endosymbionts) may also play an important role in the health and productivity of certain plants in certain environments. To help identify rare/specialized bacteria and fungi in the most important angiosperm plants, we contrasted microbiomes of the shoots, roots and rhizospheres of Arabidopsis, Brachypodium, maize, wheat, sugarcane, rice, tomato, coffee, common bean, cassava, soybean, switchgrass, sunflower, Brachiaria, barley, sorghum, and pea. Plants were grown inside sealed jars on sterile sand or field soil. About 95% and 86% of fungal and bacterial diversity inside plants was uncommon, however judging by read abundance, up to half of the mycobiome consists of uncommon fungal cells, while less than 11% of bacterial endophytes are rare. Uncommon seed transmitted microbiomes consisted mostly of Proteobacteria, Firmicutes, Bacteriodetes, Ascomycetes and Basidiomycetes that most heavily colonized shoots, to a lesser extent roots and least of all rhizospheres. Soil served as a more diverse source of rare microbes than seeds, replacing or excluding the majority of the uncommon seed transmitted microbiome. With the rarest microbes, their colonization pattern could either be the result of stringent biotic filtering by most plants, or uneven/stochastic inoculum distribution in seeds or soil. Several strong plant-microbe associations were observed such as seed transmission to shoots, roots and/or rhizospheres of Sarocladium zeae (maize), Penicillium (pea and Phaseolus), and Curvularia (sugarcane), while robust bacterial colonization from cassava field soil occurred with the cyanobacteria Leptolyngbya into Arabidopsis and Panicum roots, and Streptomyces into cassava roots. Some abundant microbes such as Sakaguchia in rice shoots or Vermispora in Arabidopsis roots appeared in no other samples, suggesting they were infrequent, stochastically deposited propagules from either soil or seed (impossible to know based on the available data). Future experiments with culturing and cross inoculation of these microbes between plants may help us better understand host preferences and their role in plant productivity, perhaps leading to their use in crop microbiome engineering and enhancement of agricultural production.

In this paper, we propose a new mathematical model based on the association between susceptible and recovered individual, where the association between susceptible and recovered individual is disturbed by white noise. This model is based on demographic changes and is used for long term behavior. We study the stability of equilibria of the deterministic model and prove the conditions for the extinction of diseases. Then, we investigate and obtain the critical condition of the stochastic epidemic model for the extinction and the permanence in mean of the disease with the white noise. To verify our results, we present some numerical simulations for real data related to disease.

The vibration feature of weak gear fault is often covered in strong background noise, which makes it necessary to establish weak feature enhancement methods. Among the enhancement methods, stochastic resonance (SR) has the unique advantage of transferring noise energy to weak signals and has a great application prospection in weak signal extraction. But the traditional SR potential model cannot form a richer potential structure and may lead to system instability when the noise is too great. To overcome these shortcomings, the article presents a periodic potential underdamping stochastic resonance (PPUSR) method after investigating the potential function and system signal-to-noise ratio (SNR). In addition, system parameters are further optimized by using ant colony algorithm. Through simulation and gear experiments, the effectiveness of the proposed method was verified. We concluded that compared with the traditional underdamped stochastic resonance (TUSR) method, the PPUSR method had a higher recognition degree and better frequency response capability.

Dynamic hedging has been adopted by many insurance companies to mitigate the financial risks associated with variable annuity guarantees. In order to simulate the performance of dynamic hedging for variable annuity products, insurance companies rely on nested stochastic projections, which is highly computationally intensive and often prohibitive for large variable annuity portfolios. Metamodeling techniques have recently been proposed to address the computational issues. However, it is difficult for researchers to obtain real datasets from insurance companies to test metamodeling techniques and publish the results in academic journals. In this paper, we create synthetic datasets that can be used for the purpose of addressing the computational issues associated with the nested stochastic valuation of large variable annuity portfolios. The runtime used to create these synthetic datasets would be about 3 years if a single CPU were used. These datasets are readily available to researchers and practitioners so that they can focus on testing metamodeling techniques.

An approach to estimating the objective function value of minimization maximum lateness problem is proposed. It is shown how to use transformed instances to define a new continuous objective function. After that, using this new objective function, the approach itself is formulated. We calculate the objective function value for some polynomially solvable transformed instances and use them as interpolation nodes to estimate the objective function of the initial instance. What is more, two new polynomial cases, that are easy to use in the approach, are proposed. In the end of the paper numeric experiments are described and their results are provided.

In this work, optimal consumption and investment strategies are derived for risk-averse investors under the 4/2 stochastic volatility class of models. We work under an expected utility (EUT) framework and consider a Constant Relative Risk Aversion (CRRA) investor, who might also be ambiguity-averse. The corresponding Hamilton-Jacobi-Bellman (HJB) and HJB-Isaacs (HJBI) equations are solved in closed-form for a subset of the parametric space and under some restrictions on the portfolio setting, for complete markets. Conditions for proper changes of measure and well defined solutions are provided. These are the first analytical solutions for the 4/2 stochastic volatility model and the embedded 3/2 model for the type of excess returns established in .

Paratemporal methods based on tree expansion have proven to be effective in efficiently generating trajectories of stochastic systems. However, combinatorial explosion of branching arising from multiple choice points presents a major hurdle that must be overcome to implement such techniques. In this paper we tackle this scalability problem by developing a systems theory-based framework covering both conventional and proposed tree expansion algorithms for speeding up discrete event system stochastic simulations while preserving desired accuracy. An example is discussed to illustrate the tree expansion framework in which a discrete event system specification (DEVS) Markov stochastic model takes the form of a tree isomorphic to a free monoid over the branching alphabet. We derive the computation times for baseline, non-merging, and merging tree expansion algorithms to compute the distribution of output values at any given depth. The results show the remarkable reduction from exponential to polynomial dependence on depth effectuated by node merging. We relate these results to the similarly reduced computation time of binomial coefficients underlying Pascal’s triangle. Finally, we discuss application of tree expansion to estimating temporal distributions in stochastic simulations involving serial and parallel compositions with potential real world use cases.

In a particle swarm with a fixed kinetic energy, the constituent particles move randomly and their velocity magnitudes follow a Maxwell distribution characterized by specific parameters. Within a certain time frame, particles in a sub-swarm may exhibit a bias in their movement direction. Thus, the particles may exhibit a special characteristic (biased random motion), where the group velocity in a particular direction remains constant u. For this biased particle swarm, composed of randomly moving particles with a fixed average speed c, particles have a higher probability of moving in a specific direction. Conversely, their likelihood of moving in all other directions is equally reduced. This study starts from the perspective of biased random walks and using Mathematica software proves that the diffusion rate of particles (in the reference frame R_u observed from R₀) in all directions is slower than that of an unbiased particle swarm with the same average speed c. Specifically, the degree of reduction is determined by the Lorentz-like factor. Lastly, we present the Ito equation for this biased random motion and provide associated verification examples. This study aims to serve as a reference for comprehending the underlying principles of the special relativity effect.

The purpose of this paper is to find the time optimal control to a target set for semilinear stochastic functional differential equations involving time delays or memories under general conditions on a target set and nonlinear terms even though the equations contain unbounded principal operators. Our research approach is construct a fundamental solution for corresponding linear systems and establish variations of constant formula of solutions for given stochastic equations. The existence result of time optimal controls for one point target set governed by the given semilinear stochastic equation is also established.

The research develops the vehicular routing problem for the distribution of refrigerated products in a multinational company. The mathematical model corresponds to the vehicle routing problem with hard time windows and stochastic service (VRPTW-ST) time model applied in the city of Santiago de Chile. For optimization of the model, optimization methods tabu search, chaotic search and general algebraic modeling were used. The results allow to validate the meta-heuristic chaotic search to solve the VRPTW-ST; chaotic search method obtains superior performance than tabu search method for solving a real problem in a large city.

Based on a previous risk calculation study along a road corridor, risk is recalculated using stochastic simulation by introducing variability for most of the parameters in the risk equation. This leads to an exceedance curve comparable to that of catastrophe models. This approach introduces uncertainty into the risk calculation in a simple way, which can be used for poorly documented cases to fulfil lack of data. This approach seems to tend to minimize risk or to question risk calculations.

Climate justice as a commons is conceived as the intertemporal climate equity and equal- ity exchange amongst generations. Sustainability, intended as the interplay amongst the economy, the society, the environment, and the governance, is essential to forge the climate justice theoretical framework. On this base, the study attempts to model intertemporal choice amongst generations in these four domains, making use of an over- lapping generations (OLG) model. The proxies detected are GDP growth (economy), environmental quality (environment), and labor growth, and environmental investment (society) as assumptions. The governance dimension is captured by the di¤erence in wealth between young and old generations. The work aims at replying to the follow- ing research question: Which are the conditions for sustainable development such that climate justice holds? The intra-intergenerational exchange is de…ned in two periods, while the individual provides their preferred economic and environmental choice mix as consumption-saving. This study shows that sustainable growth is achievable only with increased young e¤ort and less leisure and consumption.

The purpose of this study is to estimate the mean transitioning probabilities from a Healthy state to malaria positive uncomplicated state or to malaria positive severe state. It also classifies the various transitioning probabilities of moving through the various states based on some baseline characteristics. Malaria test results for 2019 over a 12-month period were collected from the University of Ghana school clinic. An H-U model for the study was developed and the transition rates from the cross-sectional data are indicated. With two states Healthy (H) and Uncomplicated (U) forming a state space, there were four possible transitions. The results show that the probability of transitioning from a Healthy state to a malaria positive state is 0.03% while the probability that an individual will remain at Healthy state (H) after the test is 99.73%. It was found that if an individual is already positive and has taken medication the probability that its second test came out negative is 6.45% while the chances that it will remain positive but uncomplicated is 93.55%. The study also showed that in the long run, about 95.98% of persons who visited the student clinic with malaria symptoms recorded negative tests for malaria parasite while about 4% recorded positive for malaria. In terms of disaggregation by gender, it was realized that the number of reported negative test results were higher for females (97.08%) than for males (96.13%). However, the infection rate is higher for males (3.87%) than females (2.92%). It is recommended that in as much as the University of Ghana has two health centers (a clinic and hospital), there should be a centralized system to track students’ health so research done would not be biased.

In this paper, we obtain the existence and uniqueness theorem for solutions of Caputo type fractional stochastic delay differential systems (FSDDSs) with Poisson jumps by utilizing delayed perturbation of Mittag-Leffler function. Moreover, by using Burkholder-Davis-Gundy's inequality, Doob's martingale inequality and Holder inequality, we prove that the solution of the averaged FSDDSs converges to that of the standard FSDDSs in the sense of Lp. Some known results in the literature are extended.

The paper describes a method of solving some stochastic processes using generating functions. A general theorem of generating functions of a particular type is derived. A generating function of this type is applied to a stochastic process yielding polynomial time algorithms for certain partitions. The method is generalized to a stochastic process describing a rather general linear transform. Finally, the main idea of the method is used in deriving a theoretical polynomial time algorithm to the knapsack problem.

Research within the field of hydrology often focuses on comparing stochastic to machine learning (ML) forecasting methods. The comparisons performed are all based on case studies, while an extensive study aiming to provide generalized results on the subject is missing. Herein, we compare 11 stochastic and 9 ML methods regarding their multi-step ahead forecasting properties by conducting 12 large-scale computational experiments based on simulations. Each of these experiments uses 2 000 time series generated by linear stationary stochastic processes. We conduct each simulation experiment twice; the first time using time series of 100 values and the second time using time series of 300 values. Additionally, we conduct a real-world experiment using 405 mean annual river discharge time series of 100 values. We quantify the performance of the methods using 18 metrics. The results indicate that stochastic and ML methods perform equally well.

Using multi-color flow cytometry analysis, we studied the immunophenotypical differences between leukemic blasts/LSCs from patients with AML/MDS and hematopoietic stem and progenitor cells (HSPCs) from patients in complete remission (CR) following their successful treatment. The panel of markers included CD34, CD38, CD45RA, CD123 as representatives for a hierarchical hematopoietic stem and progenitor cell (HSPC) classification as well as programmed death ligand 1 (PD-L1). Rather than restricting the evaluation on a 2- or 3-dimensional analysis, we applied a t-distributed stochastic neighbor embedding (t-SNE) approach to get a deeper insight and segregation between leukemic and normal HPSCs. For that purpose, we created a t-SNE map, which resulted in the visualization of 27 cell clusters based on their similarity concerning the composition and intensity of antigen expression. Two of these clusters were “leukemia-related” containing a great proportion of CD34⁺/CD38^- hematopoietic stem cells (HSCs) or CD34⁺ cells with a strong coexpression of CD45RA/CD123, respectively. CD34⁺ cells within the latter cluster were also highly positive for PD-L1 reflecting their immunosuppressive capacity. Beyond this proof of principle study, the inclusion of additional markers will be helpful to refine the differentiation between normal HSPCs and leukemic cells, particularly in the context of minimal disease detection and antigen-targeted therapeutic interventions. Furthermore, we suggest a protocol for the assignment of new cell ensembles in quantitative terms, via a numerical value, the Pearson coefficient, based on a similarity comparison of the t-SNE pattern with a reference.