Parallel Monte Carlo sampling scheme for sphere and hemisphere

The sampling of certain solid angle is a fundamental operation in realistic image synthesis, where the rendering equation describing the light propagation in closed domains is solved. Monte Carlo methods for solving the rendering equation use sampling of the solid angle subtended by unit hemisphere or unit sphere in order to perform the numerical integration of the rendering equation. In this work we consider the problem for generation of uniformly distributed random samples over hemisphere and sphere. Our aim is to construct and study the parallel sampling scheme for hemisphere and sphere. First we apply the symmetry property for partitioning of hemisphere and sphere. The domain of solid angle subtended by a hemisphere is divided into a number of equal sub-domains. Each sub-domain represents solid angle subtended by orthogonal spherical triangle with fixed vertices and computable parameters. Then we introduce two new algorithms for sampling of orthogonal spherical triangles. Both algorithms are based on a transformation of the unit square. Similarly to the Arvo's algorithm for sampling of arbitrary spherical triangle the suggested algorithms accommodate the stratified sampling. We derive the necessary transformations for the algorithms. The first sampling algorithm generates a sample by mapping of the unit square onto orthogonal spherical triangle. The second algorithm directly compute the unit radius vector of a sampling point inside to the orthogonal spherical triangle. The sampling of total hemisphere and sphere is performed in parallel for all sub-domains simultaneously by using the symmetry property of partitioning. The applicability of the corresponding parallel sampling scheme for Monte Carlo and Quasi-D/lonte Carlo solving of rendering equation is discussed.

A new solution of the rendering equation with stratified Monte Carlo approach

This paper is turned to the advanced Monte Carlo methods for realistic image creation. It offers a new stratified approach for solving the rendering equation. We consider the numerical solution of the rendering equation by separation of integration domain. The hemispherical integration domain is symmetrically separated into 16 parts. First 9 sub-domains are equal size of orthogonal spherical triangles. They are symmetric each to other and grouped with a common vertex around the normal vector to the surface. The hemispherical integration domain is completed with more 8 sub-domains of equal size spherical quadrangles, also symmetric each to other. All sub-domains have fixed vertices and computable parameters. The bijections of unit square into an orthogonal spherical triangle and into a spherical quadrangle are derived and used to generate sampling points. Then, the symmetric sampling scheme is applied to generate the sampling points distributed over the hemispherical integration domain. The necessary transformations are made and the stratified Monte Carlo estimator is presented. The rate of convergence is obtained and one can see that the algorithm is of super-convergent type.

Quasi Monte Carlo approach with uniform separation for solving of the rendering equation, ICNAAM 2007

This paper is directed to the advanced parallel Quasi Monte Carlo (QMC) methods for realistic image synthesis. We propose and consider a new QMC approach for solving the rendering equation with uniform separation. First, we apply the symmetry property for uniform separation of the hemispherical integration domain into 24 equal sub-domains of solid angles, subtended by orthogonal spherical triangles with fixed vertices and computable parameters. Uniform separation allows to apply parallel sampling scheme for numerical integration. Finally, we apply the stratified QMC integration method for solving the rendering equation. The superiority our QMC approach is proved.

Global asymptotic stability in a class of Putnam-type equations

In this paper, we initiate the study of a class of Putnam-type equation of the form x(n-1) = A(1)x(n) + A(2)x(n-1) + A(3)x(n-2)x(n-3) + A(4)/B(1)x(n)x(n-1) + B(2)x(n-2) + B(3)x(n-3) + B-4 n = 0, 1, 2,..., where A(1), A(2), A(3), A(4), B-1, B-2, B-3, B-4 are positive constants with A(1) + A(2) + A(3) + A(4) = B-1 + B-2 + B-3 + B-4, x(-3), x(-2), x(-1), x(0) are positive numbers. A sufficient condition is given for the global asymptotic stability of the equilibrium point c = 1 of such equations. (c) 2005 Elsevier Ltd. All rights reserved.

On two rational difference equations

In this paper, we study the oscillating property of positive solutions and the global asymptotic stability of the unique equilibrium of the two rational difference equations [GRAPHICS] and [GRAPHICS] where a is a nonnegative constant. (c) 2005 Elsevier Inc. All rights reserved.

An efficient flux difference splitting algorithm for unsteady duct flows of a real gas

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a real gas in a single spatial coordinate. This includes flow in a duct of variable cross-section, as well as flow with slab, cylindrical or spherical symmetry, as well as the case of an ideal gas, and can be useful when testing codes for the two-dimensional equations governing compressible flow of a real gas. The resulting scheme requires an average of the flow variables across the interface between cells, and this average is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual “square root” averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and for a number of equations of state. The results compare favourably with the results from other schemes.

