Assessment of Wet Season Precipitation in the Central United States by the Regional Climate Simulation of the WRFG Member in NARCCAP and Its Relationship with Large-Scale Circulation Biases

Assessment of past-climate simulations of regional climate models(RCMs)is important for understanding the reliability of RCMs when used to project future regional climate.Here,we assess the performance and discuss possible causes of biases in a WRF-based RCM with a grid spacing of 50 km,named WRFG,from the North American Regional Climate Change Assessment Program(NARCCAP)in simulating wet season precipitation over the Central United States for a period when observational data are available.The RCM reproduces key features of the precipitation distribution characteristics during late spring to early summer,although it tends to underestimate the magnitude of precipitation.This dry bias is partially due to the model’s lack of skill in simulating nocturnal precipitation related to the lack of eastward propagating convective systems in the simulation.Inaccuracy in reproducing large-scale circulation and environmental conditions is another contributing factor.The too weak simulated pressure gradient between the Rocky Mountains and the Gulf of Mexico results in weaker southerly winds in between,leading to a reduction of warm moist air transport from the Gulf to the Central Great Plains.The simulated low-level horizontal convergence fields are less favorable for upward motion than in the NARR and hence,for the development of moist convection as well.Therefore,a careful examination of an RCM’s deficiencies and the identification of the source of errors are important when using the RCM to project precipitation changes in future climate scenarios.

Interaction between systemic iron parameters and left ventricular structure and function in the preserved ejection fraction population:a two-sample bidirectional Mendelian randomization study

BACKGROUND Left ventricular(LV)remodeling and diastolic function in people with heart failure(HF)are correlated with iron status;however,the causality is uncertain.This Mendelian randomization(MR)study investigated the bidirectional causal relationship between systemic iron parameters and LV structure and function in a preserved ejection fraction population.METHODS Transferrin saturation(TSAT),total iron binding capacity(TIBC),and serum iron and ferritin levels were extracted as instrumental variables for iron parameters from meta-analyses of public genome-wide association studies.Individuals without myocardial infarction history,HF,or LV ejection fraction(LVEF)<50%(n=16,923)in the UK Biobank Cardiovascular Magnetic Resonance Imaging Study constituted the outcome dataset.The dataset included LV end-diastolic volume,LV endsystolic volume,LV mass(LVM),and LVM-to-end-diastolic volume ratio(LVMVR).We used a two-sample bidirectional MR study with inverse variance weighting(IVW)as the primary analysis method and estimation methods using different algorithms to improve the robustness of the results.RESULTS In the IVW analysis,one standard deviation(SD)increased in TSAT significantly correlated with decreased LVMVR(β=-0.1365;95%confidence interval[CI]:-0.2092 to-0.0638;P=0.0002)after Bonferroni adjustment.Conversely,no significant relationships were observed between other iron and LV parameters.After Bonferroni correction,reverse MR analysis showed that one SD increase in LVEF significantly correlated with decreased TSAT(β=-0.0699;95%CI:-0.1087 to-0.0311;P=0.0004).No heterogeneity or pleiotropic effects evidence was observed in the analysis.CONCLUSIONS We demonstrated a causal relationship between TSAT and LV remodeling and function in a preserved ejection fraction population.

A Two-Layer Encoding Learning Swarm Optimizer Based on Frequent Itemsets for Sparse Large-Scale Multi-Objective Optimization

Traditional large-scale multi-objective optimization algorithms(LSMOEAs)encounter difficulties when dealing with sparse large-scale multi-objective optimization problems(SLM-OPs)where most decision variables are zero.As a result,many algorithms use a two-layer encoding approach to optimize binary variable Mask and real variable Dec separately.Nevertheless,existing optimizers often focus on locating non-zero variable posi-tions to optimize the binary variables Mask.However,approxi-mating the sparse distribution of real Pareto optimal solutions does not necessarily mean that the objective function is optimized.In data mining,it is common to mine frequent itemsets appear-ing together in a dataset to reveal the correlation between data.Inspired by this,we propose a novel two-layer encoding learning swarm optimizer based on frequent itemsets(TELSO)to address these SLMOPs.TELSO mined the frequent terms of multiple particles with better target values to find mask combinations that can obtain better objective values for fast convergence.Experi-mental results on five real-world problems and eight benchmark sets demonstrate that TELSO outperforms existing state-of-the-art sparse large-scale multi-objective evolutionary algorithms(SLMOEAs)in terms of performance and convergence speed.

