Large-Scale Scientific Computing: 7th International Conference, LSSC 2009, Sozopol, Bulgaria, June 4-8, 2009 Revised Papers (Lecture Notes in Computer Science, 5910)

The 7th International Conference on Large-Scale Scienti?c Computations (LSSC 2009) was held in Sozopol, Bulgaria, June 4–8, 2009. The conference was organized and sponsored by the Institute for Parallel Processing at the B- garian Academy of Sciences. The conference was devoted to the 70th birthday anniversary of Professor Zahari Zlatev. The Bulgarian Academy of Sciences awarded him the Marin Drinov medal on ribbon for his outstanding results in environmental mat- matics and for his contributions to the Bulgarian mathematical society and the Academy of Sciences. The plenary invited speakers and lectures were: – P. Arbenz, “?Finite Element Analysis of Human Bone Structures” – Y. Efendiev, “Mixed Multiscale Finite Element Methods Using Limited Global Information” – U. Langer, “Fast Solvers for Non-Linear Time-Harmonic Problems” – T. Manteu?el, “First-Order System Least-Squares Approach to Resistive Magnetohydrodynamic Equations” – K. Sabelfeld, “Stochastic Simulation for Solving Random Boundary Value Problems and Some Applications” – F. Tro ̈ltzsch,“OnFinite ElementErrorEstimatesforOptimalControlPr- lems with Elliptic PDEs” – Z. Zlatev, “On Some Stability Properties of the Richardson Extrapolation Applied Together with the ?-method” The success of the conference and the present volume in particular are an outcome of the joint e?orts of many partnersfrom various institutions and or- nizations. Firstwe wouldlike to thank allthe membersofthe Scienti?c Comm- tee for their valuable contribution forming the scienti?c face of the conference, as well as for their help in reviewing contributed papers. We especially thank the organizers of the special sessions. Title Page Preface Table of Contents I Plenary and Invited Papers An Efficiency-Based Adaptive Refinement Scheme Applied to Incompressible, Resistive Magnetohydrodynamics Introduction The MHD Equations and FOSLS Formulation Solution Algorithm Adaptive Refinement Numerical Results Test Problem: Island Coalescence Discussion References Discontinuous Galerkin Subgrid Finite Element Method for Heterogeneous Brinkman’s Equations Introduction Preliminaries Subgrid Method for Brinkman’s Equations  Subgrid Method and Alternating Schwarz Iterations Numerical Results References Stochastic Simulation for Solving Random Boundary Value Problems and Some Applications Introduction Spectral Representations of Homogeneous Random Fields Random Walks and Double Randomization Method Random Loads Respond to Boundary Excitations Inhomogeneous Random Fields and Karhunen-Lo`eve Expansion Fractal Density of the Universe References On Finite Element Error Estimates for Optimal Control Problems with Elliptic PDEs Introduction Optimal Control Problem with Control Constraints The Unconstrained Optimal Control Problem Constraints on the Control Problems with State Constraints Available Error Estimates Control in IR$_m$ and State Constraints in Finitely Many Points References On Some Stability Properties of the Richardson Extrapolation Applied Together with the $\theta$-Method Richardson Extrapolation Combining the Richardson Extrapolation and the $\theta$-Method Studying the Stability Properties Introduction of an Atmospheric Chemistry Scheme Organization of the Numerical Experiments Results from the Numerical Experiments Achieving Second Order of Accuracy Comparing Results Obtained by the Backward Euler Formula and the $\theta$-Method Comparing Computing Times Needed to Obtain Prescribed Accuracy Concluding Remark References II Multilevel and Multiscale Preconditioning Methods On the Use of Aggregation-Based Parallel Multilevel Preconditioners in the LES of Wall-Bounded Turbulent Flows Introduction Pressure Systems in the LES of Turbulent Channel Flows Algebraic Multilevel Domain Decomposition Preconditioners Numerical Experiments Conclusions and Future Work References An Additive Matrix Preconditioning Method with Application for Domain Decomposition and Two-Level Matrix Partitionings Introduction Error Matrix for an Additive Preconditioning Method Condition Number Bound for a Perturbed Preconditioner Concluding Remarks References Numerical Study of AMLI Methods for Weighted Graph-Laplacians Introduction AMLI Method Preconditioner Construction for Weighted Graph-Laplacian Numerical Results Concluding Remarks References A Scalable TFETI Based Algorithm for 2D and 3D Frictionless Contact Problems Introduction TFETI and Contact Problems Optimality Numerical Experiments Comments and Conclusions References Multilevel Preconditioning of Crouzeix-Raviart 3D Pure Displacement Elasticity Problems Introduction Composite FR Algorithm Numerical Tests Concluding Remarks References Element-by-Element Schur Complement Approximations for General Nonsymmetric Matrices of Two-by-Two Block Form Introduction The Physical Problem and Mathematical Models Used Preconditioning Strategy and Estimates Numerical Illustrations Conclusions and Some Open Problems References Numerical Simulation of Fluid-Structure Interaction Problems on Hybrid Meshes with Algebraic Multigrid Methods Problem Setting of the Fluid-Structure Interaction Geometrical Description The Physical Model in Strong Form Weak Formulations Discretization Iterative Methods for the Interface Equation Numerical Results References III Multilevel and Multiscale Methods for Industrial Applications Recent Developments in the Multi-Scale-Finite-Volume Procedure Introduction The IMSFV-Procedure The Galerkin Approach Numerical Tests Conclusions References Boundary Element Simulation of Linear Water Waves in a Model Basin Introduction Modeling The Boundary Element Method Data-Sparse Approximation Numerical Results Conclusion References Numerical Homogenization of Bone Microstructure Introduction Locking-Free Nonconforming FEM Discretization Numerical Homogenization NumericalExperiments Model Problem Homogenization of Bone Microstructures References Multiscale Modeling and Simulation of Fluid Flows in Highly Deformable Porous Media Introduction Multilevel Algorithm Numerical Examples References IV Environmental Modelling Assimilation of Chemical Ground Measurements in Air Quality Modeling Introduction Methodology Background Fields Algorithms Experiment Setup Results and Discussions References Numerical Simulations with Data Assimilation Using an Adaptive POD Procedure Introduction Numerical Model and Adaptive POD 4D-Var Assimilation Procedure Numerical Results Conclusions References Game-Method Model for Field Fires Introduction Description of the Game-Method Model Experimental Results Conclusion References Joint Analysis of Regional Scale Transport and Transformation of Air Pollution from Road and Ship Transport Introduction Modelling Tools Methodology Results and Discussions Conclusion References Runs of UNI–DEM Model on IBM Blue Gene/P Computer and Analysis of the Model Performance Introduction Mathematical Description of the Danish Eulerian Model Need for Efficient Implementation of the Computations Parallel Implementation of DEM IBM Blue Gene/P Computer Parallelization Strategy Numerical Results Conclusions References Sensitivity Analysis of a Large-Scale Air Pollution Model: Numerical Aspects and a Highly Parallel Implementation Introduction Sensitivity Analysis Concept Definition of the Total Sensitivity Indices (TSI) Sobol’ Approach, Based on HDMR and ANOVA The Danish Eulerian Model DEM, UNI-DEM, and SA-DEM Performance Results Sensitivity Analysis Experiments Conclusions and Plans for Future Work References Advanced Results of the PM10 and PM2.5 Winter 2003 Episode by Using MM5-CMAQ and WRF/CHEM Models Introduction PM10 and PM2.5 Episode Emission Data MM5-CMAQ and WRF-CHEM Architectures and Configurations Changes in Model Configurations Model Results Conclusions References Four-Dimensional Variational Assimilation of Atmospheric Chemical Data – Application to Regional Modelling of Air Quality Introduction Four-Dimensional Variational Data Assimilation The EURAD-IM Chemistry Transport Model Applications and Results Nested 4d-Var Application to the BERLIOZ Campaign Operational 3d-Var Application for PROMOTE Summary References Numerical Study of Some High PM10 Levels Episodes Introduction Modeling Tools Model Configuration and Brief Description of the Simulations Comparison of the PM10 Simulations with Measurments Process Analysis Conclusions References V Control and Uncertain Systems Hausdorff Continuous Viscosity Solutions of Hamilton-Jacobi Equations Introduction The Space of Hausdorff Continuous Functions The Envelope Viscosity Solutions and Hausdorff Continuous Viscosity Solutions Existence and Uniqueness Conclusion References Stochastic Skiba Sets: An Example from Models of Illicit Drug Consumption Introduction The Model Analysis Discussion, Extensions, and Concluding Remarks References Classical and Relaxed Optimization Methods for Nonlinear Parabolic Optimal Control Problems Classical and Relaxed Optimal Control Problems Classical and Relaxed Optimization Methods Numerical Examples References Exponential Formula for Impulsive Differential Inclusions Preliminaries Euler Approximation of the Solution Set L$_1$–Metric for Piecewise Lipschitz Functions and of the Reachable Sets of Impulsive ODI References Directional Sensitivity Differentials for Parametric Bang-Bang Control Problems Introduction Regular Bang-Bang Controls Strong Stability of Bang-Bang Structure Second-Order Condition Directional Sensitivity Differentials References Estimates of Trajectory Tubes of Uncertain Nonlinear Control Systems Introduction Preliminaries Basic Notations Trajectory Tubes in Uncertain Systems with Quadratic Nonlinearity Results Differential System with Uncertainty in Initial States Control System under Uncertainty Conclusion References Asymptotics for Singularly Perturbed Reachable Sets Problem Statement Splitting of Dynamic System Asymptotics of the Support Functions to Reachable Sets Concluding Remarks References On Optimal Control Problem for the Bundle of Trajectories of Uncertain System Introduction Linear Models Input Design Procedure for Nonlinear Control System Input Synthesis Procedure References High-Order Approximations to Nonholonomic Affine Control Systems Introduction A Second Order Approximation: Basic Idea An Auxiliary Problem A Constructive Second Order Scheme for Two Non-commutative Fields References VI Application of Metaheuristics to Large Scale Problems A Parametric Multi-start Algorithm for Solving the Response Time Variability Problem Introduction The Response Time Variability Problem (RTVP) The Multi-start Algorithm for Solving the RTVP Computational Experiment Final Conclusions References Enhancing the Scalability of Metaheuristics by Cooperative Coevolution Introduction Cooperative Coevolution for Large Scale Optimization Cooperative Coevolutionary Differential Evolution and Harmony Search Numerical Results Conclusions References Hybrid ACO Algorithm for the GPS Surveying Problem Introduction Problem Description Hybrid ACO Algorithm Experimental Results Conclusions and Future Work References Generalized Nets as Tools for Modeling of the Ant Colony Optimization Algorithms Introduction Short Remarks on Generalized Nets GNs as Tools for Modelling of ACO Algorithms A GN-Model of ACO for Multiple Knapsack Problem Conclusion References A Scatter Search Approach for Solving the Automatic Cell Planning Problem Introduction The ACP Problem Setting Up AbYSS Solution Representation The Template Methods Experimentation ACP Instances Algorithm Settings Methodology Results Conclusions and Future Work References Some Aspects Regarding the Application of the Ant Colony Meta-heuristic to Scheduling Problems Introduction Resource-Constrained Project Scheduling Problem Solution ACO (Ant Colony Optimization) Building the Solution Updating the Pheromone Dynamic Evaporation Approach Results Conclusion References High-Performance Heuristics for Optimization in Stochastic Traffic Engineering Problems Introduction Problem Statement and Preliminaries Performance of Exact Algorithm High-Performance Heuristics Conclusion References Optimization of Complex SVM Kernels Using a Hybrid Algorithm Based on Wasp Behaviour Introduction Kernel Substitution Method The Co-mutation $LR$ − $M_ijn$ Operator Wasp Behaviour Computational Models Main Results: The Model for Constructing Complex SVM Kernels Basic Construction Optimization of the $LR$ − $M_ijn$ Operator Using a Wasp-Based Computational Scheme SVM Algorithm Experimental Results Conclusions and Further Directions of Study References VII Monte Carlo: Methods, Applications, Distributed Computing Enabling Cutting-Edge Semiconductor Simulation through Grid Technology Introduction to NanoCMOS Challenges NanoCMOS E-Infrastructure and Simulation Methodology NanoCMOS Scientific Results and Discussion Conclusions and Future Work References Incremental Reuse of Paths in Random Walk Radiosity Introduction Previous Work: Reuse of Shooting Paths in Radiosity Our Contribution: Incremental Reuse of the Shooting Paths Applications Results Conclusions and Future Work References Monte Carlo Adaptive Technique for Sensitivity Analysis of a Large-Scale Air Pollution Model Introduction Background Studies Sobol’ Global Sensitivity Indices Concept Monte Carlo Approach for Small Sensitivity Indices Description of the Algorithms Analysis of Numerical Results and Discussion References Parallel Implementation of the Stochastic Radiosity Method Introduction Stochastic Iteration Method Sequential Implementation Parallel Implementation Results Conclusions References Monte-Carlo Modeling of Electron Kinetics in Room Temperature Quantum-Dot Photodetectors Introduction Simulation Models Results of Simulation Analytical Results Conclusion References Particle Model of the Scattering-Induced Wigner Function Correction Introduction Scattering-Induced Wigner Function Correction Classical Limit Monte Carlo Approach Simulation Results References Analysis of the Monte Carlo Image Creation by Uniform Separation Introduction Uniform Separation of Integration Domain Error Analysis Conclusions References The Role of the Boundary Conditions on the Current Degradation in FD-SOI Devices Introduction Convergence of the Scheme The Role of the Thermal Boundary Conditions Non-stationary Electron and Phonon Transports Conclusions References Gamma Photon Transport on the GPU for PET Introduction New Algorithm Free Path Sampling Termination Handling Scattering Direction Deterministic Detector Response Calculation Results Conclusions References Transport in Nanostructures: A Comparative Analysis Using Monte Carlo Simulation, the Spherical Harmonic Method, and Higher Moments Models Introduction Simulation Methods Monte Carlo Solution of the Boltzmann Equation Higher Moments Models Spherical Harmonic Expansion Results Conclusion References Thermal Modeling of GaN HEMTs Theoretical Model Simulation Results Conclusions References VIII Grid and Scientific and Engineering Applications Tuning the Generation of Sobol Sequence with Owen Scrambling Introduction Preliminaries Sobol Sequence Owen Scrambling Algorithm for GPU Implementation and Computational Tests The Graphics Processing Unit The Algorithm Numerical Tests Conclusions References Applying the Improved Saleve Framework for Modeling Abrasion of Pebbles Introduction Grid Authentication Sequence in Saleve Authentication towards the gLite Middleware Involving the Client in the Delegation Process Simulating the Abrasion of Pebbles Performance of the Application Concluding Remarks References Information Flow and Mirroring in an Agent-Based Grid Resource Brokering System Introduction SystemOverview Need for Information Mirroring What Should Be Mirrored, When and Where? Team Member Data Job Contracts and Their Execution Trust Information Volume Data Collection and Storage Concluding Remarks References Scatter Search and Grid Computing to Improve Nuclear Fusion Devices Introduction NuclearFusion Scatter Search Algorithm Concurrent Implementation Workflow to Measure the Equilibrium Measuring the Equilibrium Workflow Application Grid-Oriented SS Implementation Python Scripts SS Implementation Results Computational Results Evaluation Results Conclusions and Future Work References Service-Oriented Integration of Grid Applications in Heterogeneous Grids Introduction Problem Definition Service-Oriented Approach for Integration Applying Approach Experience Conclusion and Future Work References Grid Based Environment Application Development Methodology Introduction Related Works Grid Application Development Methodology Algorithm Identification and Analysis Data Model Definition Identify Atomic Parts of the Computation Algorithm Implementation Graph Based Process Description Pattern Process Description Instantiated Process Description Graph Process Execution User Interface Development Conclusions References User Level Grid Quality of Service Introduction The Need for Grid QoS Current Grid QoS User Level QoS System ULQoSOperation Usage Results Future Work Conclusion References Dictionary Compression and Information Source Correction Introduction Word Group Compression Storage Formats Binary Non-indexed Tree Representation Binary Indexed Tree Representation Statistical Analysis Word Group Cutting DAG Transformation Towards Global Optimality Conclusions References A Parallelization of Finite Volume Method for Calculation of Gas Microflows by Domain Decomposition Methods Introduction Parallel Organization Speedup Analysis Conclusions References Background Pollution Forecast over Bulgaria Introduction Operational Design of BGCW Emission Modelling CMAQ Calculations and Validation Conclusion References Climate Change Impact Assessment of Air Pollution Levels in Bulgaria Introduction Meteorological Data Base Air Quality Models Used, Computational Domains and Data Flow Modelling Validation of Model Results Long Term Simulations: Results for 1991–2000 Time–Slice Conclusion References Effective Algorithm for Calculating Spatial Deformations of Pre-stressed Concrete Beams Introduction Description The Problem Solution Testing Potential Algorithm Analysis Concept Computation Need of Gradient Method Computation Need of Simplex Decomposition Memory Need Parallelization Results Conclusion References IX Reliable Numerical Methods for Differential Equations Structurally Stable Numerical Schemes for Applied Dynamical Models Introduction Topological Dynamic Consistency SEIR(→S) Model $\mathcal{V}$-Structurally Stable Scheme by the Nonstandard Finite Difference Method References Matrix and Discrete Maximum Principles Introduction The Discrete Heat Problem and the Iterative Algebraic Problems Matrix Maximum Principles and Their Relations MMP to Linear Algebraic Systems of Special Form Qualitative Properties of IAP References On a Bisection Algorithm That Produces Conforming Locally Refined Simplicial Meshes Introduction Mesh Density Function Convergence of the Mesh-Size for GCB-Algorithm Numerical Results References A Discrete Maximum Principle for Nonlinear Elliptic Systems with Interface Conditions Introduction Some Algebraic Background The Problem Finite Element Discretization The Discrete Maximum Principle Some Related Results References Inverse Problem for Coefficient Identification in Euler-Bernoulli Equation by Linear Spline Approximation Introduction Inverse Problem Formulation Variational Imbedding Equation for $u$ Equation for $\sigma$ Difference Scheme Grid Pattern and Approximations General Construction of the Algorithm NumericalExperiments Linear Coefficient Linear Coefficient as a Piecewise Linear Function References A Fully Implicit Method for Fluid Flow Based on Vectorial Operator Splitting Introduction Problem Statement Incompressible Navier–Stokes Equations Initial and Boundary Conditions Singularity of Direct Fully Implicit Schemes Balanced Pressure Equation Difference Problem Analytical Requirements Coordinate Splitting and Spatial Discretization Numerical Results Conclusion References Discrete Maximum Principle for Finite Element Parabolic Operators Introduction Mesh Operators, the Discrete Boundary Maximum Principles and the Non-negativity Preservation Property Conditions for DWBMP/DSBMP of $\mathcal{L}_{a,b,c}$ Local Conditions for DWBMP/DSBMP of $\mathcal{L}_{a,b,c}$ Global Conditions for DWBMP/DSBMP of $\mathcal{L}_{a,b,c}$ Numerical Examples References Climate Change Scenarios for Hungary Based on Numerical Simulations with a Dynamical Climate Model Introduction Regional Climate Model PRECIS Design of the Statistical Analysis Results Conclusions References Numerical Simulations of Reaction-Diffusion Systems Arising in Chemistry Using Exponential Integrators Introduction General Framework of the Exponential Integrators Defining the Model Equations Numerical Schemes and Computational Issues Conclusions References On the Discretization Time-Step in the Finite Element Theta-Method of the Two-Dimensional Discrete Heat Equation Preliminaries Analysis of FEM Equation Numerical Experiments Conclusions and Further Works References X Novel Applications of Optimization Ideas to the Numerical Solution of PDEs A Locally Conservative Mimetic Least-Squares Finite Element Method for the Stokes Equations Introduction Notation Mimetic LSFEM for the Stokes Equations Conservative Properties References Additive Operator Decomposition and Optimization–Based Reconnection with Applications Introduction Additive Operator Decomposition Reconnection via an Optimization Formulation Existence of Optimal Solutions Reformulation Error A Solution Algorithm Numerical Results References Least-Squares Spectral Element Method on a Staggered Grid Mimetic Formulation of the Poisson Equation Dual Spectral Element Meshes The Reduction Operators $\mathcal{R}$ and InterpolationOperators $\mathcal{I}$ The Conservation Law dq(n−1) = f(n) Discretization of u(1) = dφ(0) Results References Mimetic Least-Squares Spectral/hp Finite Element Method for the Poisson Equation Introduction Mimetic Approaches for the 2D Poisson Equation Differential Geometry, a Refresher Tonti Diagram for the 2D Poisson Equation and Its Discretizations Weak Material Laws: The Role of Least-Squares Application to the 2D Poisson Equation Concluding Remarks References Adaptive Least Squares Finite Element Methods in Elasto-Plasticity Introduction Least-Squares Formulation of IncrementalElasto-Plasticity Plane Stress Model and Finite Element Approximation Generalized Gauss-Newton Method A Computational Example References The Automatic Construction and Solution of a Partial Differential Equation from the Strong Form Introduction A Calculus for B ́ezier Patches and Surfaces Basic Definition Sums, Negation, Subtraction, and Multiplication Derivatives Smoothness Symbolic Coefficients Algorithm Decomposition into a First Order System Approximating Non-B ́ezier Surfaces Constructing the Least Squares Problem Numerical Experiments Conclusion References XI Contributed Papers Zienkiewicz-Type Finite Element Applied to Fourth-Order Problems Introduction Main Result Numerical Results References Acceleration of Convergence for Eigenpairs Approximated by Means of Non-conforming Finite Element Methods Introduction ConvergenceAnalysis Superconvergent Postprocessing Technique Numerical Results References A Two-Grid Method on Layer-Adapted Meshes for a Semilinear 2D Reaction-Diffusion Problem Introduction Meshes and Linear Problem Uniform Convergence for the Nonlinear Problem Two-Grid Algorithm Numerical Experiments References Comparative Analysis of Solution Methods for Delay Differential Equations in Hematology Introduction Mathematical Model Solution Methods and Software Tools Numerical Tests and Analysis of the Results References Solving Non-linear Systems of Equations on Graphics Processing Units Introduction Preliminaries Main Results Isometric Conversion between Dimension and Resolution: A Simpler Approach Computing the Solutions of Nonlinear Operator Equations Concluding Remarks and Simple Model Examples References Computing $n$-Variate Orthogonal Discrete Wavelet Transforms on Graphics Processing Units Introduction Main Result: A New Basis-Matching Algorithm Canonical Isometries between Hilbert Spaces Discrete Wavelet Transform Problem Setting Computing B$_{d,j}$ Concluding Remarks References Wavelet Compression, Data Fitting and Approximation Based on Adaptive Composition of Lorentz-Type Thresholding and Besov-Type Non-threshold Shrinkage Introduction Preliminaries Riesz Unconditional Wavelet Bases in Besov Spaces Non-parametric Regression Wavelet Shrinkage Composite Besov-Lorentz Shrinkage Besov-Lorentz Shrinkage versus Firm Thresholding References On Interpolation in the Unit Disk Based on Both Radon Projections and Function Values Introduction and Preliminaries Interpolation Problem for Radon Projections Type of Data Interpolation Problem for Mixed Type of Data Regular Schemes for Mixed Type of Data Numerical Experiments References Comparison between the Marching-Cube and the Marching-Simplex Methods Introduction A Hierarchy of “Marching” Algorithms for Intersection Marching Square/Cube Algorithm A Vertex-Based CPU-Computable Upgrade A Pixel-Based GPU-Computable Upgrade The Marching Simplex Algorithm Comparison of the Considered Four Marching Algorithms Future Work References Transitions from Static to Dynamic State in Three Stacked Josephson Junctions Introduction Mathematical Model Numerical Methods and Algorithms Numerical Results References Factorization-Based Graph Repartitionings Introduction Basic Terminology and Our Restrictions Framework for the Graph-Based Repartitioning Removal of Vertices Insertion of Vertices Repartitioning for Generalized Decompositions Conclusion References An Efficient Algorithm for Bilinear Strict Equivalent (BSE)- Matrix Pencils Introduction Creating the Conditions for the BSE The Algorithm for the BSE Numerical Examples Conclusion-Further Research References Two-Grid Decoupling Method for Elliptic Problems on Disjoint Domains Introduction Two-Rectangles Model Problem Finite Element Method Solution Two-Grid Decoupling FEM Algorithms Numerical Examples Conclusions References A Method for Sparse-Matrix Computation of B-Spline Curves and Surfaces Introduction Bezier Curves and the De Casteljau Algorithm The B-Spline Factor Matrix T(t) Commutativity Relations between T(t) Matrices and Their Derivatives B-Splines in Matrix Notation Knot Insertion The B-Spline Matrix Complete Hermite Interpolation with B-Splines References Parallel MIC(0) Preconditioning for Numerical Upscaling of Anisotropic Linear Elastic Materials Introduction Homogenization Scheme Parallel MIC(0) Preconditioning NumericalExperiments References The Bpmpd Interior Point Solver for Convex Quadratically Constrained Quadratic Programming Problems Introduction The Infeasible Primal–Dual Method Implementation Improving the Numerical Properties of Factorizations Implementation of the Numerical Kernels Sparse Matrix Ordering Numerical Results and Conclusions References On Shape Optimization of Acoustically Driven Microfluidic Biochips Introduction The Physical Model Shape Optimization Problem References Block Preconditioners for the Incompressible Stokes Problem Introduction Block Preconditioners SIMPLE-Type Preconditioners LSC Preconditioner NumericalExperiments Isoviscos Problem Variable Viscosity Problem Conclusions References Author Index