A non-uniform mesh scheme for compressible flow

A non-uniform mesh scheme is presented for the computation of compressible flows governed by the Euler equations of gas dynamics. The scheme is based on flux-difference splitting and represents an extension of a similar scheme designed for uniform meshes. The numerical results demonstrate that little, if any, spurious oscillation occurs as a result of the non-uniformity of the mesh; and importantly, shock speeds are computed correctly.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Error estimates for a fully discrete spectral scheme for a class of nonlinear, nonlocal dispersive wave equations

We analyze a fully discrete spectral method for the numerical solution of the initial- and periodic boundary-value problem for two nonlinear, nonlocal, dispersive wave equations, the Benjamin–Ono and the Intermediate Long Wave equations. The equations are discretized in space by the standard Fourier–Galerkin spectral method and in time by the explicit leap-frog scheme. For the resulting fully discrete, conditionally stable scheme we prove an L2-error bound of spectral accuracy in space and of second-order accuracy in time.

Analyzing diffusion tensor images with ghosting artifacts: the effects of direct and indirect normalization

The current study aims to assess the applicability of direct or indirect normalization for the analysis of fractional anisotropy (FA) maps in the context of diffusion-weighted images (DWIs) contaminated by ghosting artifacts. We found that FA maps acquired by direct normalization showed generally higher anisotropy than indirect normalization, and the disparities were aggravated by the presence of ghosting artifacts in DWIs. The voxel-wise statistical comparisons demonstrated that indirect normalization reduced the influence of artifacts and enhanced the sensitivity of detecting anisotropy differences between groups. This suggested that images contaminated with ghosting artifacts can be sensibly analyzed using indirect normalization.

A molecular T-matrix approach to calculating low-energy electron diffusion intensities for ordered molecular adsorbates

We present a novel approach to calculating Low-Energy Electron Diffraction (LEED) intensities for ordered molecular adsorbates. First, the intra-molecular multiple scattering is computed to obtain a non-diagonal molecular T-matrix. This is then used to represent the entire molecule as a single scattering object in a conventional LEED calculation, where the Layer Doubling technique is applied to assemble the different layers, including the molecular ones. A detailed comparison with conventional layer-type LEED calculations is provided to ascertain the accuracy of this scheme of calculation. Advantages of this scheme for problems involving ordered arrays of molecules adsorbed on surfaces are discussed.

Comparison of algal and fish sources on the oxidative stability of poultry meat and its enrichment with omega-3 polyunsaturated fatty acids

Human consumption of long-chain n-3 polyunsaturated fatty acids (LC n-3 PUFA) is below recommendations, and enriching chicken meat (by incorporating LC n-3 PUFA into broiler diets) is a viable means of increasing consumption. Fish oil is the most common LC n-3 PUFA supplement used but is unsustainable and reduces the oxidative stability of the meat. The objective of this experiment was to compare fresh fish oil (FFO) with fish oil encapsulated (EFO) in a gelatin matrix (to maintain its oxidative stability) and algal biomass at a low (LAG, 11), medium (MAG, 22), or high (HAG, 33 g/kg of diet) level of inclusion. The C22:6n-3 contents of the FFO, EFO, and MAG diets were equal. A control (CON) diet using blended vegetable oil was also made. As-hatched 1-d-old Ross 308 broilers (144) were reared (21 d) on a common starter diet then allocated to treatment pens (4 pens per treatment, 6 birds per pen) and fed treatment diets for 21 d before being slaughtered. Breast and leg meat was analyzed (per pen) for fatty acids, and cooked samples (2 pens per treatment) were analyzed for volatile aldehydes. Concentrations (mg/100 g of meat) of C20:5n-3, C22:5n-3, and C22:6n-3 were (respectively) CON: 4, 15, 24; FFO: 31, 46, 129; EFO: 18, 27, 122; LAG: 9, 19, 111; MAG: 6, 16, 147; and HAG: 9, 14, 187 (SEM: 2.4, 3.6, 13.1) in breast meat and CON: 4, 12, 9; FFO: 58, 56, 132; EFO: 63, 49, 153; LAG: 13, 14, 101; MAG: 11, 15, 102; HAG: 37, 37, 203 (SEM: 7.8, 6.7, 14.4) in leg meat. Cooked EFO and HAG leg meat was more oxidized (5.2 mg of hexanal/kg of meat) than the other meats (mean 2.2 mg/kg, SEM 0.63). It is concluded that algal biomass is as effective as fish oil at enriching broiler diets with C22:6 LC n-3 PUFA, and at equal C22:6n-3 contents, there is no significant difference between these 2 supplements on the oxidative stability of the meat that is produced.