Enhancing Evolutionary Algorithms With Pattern Mining for Sparse Large-Scale Multi-Objective Optimization Problems

Sparse large-scale multi-objective optimization problems(SLMOPs)are common in science and engineering.However,the large-scale problem represents the high dimensionality of the decision space,requiring algorithms to traverse vast expanse with limited computational resources.Furthermore,in the context of sparse,most variables in Pareto optimal solutions are zero,making it difficult for algorithms to identify non-zero variables efficiently.This paper is dedicated to addressing the challenges posed by SLMOPs.To start,we introduce innovative objective functions customized to mine maximum and minimum candidate sets.This substantial enhancement dramatically improves the efficacy of frequent pattern mining.In this way,selecting candidate sets is no longer based on the quantity of nonzero variables they contain but on a higher proportion of nonzero variables within specific dimensions.Additionally,we unveil a novel approach to association rule mining,which delves into the intricate relationships between non-zero variables.This novel methodology aids in identifying sparse distributions that can potentially expedite reductions in the objective function value.We extensively tested our algorithm across eight benchmark problems and four real-world SLMOPs.The results demonstrate that our approach achieves competitive solutions across various challenges.

A semantic vector map-based approach for aircraft positioning in GNSS/GPS denied large-scale environment

Accurate positioning is one of the essential requirements for numerous applications of remote sensing data,especially in the event of a noisy or unreliable satellite signal.Toward this end,we present a novel framework for aircraft geo-localization in a large range that only requires a downward-facing monocular camera,an altimeter,a compass,and an open-source Vector Map(VMAP).The algorithm combines the matching and particle filter methods.Shape vector and correlation between two building contour vectors are defined,and a coarse-to-fine building vector matching(CFBVM)method is proposed in the matching stage,for which the original matching results are described by the Gaussian mixture model(GMM).Subsequently,an improved resampling strategy is designed to reduce computing expenses with a huge number of initial particles,and a credibility indicator is designed to avoid location mistakes in the particle filter stage.An experimental evaluation of the approach based on flight data is provided.On a flight at a height of 0.2 km over a flight distance of 2 km,the aircraft is geo-localized in a reference map of 11,025 km~2using 0.09 km~2aerial images without any prior information.The absolute localization error is less than 10 m.

Online identification and extraction method of regional large-scale adjustable load-aggregation characteristics

This article introduces the concept of load aggregation,which involves a comprehensive analysis of loads to acquire their external characteristics for the purpose of modeling and analyzing power systems.The online identification method is a computer-involved approach for data collection,processing,and system identification,commonly used for adaptive control and prediction.This paper proposes a method for dynamically aggregating large-scale adjustable loads to support high proportions of new energy integration,aiming to study the aggregation characteristics of regional large-scale adjustable loads using online identification techniques and feature extraction methods.The experiment selected 300 central air conditioners as the research subject and analyzed their regulation characteristics,economic efficiency,and comfort.The experimental results show that as the adjustment time of the air conditioner increases from 5 minutes to 35 minutes,the stable adjustment quantity during the adjustment period decreases from 28.46 to 3.57,indicating that air conditioning loads can be controlled over a long period and have better adjustment effects in the short term.Overall,the experimental results of this paper demonstrate that analyzing the aggregation characteristics of regional large-scale adjustable loads using online identification techniques and feature extraction algorithms is effective.

Chronic Hepatitis B Virus Infection: Biological Parameters in Patients Treated with Tenofovir Disoproxil Fumarate

Chronic hepatitis B causes a liver disease characterized by inflammation of the liver parenchyma. The aim of this study was to investigate the evolution of biological parameters in patients treated with Tenofovir for chronic B infection at the Commune V referral health center in Bamako. We obtained a prevalence of 14.15%. The most represented age group was 31 - 40 years, with 36.8%. The sex ratio was 1.44 in favour of men. Viral load was undetectable after 18 months of treatment in 25 patients (42.37%). Tenofovir, the 1st-line drug in Mali, is effective on the biological parameters monitored in patients.

Large-Scale Multi-Objective Optimization Algorithm Based on Weighted Overlapping Grouping of Decision Variables

The large-scale multi-objective optimization algorithm(LSMOA),based on the grouping of decision variables,is an advanced method for handling high-dimensional decision variables.However,in practical problems,the interaction among decision variables is intricate,leading to large group sizes and suboptimal optimization effects;hence a large-scale multi-objective optimization algorithm based on weighted overlapping grouping of decision variables(MOEAWOD)is proposed in this paper.Initially,the decision variables are perturbed and categorized into convergence and diversity variables;subsequently,the convergence variables are subdivided into groups based on the interactions among different decision variables.If the size of a group surpasses the set threshold,that group undergoes a process of weighting and overlapping grouping.Specifically,the interaction strength is evaluated based on the interaction frequency and number of objectives among various decision variables.The decision variable with the highest interaction in the group is identified and disregarded,and the remaining variables are then reclassified into subgroups.Finally,the decision variable with the strongest interaction is added to each subgroup.MOEAWOD minimizes the interactivity between different groups and maximizes the interactivity of decision variables within groups,which contributed to the optimized direction of convergence and diversity exploration with different groups.MOEAWOD was subjected to testing on 18 benchmark large-scale optimization problems,and the experimental results demonstrate the effectiveness of our methods.Compared with the other algorithms,our method is still at an advantage.

A Large-Scale Group Decision Making Model Based on Trust Relationship and Social Network Updating

With the development of big data and social computing,large-scale group decisionmaking(LGDM)is nowmerging with social networks.Using social network analysis(SNA),this study proposes an LGDM consensus model that considers the trust relationship among decisionmakers(DMs).In the process of consensusmeasurement:the social network is constructed according to the social relationship among DMs,and the Louvain method is introduced to classify social networks to form subgroups.In this study,the weights of each decision maker and each subgroup are computed by comprehensive network weights and trust weights.In the process of consensus improvement:A feedback mechanism with four identification and two direction rules is designed to guide the consensus of the improvement process.Based on the trust relationship among DMs,the preferences are modified,and the corresponding social network is updated to accelerate the consensus.Compared with the previous research,the proposedmodel not only allows the subgroups to be reconstructed and updated during the adjustment process,but also improves the accuracy of the adjustment by the feedbackmechanism.Finally,an example analysis is conducted to verify the effectiveness and flexibility of the proposed method.Moreover,compared with previous studies,the superiority of the proposed method in solving the LGDM problem is highlighted.

Study on the clogging mechanism of punching screen in sand control by the punching structure parameters

As an independent sand control unit or a common protective shell of a high-quality screen,the punching screen is the outermost sand retaining unit of the sand control pipe which is used in geothermal well or oil and gas well.However,most screens only consider the influence of the internal sand retaining medium parameters in the sand control performance design while ignoring the influence of the plugging of the punching screen on the overall sand retaining performance of the screen.To explore the clogging mechanism of the punching screen,this paper established the clogging mechanism calculation model of a single punching screen sand control unit by using the computational fluid mechanics-discrete element method(CFD-DEM)combined method.According to the combined motion of particles and fluids,the influence of the internal flow state on particle motion and accumulation was analyzed.The results showed that(1)the clogging process of the punching sand control unit is divided into three stages:initial clogging,aggravation of clogging and stability of clogging.In the initial stage of blockage,coarse particles form a loose bridge structure,and blockage often occurs preferentially at the streamline gathering place below chamfering inside the sand control unit.In the stage of blockage intensification,the particle mass develops into a relatively complete sand bridge,which develops from both ends of the opening to the center of the opening.In the stable plugging stage,the sand deposits show a“fan shape”and form a“V-shaped”gully inside the punching slot element.(2)Under a certain reservoir particle-size distribution,The slit length and opening height have a large influence on the permeability and blockage rate,while the slit width size has little influence on the permeability and blockage rate.The microscopic clogging mechanism and its law of the punching screen prevention unit are proposed in this study,which has some field guidance significance for the design of punching screen and sand prevention selection.

Quaternion-Based Adaptive Trajectory Tracking Control of a Rotor-Missile with Unknown Parameters Identification