Using Doppler radar with a simple explicit microphysics model to diagnose problems with ice sublimation depth-scales in forecast models

Several previous studies have attempted to assess the sublimation depth-scales of ice particles from clouds into clear air. Upon examining the sublimation depth-scales in the Met Office Unified Model (MetUM), it was found that the MetUM has evaporation depth-scales 2–3 times larger than radar observations. Similar results can be seen in the European Centre for Medium-Range Weather Forecasts (ECMWF), Regional Atmospheric Climate Model (RACMO) and Météo-France models. In this study, we use radar simulation (converting model variables into radar observations) and one-dimensional explicit microphysics numerical modelling to test and diagnose the cause of the deep sublimation depth-scales in the forecast model. The MetUM data and parametrization scheme are used to predict terminal velocity, which can be compared with the observed Doppler velocity. This can then be used to test the hypothesis as to why the sublimation depth-scale is too large within the MetUM. Turbulence could lead to dry air entrainment and higher evaporation rates; particle density may be wrong, particle capacitance may be too high and lead to incorrect evaporation rates or the humidity within the sublimating layer may be incorrectly represented. We show that the most likely cause of deep sublimation zones is an incorrect representation of model humidity in the layer. This is tested further by using a one-dimensional explicit microphysics model, which tests the sensitivity of ice sublimation to key atmospheric variables and is capable of including sonde and radar measurements to simulate real cases. Results suggest that the MetUM grid resolution at ice cloud altitudes is not sufficient enough to maintain the sharp drop in humidity that is observed in the sublimation zone.

The effect of the year of wheat variety release on productivity and stability of performance on two organic and two non-organic farms

Nineteen wheat cultivars, released from 1934 to 2000, were grown at two organic and two non-organic sites in each of 3 years. Assessments included grain yield, grain protein concentration, protein yield, disease incidence and green leaf area. The superiority of each cultivar (the sum of the squares of the differences between its mean in each environment and the mean of the best cultivar there, divided by twice the number of environments; CS) was calculated for yield, grain protein concentration and protein yield, and ranked in each environment. The yield and grain protein concentration CS were more closely correlated with cultivar release date at the non-organic sites than at organic sites. This difference may be attributed to higher yield levels with larger differences among cultivars at the non-organic sites, rather than to improved stability (i.e. similar ranks) across sites. The significant difference in the correlation of protein yield CS and cultivar age between organic and non-organic sites would support evidence that the ability to take up mineral nitrogen (N) compared to soil N has been a component of the selection conditions of more modern cultivars (released after 1989). This is supported by assessment of green leaf area (GLA), where more modern cultivars in the non-organic systems had greater late-season GLA, a trend that was not identified in organic conditions. This effect could explain the poor correlation between age and protein yield CS in organic compared to non-organic conditions where modern cultivars are selected to benefit from later nitrogen (N) availability which includes the spring nitrogen applications tailored to coincide with peak crop demand. Under organic management, N release is largely based on the breakdown of fertility-building crops incorporated (ploughed-in) in the previous autumn. The release of nutrients from these residues is dependent on the soil conditions, which includes temperature and microbial populations, in addition to the potential leaching effect of high winter rainfall in the UK. In organic cereal crops, early resource capture is a major advantage for maximizing the utilization of nutrients from residue breakdown. It is concluded that selection of cultivars under conditions of high agrochemical inputs selects for cultivars that yield well under maximal conditions in terms of nutrient availability and pest, disease and weed control. The selection conditions for breeding have a tendency to select cultivars which perform relatively better in non-organic compared to organic systems.

A hybrid data assimilation scheme for model parameter estimation: application to morphodynamic modelling

We present a novel algorithm for joint state-parameter estimation using sequential three dimensional variational data assimilation (3D Var) and demonstrate its application in the context of morphodynamic modelling using an idealised two parameter 1D sediment transport model. The new scheme combines a static representation of the state background error covariances with a flow dependent approximation of the state-parameter cross-covariances. For the case presented here, this involves calculating a local finite difference approximation of the gradient of the model with respect to the parameters. The new method is easy to implement and computationally inexpensive to run. Experimental results are positive with the scheme able to recover the model parameters to a high level of accuracy. We expect that there is potential for successful application of this new methodology to larger, more realistic models with more complex parameterisations.