This paper investigates the adaptive trajectory tracking control problem and the unknown parameter identification problem of a class of rotor-missiles with parametric system uncertainties.First,considering the uncertainty of structural and aerodynamic parameters,the six-degree-of-freedom(6Do F) nonlinear equations describing the position and attitude dynamics of the rotor-missile are established,respectively,in the inertial and body-fixed reference frames.Next,a hierarchical adaptive trajectory tracking controller that can guarantee closed-loop stability is proposed according to the cascade characteristics of the 6Do F dynamics.Then,a memory-augmented update rule of unknown parameters is proposed by integrating all historical data of the regression matrix.As long as the finitely excited condition is satisfied,the precise identification of unknown parameters can be achieved.Finally,the validity of the proposed trajectory tracking controller and the parameter identification method is proved through Lyapunov stability theory and numerical simulations.

Investigating the Evolution of Amati Parameters with Redshift

Gamma-ray bursts(GRBs) are among the brightest objects in the Universe and, hence, can be observed up to a very high redshift. Properly calibrated empirical correlations between intensity and spectral correlations of GRBs can be used to estimate the cosmological parameters. However, the possibility of the evolution of GRBs with redshift is a long-standing puzzle. In this work, we used 162 long-duration GRBs to determine whether GRBs below and above a certain redshift have different properties. The GRBs are split into two groups, and we fit the Amati relation for each group separately. Our findings demonstrate that estimations of the Amati parameters for the two groups are substantially dissimilar. We perform simulations to investigate whether the selection effects could cause the difference. Our analysis shows that the differences may be either intrinsic or due to systematic errors in the data, and the selection effects are not their true origin. However, in-depth analysis with a new data set comprised of 119 long GRBs shows that intrinsic scatter may partly be responsible for such effects.

Comparison of ocular biometric parameters between two swept-source optical coherence tomography devices and Scheimpflug tomography in patients with cataract

AIM:To assess and compare the variations and agreements across different ocular biometric parameters using swept-source optical coherence tomography(SS-OCT)and Scheimpflug tomography in patients diagnosed with cataract.METHODS:This prospective case series was conducted at Tianjin Medical University Eye Hospital.In total,212 eyes from 212 patients scheduled for phacoemulsification were included.Eyes were evaluated preoperatively using two SSOCT devices(IOLMaster700 and CASIA2)and Scheimpflug tomography(Pentacam).Central corneal thickness(CCT),anterior chamber depth(ACD),aqueous depth(AQD),white-to-white distance(WTW),flat simulated keratometry(Kf),steep simulated keratometry(Ks),mean keratometry(Km),and total corneal keratometry(TKm)were measured.Intraclass correlation coefficient(ICC),95%confidence intervals(CI)and limits of agreement(LoA)widths were conducted to assess differences and correlations between devices.RESULTS:All parameters,except for Ks,were significantly different.Pairwise comparison revealed no significant differences between keratometry obtained by IOLMaster 700 and Pentacam.LoA widths of all paired comparisons for Ks were>0.80 D.Except for WTW between IOLMaster 700 and CASIA2 and between CASIA2 and Pentacam,other Pearson’s coefficients between devices showed a strong correlation(all r>0.95).The ICC of WTW(ICC=0.438,95%CI 0.167-0.625)showed poor reliability.The reliability of CCT,ACD,and AQD was excellent(all ICC>0.95),whereas that of TKm was good(ICC=0.827,95%CI 0.221-0.939).A significant linear correlation was also observed among devices.CONCLUSION:The ocular parameters derived from the use of IOLMaster700,CASIA2,and Pentacam exhibit significant discrepancies;as such,measurements from these devices should not be deemed as interchangeable.

Digital measuring the ocular morphological parameters of guinea pig eye in vivo with Python

AIM:To quantitatively measure ocular morphological parameters of guinea pig with Python technology.METHODS:Thirty-six eyeballs of eighteen 3-weekold guinea pigs were measured with keratometer and photographed to obtain the horizontal,coronal,and sagittal planes respectively.The corresponding photo pixels-actual length ratio was acquired by a proportional scale.The edge coordinates were identified artificially by ginput function.Circle and conic curve fitting were applied to fit the contour of the eyeball in the sagittal,coronal and horizontal view.The curvature,curvature radius,eccentricity,tilt angle,corneal diameter,and binocular separation angle were calculated according to the geometric principles.Next,the eyeballs were removed,canny edge detection was applied to identify the contour of eyeball in vitro.The results were compared between in vivo and in vitro.RESULTS:Regarding the corneal curvature and curvature radius on the horizontal and sagittal planes,no significant differences were observed among results in vivo,in vitro,and the keratometer.The horizontal and vertical binocular separation angles were 130.6°±6.39°and 129.8°±9.58°respectively.For the corneal curvature radius and eccentricity in vivo,significant differences were observed between horizontal and vertical planes.CONCLUSION:The Graphical interface window of Python makes up the deficiency of edge detection,which requires too much definition in Matlab.There are significant differences between guinea pig and human beings,such as exotropic eye position,oblique oval eyeball,and obvious discrepancy of binoculus.This study helps evaluate objectively the ocular morphological parameters of small experimental animals in emmetropization research.

Analysis and comparison of retinal vascular parameters under different glucose metabolic status based on deep learning

AIM:To develop a deep learning-based model for automatic retinal vascular segmentation,analyzing and comparing parameters under diverse glucose metabolic status(normal,prediabetes,diabetes)and to assess the potential of artificial intelligence(AI)in image segmentation and retinal vascular parameters for predicting prediabetes and diabetes.METHODS:Retinal fundus photos from 200 normal individuals,200 prediabetic patients,and 200 diabetic patients(600 eyes in total)were used.The U-Net network served as the foundational architecture for retinal arteryvein segmentation.An automatic segmentation and evaluation system for retinal vascular parameters was trained,encompassing 26 parameters.RESULTS:Significant differences were found in retinal vascular parameters across normal,prediabetes,and diabetes groups,including artery diameter(P=0.008),fractal dimension(P=0.000),vein curvature(P=0.003),C-zone artery branching vessel count(P=0.049),C-zone vein branching vessel count(P=0.041),artery branching angle(P=0.005),vein branching angle(P=0.001),artery angle asymmetry degree(P=0.003),vessel length density(P=0.000),and vessel area density(P=0.000),totaling 10 parameters.CONCLUSION:The deep learning-based model facilitates retinal vascular parameter identification and quantification,revealing significant differences.These parameters exhibit potential as biomarkers for prediabetes and diabetes.

Effects of Fresh Cupressus sempervirens Leaves Infusion on Growth Performance, Intestinal Microbiota and Haemato-Biochemical Parameters in Broilers

The accumulation of growth-promoting antibiotic residues in animal products and the resistance developed by bacteria in poultry farms has led to a search for natural compounds derived from plants. This study was designed to promote the production performance of broiler chickens using fresh Cupressus sempervirens leaves infusion. Fresh Cupressus sempervirens leaves were harvested, washed, chopped and ground to a paste using a blender and fermented for three days in a closed container at a rate of 500 g/L of water. The solution obtained was filtered and added at the rate of 2, 4, 6, 8 and 10 ml/L of drinking water. The chickens fed on the graded level of the solution were compared to a control ration without an additive and positive control group supplemented with 1 g antibiotic/kg feed. At the finisher phase and throughout the study period, water intake increased significantly (P < 0.05) with increasing levels of infusion. Feed intake decreased significantly (P < 0.05) with 2 and 4 ml of infusion/L drinking water. Live weight and weight gain were significantly (P < 0.05) higher with 6 ml/L, while feed conversion significantly (P < 0.05) decreased with the same treatment compared with the control treatment without additives (T0). Carcass characteristics were not significantly (P > 0.05) affected by the inclusion of Cupressus sempervirens infusion. Haematological parameters significantly (P < 0.05) increase independently of the rate of incorporation of the infusion into the drinking water, with the exception of RBCs, MCHT and PCT. Serum content in total protein, globulins, LDL cholesterol and triglycerides were significantly (P < 0.05) high with 8 and 10 ml Cupressus sempervirens infusion/litre drinking water as compared to all other treatments. AST, ALT, urea, creatine, albumin, total cholesterol and HDL-cholesterol were not significantly affected. The lactic acid bacteria load increased significantly (P E. coli and salmonella counts decreased significantly (P < 0.05) with infusion compared to the control without additive. In conclusion, 6 ml of Cupressus sempervirens infusion can be used as an alternative to antibiotic feed additives to promote growth performance in broilers.

Determination of uncertainties of geomechanical parameters of metamorphic rocks using petrographic analyses

Geomechanical parameters of intact metamorphic rocks determined from laboratory testing remain highly uncertain because of the great intrinsic variability associated with the degrees of metamorphism.The aim of this paper is to develop a proper methodology to analyze the uncertainties of geomechanical characteristics by focusing on three domains,i.e.data treatment process,schistosity angle,and mineralogy.First,the variabilities of the geomechanical laboratory data of Westwood Mine(Quebec,Canada)were examined statistically by applying different data treatment techniques,through which the most suitable outlier methods were selected for each parameter using multiple decision-making criteria and engineering judgment.Results indicated that some methods exhibited better performance in identifying the possible outliers,although several others were unsuccessful because of their limitation in large sample size.The well-known boxplot method might not be the best outlier method for most geomechanical parameters because its calculated confidence range was not acceptable according to engineering judgment.However,several approaches,including adjusted boxplot,2MADe,and 2SD,worked very well in the detection of true outliers.Also,the statistical tests indicate that the best-fitting probability distribution function for geomechanical intact parameters might not be the normal distribution,unlike what is assumed in most geomechanical studies.Moreover,the negative effects of schistosity angle on the uniaxial compressive strength(UCS)variabilities were reduced by excluding the samples within a specific angle range where the UCS data present the highest variation.Finally,a petrographic analysis was conducted to assess the associated uncertainties such that a logical link was found between the dispersion and the variabilities of hard and soft minerals.

A Novel On-Site-Real-Time Method for Identifying Characteristic Parameters Using Ultrasonic Echo Groups and Neural Network

On-site and real-time non-destructive measurement of elastic constants for materials of a component in a in-service structure is a challenge due to structural complexities,such as ambiguous boundary,variable thickness,nonuniform material properties.This work develops for the first time a method that uses ultrasound echo groups and artificial neural network(ANN)for reliable on-site real-time identification of material parameters.The use of echo groups allows the use of lower frequencies,and hence more accommodative to structural complexity.To train the ANNs,a numerical model is established that is capable of computing the waveform of ultrasonic echo groups for any given set of material properties of a given structure.The waveform of an ultrasonic echo groups at an interest location on the surface the structure with material parameters varying in a predefined range are then computed using the numerical model.This results in a set of dataset for training the ANN model.Once the ANN is trained,the material parameters can be identified simultaneously using the actual measured echo waveform as input to the ANN.Intensive tests have been conducted both numerically and experimentally to evaluate the effectiveness and accuracy of the currently proposed method.The results show that the maximum identification error of numerical example is less than 2%,and the maximum identification error of experimental test is less than 7%.Compared with currently prevailing methods and equipment,the proposefy the density and thickness,in addition to the elastic constants.Moreover,the reliability and accuracy of inverse prediction is significantly improved.Thus,it has broad applications and enables real-time field measurements,which has not been fulfilled by any other available methods or equipment.

Surface Water Quality Profiling Using Physicochemical Parameters in Open Defecation Free and Non-Open Defecation Free Local Government Areas in Benue State, Nigeria

Physicochemical parameters of surface water sources in the study of local government areas (LGAs) were assessed using standard procedures. The mean physicochemical parameters for pH (5.49), NO2 (0.23 mg/L), SO4 (0.77 mg/L), Na (28.72 mg/L), Ca (28.94 mg/L), Mg (17.50 mg/L), Cl (11.65 mg/L), TSS (6.27 mg/L), TDS (104.23 mg/L), BOD (2.83 mg/L) and F (0.87 mg/L) were below WHO standards irrespective of their defecation status. The values for electrical conductivity (EC) (2770.50 µs/cm, turbidity (481.24 NTU), dissolved oxygen (DO) (5.32 mg/L), chemical oxygen demand (COD) (445.50 mg/L), K (125.06 mg/L), PO4 (0.78 mg/L) and Fe (0.57 mg/L) were above the WHO limits for safe water. Higher EC and COD values obtained in the study is evidence of pollution of the water sources by organic matter.

Factorized Smith Method for A Class of High-Ranked Large-ScaleТ-Stein Equations

We introduce a factorized Smith method(FSM)for solving large-scale highranked J-Stein equations within the banded-plus-low-rank structure framework.To effectively reduce both computational complexity and storage requirements,we develop techniques including deflation and shift,partial truncation and compression,as well as redesign the residual computation and termination condition.Numerical examples demonstrate that the FSM outperforms the Smith method implemented with a hierarchical HODLR structured toolkit in terms of CPU time.