Experience Seamless Full Ebook Downloads for Every Genre at textbookfull.

com

Fundamentals of RCS Prediction Methodology using

Parallelized Numerical Electromagnetics Code NEC
and Finite Element Pre processor Vineetha Joy

https://textbookfull.com/product/fundamentals-of-rcs-
prediction-methodology-using-parallelized-numerical-
electromagnetics-code-nec-and-finite-element-pre-processor-
vineetha-joy/

OR CLICK BUTTON

DOWNLOAD NOW

Explore and download more ebook at https://textbookfull.com

 Recommended digital products (PDF, EPUB, MOBI) that
you can download immediately if you are interested.

Electromagnetics through the Finite Element Method: A

Simplified Approach Using Maxwell's Equations 1st Edition
José Roberto Cardoso
https://textbookfull.com/product/electromagnetics-through-the-finite-
element-method-a-simplified-approach-using-maxwells-equations-1st-
edition-jose-roberto-cardoso/
textboxfull.com

Numerical Solution of Differential Equations: Introduction

to Finite Difference and Finite Element Methods 1st
Edition Zhilin Li
https://textbookfull.com/product/numerical-solution-of-differential-
equations-introduction-to-finite-difference-and-finite-element-
methods-1st-edition-zhilin-li/
textboxfull.com

Harmonic Balance Finite Element Method: Applications in

Nonlinear Electromagnetics and Power Systems 1st Edition
Junwei Lu
https://textbookfull.com/product/harmonic-balance-finite-element-
method-applications-in-nonlinear-electromagnetics-and-power-
systems-1st-edition-junwei-lu/
textboxfull.com

Fundamentals of the Finite Element Method for Heat and

Mass Transfer 2nd Edition Perumal Nithiarasu

https://textbookfull.com/product/fundamentals-of-the-finite-element-
method-for-heat-and-mass-transfer-2nd-edition-perumal-nithiarasu/

textboxfull.com
Finite Element Analysis of Weld Thermal Cycles Using ANSYS
1st Edition G. Ravichandran (Author)

https://textbookfull.com/product/finite-element-analysis-of-weld-
thermal-cycles-using-ansys-1st-edition-g-ravichandran-author/

textboxfull.com

The finite element method: fundamentals and applications

in civil, hydraulic, mechanical and aeronautical
engineering Zhu B.
https://textbookfull.com/product/the-finite-element-method-
fundamentals-and-applications-in-civil-hydraulic-mechanical-and-
aeronautical-engineering-zhu-b/
textboxfull.com

Extended finite element and meshfree methods Timon Rabczuk

https://textbookfull.com/product/extended-finite-element-and-meshfree-
methods-timon-rabczuk/

textboxfull.com

Finite Element Analysis A Primer 2nd Edition Musa

https://textbookfull.com/product/finite-element-analysis-a-primer-2nd-
edition-musa/

textboxfull.com

Numerical Methods for Engineering: An introduction using

MATLAB® and computational electromagnetics examples
(Electromagnetic Waves) 2nd Edition Karl F. Warnick
https://textbookfull.com/product/numerical-methods-for-engineering-an-
introduction-using-matlab-and-computational-electromagnetics-examples-
electromagnetic-waves-2nd-edition-karl-f-warnick/
textboxfull.com
SPRINGER BRIEFS IN ELEC TRIC AL AND COMPUTER
ENGINEERING  COMPUTATIONAL ELECTROMAGNETICS

Vineetha Joy
G. L. Rajeshwari
Hema Singh
Raveendranath U. Nair

Fundamentals
of RCS Prediction
Methodology using
Parallelized Numerical
Electromagnetics
Code (NEC) and Finite
Element Pre-processor

123
SpringerBriefs in Electrical and Computer
Engineering

SpringerBriefs in Computational Electromagnetics

Series Editors
K. J. Vinoy, Indian Institute of Science, Professor, Electrical Communication
Engineering, Indian Institute of Science, Bangalore, Karnataka, India
Rakesh Mohan Jha (Late), Centre for Electromagnetics, CSIR-National Aerospace
Laboratories, Bangalore, Karnataka, India
The phenomenon of electric and magnetic field vector (wave) propagation through the
free-space, or any other medium is considered within the ambit of electromagnetics. The
media themselves, in general, could be of diverse type, such as linear/non-linear,
isotropic/non-isotropic, homogeneous/inhomogeneous, reciprocal/non-reciprocal, etc. Such
electromagnetic wave propagation problems are formulated with the set of Maxwell’s
equations. Computational Electromagnetics endeavors to provide the solution to the
Maxwell’s equations for a given formulation. It is often difficult to find closed forms
solutions to the Maxwell’s equation formulations. The advent of computers, and particularly
the initial developments of efficient coding for numerical analysis, encouraged the
development of numerical electromagnetics. A second motivation came from the interaction
of the electromagnetic wave with the matter. This could be visualized as scattering bodies,
which required incorporation of the phenomena of reflection, refraction, diffraction and
polarization. The finite/large nature of the scatterer required that problem of electromagnetics
is considered with respect to the operation wavelength leading to the classification of
low-frequency, high-frequency and resonance region problems. This also inspired various
asymptotic and grid-based finite-method techniques, for solving specific electromagnetic
problems. Surface modeling and ray tracing are also considered for such electromagnetic
problems. Further, design optimization towards hardware realization have led to the recourse
to various soft computing algorithms. Computational Electromagnetics is deemed to
encompass the numerical electromagnetics along with all other above developments. With
the wide availability of massively parallel high performance parallel computing platforms,
new possibilities have emerged for reducing the computation time and developing macro
models that can even be employed for several practical multi-physics scenarios. Both volume
and surface discretization methods have been given a new boost, and several acceleration
techniques including GPU based computation, learning based approaches, and model order
reduction have been attempted. Limitations of generating meshes and modifying these for
parametric estimation have been addressed by statistical approaches and smart solvers. Many
nature-inspired algorithms and other soft computing approaches have been employed for
electromagnetic synthesis problems. One of the recent additions is Game Theoretic
optimization.
Finally, the emergence of Computational Electromagnetics has been motivated by myriad
applications. These diverse application include but are not restricted to those in Electronics
and Communication, Wireless Propagation, Computer Hardware, Aerospace Engineering,
Biomedical Engineering, Radio-astronomy, Terahertz Technology, Photonics, etc for
modelling of devices, components, systems and even large structures. Some of the
well-known applications are in Analysis and design of radio frequency (RF) circuit, antennas
and systems, Analysis of antenna on structures, radar imaging, radar cross section
(RCS) computation and reduction, and analysis of electromagnetic wave-matter interactions
at discrete, random and periodic geometries including metamaterials. Authors are encouraged
to submit original research work in the area of Computational Electromagnetics. The content
could be either theoretical development, or specific to particular applications. This Series also
encourages state-of-the-art reviews and easy to comprehend tutorials.

More information about this subseries at http://www.springer.com/series/13885

Vineetha Joy G. L. Rajeshwari
• •

Hema Singh Raveendranath U. Nair

•

Fundamentals of RCS
Prediction Methodology
using Parallelized Numerical
Electromagnetics Code
(NEC) and Finite Element
Pre-processor

123
Vineetha Joy G. L. Rajeshwari
Centre for Electromagnetics (CEM) Centre for Electromagnetics (CEM)
CSIR-National Aerospace Laboratories CSIR-National Aerospace Laboratories
Bengaluru, Karnataka, India Bengaluru, Karnataka, India

Hema Singh Raveendranath U. Nair

Centre for Electromagnetics (CEM) Centre for Electromagnetics (CEM)
CSIR-National Aerospace Laboratories CSIR-National Aerospace Laboratories
Bengaluru, Karnataka, India Bengaluru, Karnataka, India

ISSN 2191-8112 ISSN 2191-8120 (electronic)

SpringerBriefs in Electrical and Computer Engineering
ISSN 2365-6239 ISSN 2365-6247 (electronic)
SpringerBriefs in Computational Electromagnetics
ISBN 978-981-15-7163-3 ISBN 978-981-15-7164-0 (eBook)
https://doi.org/10.1007/978-981-15-7164-0
© The Author(s), under exclusive license to Springer Nature Singapore Pte Ltd. 2021
This work is subject to copyright. All rights are solely and exclusively licensed by the Publisher, whether
the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of
illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and
transmission or information storage and retrieval, electronic adaptation, computer software, or by similar
or dissimilar methodology now known or hereafter developed.
The use of general descriptive names, registered names, trademarks, service marks, etc. in this
publication does not imply, even in the absence of a specific statement, that such names are exempt from
the relevant protective laws and regulations and therefore free for general use.
The publisher, the authors and the editors are safe to assume that the advice and information in this
book are believed to be true and accurate at the date of publication. Neither the publisher nor the
authors or the editors give a warranty, expressed or implied, with respect to the material contained
herein or for any errors or omissions that may have been made. The publisher remains neutral with regard
to jurisdictional claims in published maps and institutional affiliations.

This Springer imprint is published by the registered company Springer Nature Singapore Pte Ltd.
The registered company address is: 152 Beach Road, #21-01/04 Gateway East, Singapore 189721,
Singapore
To
Our Families
Preface

With the advent of stealth technology, precise computation of radar cross section
(RCS) has become an inevitable component in the efficient design and development
of military vehicles. However, the scattering characterization of stealth platforms is
undoubtedly the most demanding problem in modern applied electromagnetics.
When it comes to the computation of scattered fields from conducting structures,
the accuracy and efficacy of a method of moments (MoM)-based integral equation
formulation are difficult to be surpassed. Among MoM-based solvers, Numerical
Electromagnetics Code (NEC) is a versatile open-source computer program used
for the electromagnetic analysis of metallic structures in the presence of sources or
incident fields. Since its development, NEC has continued to be one of the most
widely used electromagnetic simulation codes even in the presence of commercial
MoM-based solvers. This can be attributed to the availability of the well-
documented computational engine. However, a serious drawback of NEC is the
absence of an appropriate module for the wire-grid-based meshing of geometries
and the consequent generation of segmentation data in an NEC compatible format.
This is one of the most intricate and time-consuming steps in the computation of
RCS using NEC. In this regard, this brief presents a detailed methodology for the
computation of RCS of metallic structures using a parallelized version of NEC in
conjunction with a finite element pre-processor, which has been strategically
incorporated for geometry modelling catering to NEC guidelines. It includes a
thorough overview of the theoretical background of NEC including all relevant
aspects of formulation and modelling. This book will serve as a valuable resource
for students, researchers, scientists, and engineers working in the field of RCS
predictions and measurements.

Bengaluru, India Vineetha Joy

G. L. Rajeshwari
Hema Singh
Raveendranath U. Nair

vii
Acknowledgements

First and foremost, we would like to thank the God Almighty for being the beacon
of hope during the course of this research work.
Further, we would like to thank Shri. Jitendra J. Jadhav, Director, CSIR-National
Aerospace Laboratories, Bengaluru, for the permission to write this SpringerBrief.
We would also like to acknowledge the valuable suggestions from our col-
leagues at the Centre for Electromagnetics during the course of writing this book.
But for the concerted support and encouragement of Ms. Swati Meherishi,
Editorial Director, Applied Sciences and Engineering, and Ms. Muskan Jaiswal of
Springer Nature India Private Limited, it would not have been possible to bring out
this book within such a short span of time.
We are also forever beholden to our family members for their incessant inspi-
ration which supported us to stay at rough tracks. Vineetha Joy would like to
specially thank her husband Johnu George for his constant support and encour-
agment during the course of this work.

ix
Contents

1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
2 Theoretical Background for the Computation of Radar
Cross-Section (RCS) Using NEC . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.1 Formulation of Integral Equations (IE) . . . . . . . . . . . . . . . . . . . . . 5
2.1.1 Electric-Field Integral Equation (EF-IE) . . . . . . . . . . . . . . 5
2.1.2 Magnetic-Field Integral Equation (MF-IE) . . . . . . . . . . . . . 6
2.1.3 EF-IE/MF-IE Combined Equation . . . . . . . . . . . . . . . . . . . 7
2.2 Computation of Scattered Fields . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.3 Numerical Solution of Integral Equations . . . . . . . . . . . . . . . . . . . 9
2.3.1 Selection of Basis Functions and Weight
Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3.2 Solution of Matrix Equation . . . . . . . . . . . . . . . . . . . . . . . 11
2.4 Strategy for Parallelization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
3 Modelling Guidelines and Input/Output Formats in NEC . . . . . . . . 17
3.1 Modelling of Wires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.2 Modelling of Surfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.3 Format of NEC Input File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.4 Format of NEC Output File . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
Reference . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4 Methodology for the Computation of RCS using Parallelized NEC
and Finite Element Pre-processor . . . . . . . . . . . . . . . . . . . . . . . . . .. 23
4.1 Overview of Finite Element Pre-processor . . . . . . . . . . . . . . . . .. 23
4.2 Revised Methodology for the Computation of RCS using
Parallelized NEC and Finite Element Pre-processor . . . . . . . . . .. 25
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 27

xi
xii Contents

5 Studies on Pre-selected Structures . . . . . . . . . . . . . . . . . . . . . . . . . . 29

5.1 Monostatic RCS of a PEC Dihedral . . . . . . . . . . . . . . . . . . . . . . . 29
5.2 Monostatic RCS of NASA Almond . . . . . . . . . . . . . . . . . . . . . . . 30
5.3 Monostatic RCS of Metallic Single Ogive . . . . . . . . . . . . . . . . . . 33
5.4 Monostatic RCS of Rocket-Shaped Model . . . . . . . . . . . . . . . . . . 35
5.5 Monostatic RCS of Micro Aerial Vehicle (MAV) . . . . . . . . . . . . . 37
5.6 Bi-Static RCS of Unmanned Aerial Vehicle (UAV) . . . . . . . . . . . 39
5.7 Monostatic RCS of Automobile . . . . . . . . . . . . . . . . . . . . . . . . . . 40
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
6 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53

Appendix A: Libraries used for Parallel Computations . . . . . . . . . . . . . . 55

Appendix B: Format of Command Lines in NEC . . . . . . . . . . . . . . . . . . 61
Appendix C: Format of Various Output Files . . . . . . . . . . . . . . . . . . . . . 67
Author Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Subject Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
About the Authors

Mrs. Vineetha Joy is working as Scientist at Centre for Electromagnetics of

CSIR-National Aerospace Laboratories, Bangalore, India since March, 2016. She
obtained B. Tech in Electronics and Communication Engineering with Second
Rank from University of Calicut in 2011 and received M. Tech degree with First
Rank in RF & Microwave Engineering from Indian Institute of Technology (IIT),
Kharagpur in 2014. She is actively involved in various R&D programs spanning
domains like Computational Electromagnetics, Electromagnetic Design and
Performance Analysis of Radomes, Design and Development of broadband radar
absorbing structures (RAS), Hybrid Numerical Techniques for Scattering Analysis,
and Electromagnetic Characterization of Potential Materials for Airborne
Structures. She has authored/co-authored several research publications including
peer reviewed journal papers, symposium papers, technical documents and test
reports.

Ms. G. L. Rajeshwari obtained M. Tech. in Digital Electronics and Communi-

cation in 2019 from Department of Electronics & Communication Engineering,
Bangalore Institute of Technology, Bengaluru, India and B.E. in Electronics &
Communication Engineering, in 2017 from Nitte Meenakshi Institute of Tech-
nology, Bengaluru. Her research interest includes Computational Electromagnetics,
Hybrid Numerical Techniques for Scattering Analysis etc.

Dr. Hema Singh is working as Senior Principal Scientist in Centre for Electro-
magnetics, National Aerospace Laboratories (CSIR-NAL), Bangalore, India. She
has received Ph.D. degree in Electronics Engineering from IIT-BHU, Varanasi
India in Feb. 2000. For the period 1999-2001, she was Lecturer in Physics at P.G.
College, Kashipur, Uttaranchal, India. She was a Lecturer in EEE of Birla Institute
of Technology & Science (BITS), Pilani, Rajasthan, India, for the period
2001-2004. She joined CSIR-NAL as Scientist in January 2005. Her active areas of
research are Computational Electromagnetics for Aerospace Applications, EM
analysis of propagation in an indoor environment, Phased Arrays, Conformal
Antennas, Radar Cross Section (RCS) Studies including Active RCS Reduction.

xiii
xiv About the Authors

She has contributed in the projects not only sponsored by DRDO on low RCS
phased array, Active RCS reduction, RAS development, but also in the project
sponsored by Boeing USA on EM analysis of RF field build-up within Boeing 787
Dreamliner, She received Best Woman Scientist Award in CSIR-NAL, Bangalore
for period of 2007-2008 for her contribution in the area of active RCS reduction.
Dr. Singh co-authored 14 books, 2 book chapters, 7 software copyrights, 370
scientific research papers and technical reports.

Dr. Raveendranath U. Nair is working as Senior Principal Scientist & Head in

Centre for Electromagnetics, National Aerospace Laboratories (CSIR-NAL),
Bangalore, India. He holds a Ph.D. degree in Microwave Electronics from MG
University, India. He has over 20 years’ experience in the field of electromagnetic
design and analysis of radomes, design and development of FSS based structures
for airborne platforms, radar cross section (RCS) studies, design and development
of artificially engineered materials etc. He has contributed significantly to the
National Radome Programs including Doppler Weather Radar (DWR) radome
installed at SHAAR Sreeharikota, Fire Control Radar (FCR) radome for Jaguar
aircraft, Astra missile ceramic radome, nosecone radome for Saras aircraft, multi-
band radomes for TU-142M aircraft etc. He has authored/co-authored over 200
research publications including peer reviewed journal papers, symposium papers,
technical reports, and two books.
List of Figures

Fig. 2.1 a 1D column distribution, b 2D block-cyclic distribution . . . .. 14

Fig. 2.2 Variation in memory requirements of NEC at a particular node
with respect to number of segments . . . . . . . . . . . . . . . . . . . . .. 15
Fig. 3.1 Geometry of a wire segment . . . . . . . . . . . . . . . . . . . . . . . . . .. 18
Fig. 3.2 Surface patch a general geometry b rectangular patch
c triangular patch d quadrilateral patch . . . . . . . . . . . . . . . . . . .. 20
Fig. 4.1 Graphical user interface (GUI) of HyperMesh . . . . . . . . . . . . .. 24
Fig. 4.2 Workflow for the computation of RCS using NEC
in conjunction with HyperMesh . . . . . . . . . . . . . . . . . . . . . . . .. 25
Fig. 4.3 Procedure for computation of RCS using NEC
and HyperMesh . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 27
Fig. 5.1 Wire grid model of PEC dihedral in HyperMesh . . . . . . . . . . .. 30
Fig. 5.2 A portion of the NEC compatible segmentation data,
corresponding to PEC dihedral, extracted from Nastran file . . .. 31
Fig. 5.3 NEC input file for the computation of monostatic RCS
of PEC dihedral at 600 MHz . . . . . . . . . . . . . . . . . . . . . . . . . .. 31
Fig. 5.4 Monostatic RCS of PEC dihedral . . . . . . . . . . . . . . . . . . . . . . .. 32
Fig. 5.5 Wire grid model of NASA almond in HyperMesh . . . . . . . . . .. 33
Fig. 5.6 A portion of the NEC compatible segmentation data,
corresponding to NASA almond, extracted from Nastran file . .. 34
Fig. 5.7 NEC input file for the computation of monostatic RCS
of NASA almond (horizontal polarization) . . . . . . . . . . . . . . . .. 34
Fig. 5.8 NEC input file for the computation of monostatic RCS
of NASA almond (vertical polarization) . . . . . . . . . . . . . . . . . .. 35
Fig. 5.9 Monostatic RCS of metallic NASA almond. a Horizontal
polarization b vertical polarization . . . . . . . . . . . . . . . . . . . . . .. 36
Fig. 5.10 Wire grid model of metallic single ogive in HyperMesh . . . . .. 37
Fig. 5.11 A portion of the NEC compatible segmentation data,
corresponding to metallic single ogive, extracted
from Nastran file . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 38

xv
xvi List of Figures

Fig. 5.12 NEC input file for the computation of monostatic

RCS of metallic single ogive (horizontal polarization) . . . . . . . . 38
Fig. 5.13 NEC input file for the computation of monostatic
RCS of metallic single ogive (vertical polarization) . . . . . . . . . . 39
Fig. 5.14 Monostatic RCS of metallic single ogive. a Horizontal
polarization b vertical polarization . . . . . . . . . . . . . . . . . . . . . . . 40
Fig. 5.15 Cross-sectional view of the metallic rocket-shaped model . . . . . 41
Fig. 5.16 Wire grid model of metallic rocket shaped geometry
in HyperMesh . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
Fig. 5.17 A portion of the NEC compatible segmentation data,
corresponding to metallic rocket shaped model, extracted
from Nastran file . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Fig. 5.18 NEC input file for the computation of monostatic
RCS of rocket-shaped model . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
Fig. 5.19 Monostatic RCS of metallic rocket shaped model. . . . . . . . . . . . 43
Fig. 5.20 CAD model of metallic MAV . . . . . . . . . . . . . . . . . . . . . . . . . . 44
Fig. 5.21 Wire grid model of Metallic MAV in HyperMesh . . . . . . . . . . . 44
Fig. 5.22 A portion of the NEC compatible segmentation data,
corresponding to metallic MAV, extracted from Nastran file . . . 45
Fig. 5.23 NEC input file for the computation of monostatic
RCS of metallic MAV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
Fig. 5.24 Monostatic RCS of metallic MAV . . . . . . . . . . . . . . . . . . . . . . . 46
Fig. 5.25 CAD model of metallic UAV . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
Fig. 5.26 Wire grid model of metallic UAV in HyperMesh . . . . . . . . . . . 47
Fig. 5.27 A portion of the NEC compatible segmentation data,
corresponding to metallic UAV, extracted from Nastran file . . . 48
Fig. 5.28 NEC input file for the computation of bi-static
RCS of metallic UAV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Fig. 5.29 Bi-static RCS of metallic UAV model . . . . . . . . . . . . . . . . . . . . 49
Fig. 5.30 CAD model of metallic car . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Fig. 5.31 Wire grid model of metallic car in HyperMesh . . . . . . . . . . . . . 50
Fig. 5.32 A portion of the NEC compatible segmentation data,
corresponding to metallic car, extracted from Nastran file . . . . . 51
Fig. 5.33 NEC input file for the computation of monostatic
RCS of metallic car . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
Fig. 5.34 Monostatic RCS of metallic car model . . . . . . . . . . . . . . . . . . . . 52
Fig. A.1 Structure of a general MPI program . . . . . . . . . . . . . . . . . . . . . . 58
Fig. A.2 Format of a typical PBS script . . . . . . . . . . . . . . . . . . . . . . . . . . 59
List of Tables

Table 3.1 List of geometry commands . . . . . . . . . . . . . . . . . . . . . . . . . .. 21

Table 3.2 List of program control commands . . . . . . . . . . . . . . . . . . . . .. 22
Table 5.1 Simulation and meshing parameters to be used for the
computation of monostatic RCS of PEC dihedral . . . . . . . . . .. 30
Table 5.2 Simulation and meshing parameters to be used for the
computation of monostatic RCS of NASA almond . . . . . . . . .. 33
Table 5.3 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic single ogive . . . .. 37
Table 5.4 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic rocket-shaped
geometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 42
Table 5.5 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic MAV . . . . . . . . .. 45
Table 5.6 Simulation and meshing parameters to be used for the
computation of bi-static RCS of metallic UAV . . . . . . . . . . . .. 48
Table 5.7 Simulation and meshing parameters to be used for the
computation of monostatic RCS of metallic car . . . . . . . . . . . . 50
Table B.1 Options available for Int-1 (EX command) . . . . . . . . . . . . . . . . 63
Table B.2 Options available for Int-1 (RP command) . . . . . . . . . . . . . . . . 64
Table B.3 Options available for P (RP command) . . . . . . . . . . . . . . . . . . 64
Table B.4 Options available for N (RP command) . . . . . . . . . . . . . . . . . . 65
Table B.5 Options available for G (RP command) . . . . . . . . . . . . . . . . . . 65
Table B.6 Options available for A (RP command) . . . . . . . . . . . . . . . . . . 65

xvii
Chapter 1
Introduction

The burgeoning research and developmental activities happening in various domains

like antenna design, wireless communication, stealth technologies, electromagnetic
interference/compatibility (EMI/EMC), etc., have made systems much more complex
than ever before. Considering the huge expenses related to fabrication, it is not
practically possible to modify a device if its performance is not within the stipulated
limits. Therefore, techniques capable of precise characterization of complex systems
are the need of the hour. Closed-form analytical solutions are available for only
particular problems, which can hardly be directly applied in real-world systems.
This gap has led to extensive research on computational electromagnetics (CEM)
where numerical methods are used for solving electromagnetic problems. Among
all the numerical techniques, method of moments (MoM) is highly accurate and
efficient especially for the analysis of scattering from conducting structures. In this
regard, Numerical Electromagnetics Code (NEC) is a versatile MoM-based open-
source code (written in Fortran) used for the electromagnetic analysis of metallic
structures in the presence of sources or incident fields (Burke and Poggio 1981;
Burke et al. 2004). The development of NEC has been spearheaded by Lawrence
Livermore Laboratory, USA, with the support of the Naval Ocean Systems Center
and the Air Force Weapons Laboratory.
Since its development, NEC has continued to be one of the most widely used
electromagnetic simulation codes even in the presence of commercial MoM solvers.
This attributes to the availability of the well-documented powerful computational
platform, which can be used as a foundation for further developments in the area of
electromagnetic analysis. NEC does not use many of the simplifying assumptions
employed by other high-frequency asymptotic techniques and thus provides a highly
accurate and resourceful tool for electromagnetic studies.
The numerical solution of integral equations, for the induced currents on metallic
structures by sources or impinging electromagnetic waves, is the essence of NEC.
The excitations can be defined in terms of voltage sources or plane waves of
linear/elliptical polarization. Various electromagnetic entities like non-radiating

© The Author(s), under exclusive license to Springer Nature Singapore Pte Ltd. 2021 1
V. Joy et al., Fundamentals of RCS Prediction Methodology using Parallelized Numerical
Electromagnetics Code (NEC) and Finite Element Pre-processor,
SpringerBriefs in Computational Electromagnetics,
https://doi.org/10.1007/978-981-15-7164-0_1
2 1 Introduction

networks/transmission lines, perfect/imperfect conductors, perfectly or imperfectly

conducting ground planes, etc., can be easily included in NEC. The output of NEC can
be tailored for specific requirements like induced current and charge densities, near
magnetic or electric fields, impedance or admittance, radiated fields, etc. Further-
more, many of the commonly encountered performance parameters such as gain and
directivity are also available. Therefore, the program is well suited for both antenna
analysis and radar cross section (RCS) predictions (Burke and Poggio 1981).
Although numerical methods like MoM can be used in a wide range of radiation
and scattering problems, the computational complexity increases with increase in the
electrical size of the structure under consideration. The memory requirements have
been found to exhibit a quadratic growth with increase in the number of segments. A
parallelized version of NEC has been introduced (Rubinstein et al. 2003, Rubinstein
2004) to tackle this problem. The two most time-consuming steps in evaluating the
response of a structure, viz. (i) computation of the interaction matrix elements and
(ii) solution of the resultant matrix equation, have been efficiently parallelized in
Rubinstein et al. (2003).
In addition to the memory constraints, a serious issue in the computation of RCS
using NEC is the wire-grid-based meshing of geometries and the generation of NEC
input file with the segmentation data included in accordance with NEC guidelines.
This is an extremely difficult and error-prone process, especially for complex struc-
tures (Hubing et al. 1994; Toit and Davidson 1995; Ross et al. 1999), and further-
more, the code will often run in the presence of such errors and produce results that
are plausible but incorrect. Therefore, an appropriate pre-processor is essential for
the accurate prediction of RCS using NEC.
The aim of this book is to provide a strong foundation for the computation of
radar cross section (RCS) of aerospace structures using a parallelized version of
NEC in conjunction with a finite element pre-processor. This book has five chapters.
It begins with a thorough overview of the theoretical background of NEC including
all relevant aspects of formulation of integral equation and its numerical solution.
The strategy used for parallelization of the code is also included in Chap. 2. The
geometrical modelling of structures and the input/output formats are described in
Chap. 3. The revised methodology for the computation of RCS incorporating all the
required steps and details is discussed elaborately in Chap. 4 followed by case studies
and validations in Chap. 5.

References

Burke, G.J., and A.J. Poggio. 1981 January. Numerical Electromagnetics Code (NEC)—Method
of Moments, Lawrence Livermore National Laboratory, Livermore, California, Technical
Document, Rep. UCID-18834, 719.
Burke, G.J., E.K. Miller, and A.J. Poggio. 2004 June. “The Numerical Electromagnetics Code
(NEC)—A Brief history.” In Proceedings of IEEE Antennas and Propagation Society Symposium,
vol. 3, 2871–2874, Monterey, California.
References 3

Hubing, T.H., C.H.H. Lim, and J. Drewniak. 1994 March.“A Geometry Description Language for
3D Electromagnetic Analysis Codes.” In Proceedings of the 10th Annual Review of Progress in
Applied Computational Electromagnetics, Monterey, CA, 417–422.
Rubinstein, A., F. Rachidi, M. Rubinstein, and B. Reusser. 2003. A Parallel Implementation of NEC
for the Analysis of Large Structures. IEEE Transactions on Electromagnetic Compatibility 45
(2): 177–188.
Ross, J.E., L.L. Nagy, and J. Szostka. 1999. “CAD tools for vehicular antennas.” In Proceedings of
the 1999 IX National Symposium of Radio Science, Poznan, Poland, March 16–17.
Rubinstein, A. 2004 January. Simulation of Electrically Large Structures in EMC Studies: Applica-
tion to Automotive EMC, Ph.D. dissertation, École polytechnique fédérale de Lausanne (EPFL),
https://doi.org/10.5075, 156.
Toit, C.D., and D.B. Davidson. 1995. Wiregrid: A NEC2 Pre-processor. Journal of Applied
Computational Electromagnetics Society 10: 31–39.
Chapter 2
Theoretical Background
for the Computation of Radar
Cross-Section (RCS) Using NEC

The entire theoretical framework of NEC starting from the definition of integral
equations to the solution of matrix equation is briefly introduced in this chapter. The
basis functions used for representing the unknown current densities and the weight
functions used for testing are also described in detail. Further, the strategy used for
parallelization of NEC (Rubinstein et al. 2003) is also included here.

2.1 Formulation of Integral Equations (IE)

NEC employs both an electric-field integral equation (EF-IE) and a magnetic-field

integral equation (MF-IE) for evaluating the electromagnetic behaviour of conducting
bodies under consideration. EF-IEs are mainly used in the case of thin-wire geome-
tries with small conductor volumes. On the other hand, MF-IE, which cannot be
applied for the thin-wire case, is more suitable for voluminous structures with smooth
surfaces. EF-IE and MF-IE are coupled together for a structure containing both wires
and surfaces. The theoretical formulation presented in the following sections is based
on the work reported in (Burke and Poggio 1981).

2.1.1 Electric-Field Integral Equation (EF-IE)

The integral representation for the electric field due to a surface current distribution
−
→
J S on a perfectly conducting body (S) is given by,
     −
−
→ →
 r ) = (− jη/4π k)
E( JS r  · Ḡ¯ r, r  dA (2.1)
S

© The Author(s), under exclusive license to Springer Nature Singapore Pte Ltd. 2021 5
V. Joy et al., Fundamentals of RCS Prediction Methodology using Parallelized Numerical
Electromagnetics Code (NEC) and Finite Element Pre-processor,
SpringerBriefs in Computational Electromagnetics,
https://doi.org/10.1007/978-981-15-7164-0_2
6 2 Theoretical Background for the Computation …

where,
     
G r, r = k 2 I + ∇∇ g r, r
 
 
  e− j k r −r 

g r, r = 
 
 
r − r 


√ μ0
k = ω μ0 ε0 ; η = ; I = x̂ x̂ + ŷ ŷ + ẑ ẑ
ε0

r and r represent observation point and source point, respectively.

−
→
The integral equation for the current induced on S by an incident field E Inc is
obtained by imposing the following boundary condition for r ∈ S:
−
→ −
→
r ) × E Sca (
n̂( r ) + E Inc (
r) = 0 (2.2)

−
→
r ) is the unit normal vector to the surface at r and E Sca is the scattered field
where n̂(
−
→ −
→Sca
due to the induced current J s . Substituting (2.1) for E in (2.2) yields the final
EF-IE as,

−
→ →   2
−   
r ) × E Inc (
−n̂( r ) = (− jη/4π k) × n̂(
r) × J s (r ) · k I + ∇∇ g r, r dA
S
(2.3)

Since NEC employs thin cylindrical wires for representing the conducting surface
(S), the vector integral equation in (2.3) will be transformed to a scalar integral
equation. The assumptions applicable for the thin-wire approximation used in most
of the cases are given below:
(i) Currents in the transverse direction can be ignored with respect to the currents
along axial direction of the wire.
(ii) Circumferential variation in the currents along axial direction can be ignored.
(iii) The unknown current may be expressed as a filament along the wire axis.
(iv) The boundary condition pertaining to the electric field is required to be imposed
only in the axial direction.

2.1.2 Magnetic-Field Integral Equation (MF-IE)

The integral representation for the magnetic field due to a surface current distribution
−
→
J s on a perfectly conducting body is given by,
2.1 Formulation of Integral Equations (IE) 7
   
−
→ −→ 
H (
r ) = (1/4π ) J s (r ) · ∇  g r, r dA (2.4)
S

where the differentiation is with respect to the integration variable r .

−
→ −
→
If an external incident field H Inc induces a current J s on S, then the total magnetic
field within the perfectly conducting body vanishes. Hence, the boundary condition
for r immediately within the surface S can be expressed as,
−
→Sca −
→
H ( r ) + H Inc (
r) = 0 (2.5)

−
→ −
→
where H Inc is the incident field without the structure and H Sca is the scattered field
−
→
due to the induced current J s given by (2.4).
−
→
Substituting H Sca in (2.5) yields the final MF-IE as,
     
−
→ −
→ 
r0 ) × H Inc (
−n̂( r0 ) = n̂(
r0 ) × (1/4π ) lim J s r × ∇  g r, r dA (2.6)
r→
r0
S

r0 ) is the normal vector directed outwards at the surface point r0 .

where n̂(
Similar to EF-IE, the vector integral Eq. (2.6) can also be resolved into two scalar
expressions in the directions of the orthogonal surface vectors tˆ1 and tˆ2 where

r0 ) × tˆ2 (
tˆ1 ( r0 ) = n̂(
r0 ) (2.7)

2.1.3 EF-IE/MF-IE Combined Equation

The EF-IE/MF-IE combined equation is used for a structure consisting of both wires
and surfaces. In such cases, r in (2.5) is restricted to wires and the integral for
−
→Sca
E ( r ) extends over the complete structure. Similarly, r0 in (2.6) is restricted to
−
→
surfaces and the integrals for H Sca (
r ) extends over the entire structure. The integral
can be simplified by using the thin-wire approximation on wires.
The final coupled integral equations (Burke and Poggio 1981) are,
    2   
−
→  ∂  dS 
−ŝ · E Inc (
r ) = (− jη/4π k) I (s  )
k 2 ŝ · s  − g 
r , r
∂s ∂s 
L
    
→     2 
− ∂ 
− ( jη/4π k) Js r k ŝ − ∇  g r, r dA ;
∂s
S1

for r on surfaces of wire (2.8)

8 2 Theoretical Background for the Computation …

  
−
→  
tˆ2 ( r ) = (−1/4π )tˆ2 (
r ) · H Inc ( r) · I s  sˆ × ∇  g r, r dS 
L

− (1/2)tˆ1 (
r ) · JS (
r ) − (1/4π ) tˆ2 (
r)
S1
−
→  
· J S (r ) × ∇  g r, r dA ;
for r on surfaces excluding wires (2.9)
  
−
→
−tˆ1 ( r ) = (1/4π )tˆ1 (
r ) · H Inc ( r) · I (s  ) sˆ × ∇  g r, r dS 
L

− (1/2)tˆ2 (
r ) · JS (
r ) + (1/4π ) tˆ1 (
r)
S1
 
· JS (r ) × ∇ g r, r dA ;


for r on surfaces excluding wires (2.10)

where s is the parameter corresponding to distance along

 the wire axis at r and ŝ
is the unit vector tangent to the wire axis at r. and represents integration over
L S1
wires and integration over surfaces excluding wires, respectively.

2.2 Computation of Scattered Fields

The scattered fields due to induced currents can be computed using a far-field approx-
imation. This is applicable when the observation points are far apart from the induced
current distribution by distances much larger compared to the wavelength as well
as the dimensions of the current distribution. The far-zone scattered fields from a
conducting body consisting of a wire region with contour L having linear current
−
→ −
→
distribution I (s) and a surface region S with surface current density J s ( r ) are
given by,
⎧
−
→ ⎨   − → 


r0 ) = ( jkη/4π ) × (e− jkro /ro ) ×
E ( (k · I (s))k − I(s) e j k·r
⎩
L
⎫
   ⎬
 −→  −
→ 
+ (k · J s (
r ))k − J s (
r ) e j k·r dA (2.11)
⎭
S
2.2 Computation of Scattered Fields 9


where k = |rroo | , k = 2π
λ
and ro is the position of the observation point.
The contour integral can be computed in closed form over each wire segment
for the cosine, sine, and constant components of the basis functions and this in turn
becomes a summation over the wire segments. Since the surface current on the patch
is represented by a delta function at the centre of the patch, the second integral
transforms to a summation over the patch elements.
−
→
For an incident plane wave with electric field E Inc , the NEC program constants are
defined such that the radar cross-section (RCS) σ/λ2 is displayed under the column
gain in the output file.
− 
→Sca 2
σ E 
= 4π 
→ 2
(2.12)
λ2 −
 E Inc 

−
→
where E Sca is given by (2.11).

2.3 Numerical Solution of Integral Equations

The integral equations described in the previous sections are solved numerically in
NEC by using the method of moments (MoM). Method of moments starts with a
conventional linear-operator equation as,

L( f ) = u (2.13)

where f is an unknown response, u is a known excitation and L is an integral operator.

The formulation proceeds by expressing the unknown function f as a sum of
known basis functions f j (defined in the domain of the integral operator L) as,


N
f = βj f j (2.14)
j=1

The unknown coefficients (β j ) are obtained by taking the inner product of (2.13)
with a set of weighting functions wi , i.e.

wi , L( f ) = wi , u ; i = 1, 2, . . . , N (2.15)

The inner product in the present case can be defined as,


 p, q = r ) q(
( p( r ))dA
S
10 2 Theoretical Background for the Computation …

where integration is over the surface of the structure.

Substituting (2.14) in (2.15) gives


N
  
β j wi , L f j = wi , u ; i = 1, 2, . . . , N (2.16)
j=1

The above equation can be expressed in matrix form as,

[G] [B] = [U ] (2.17)

  
where G i j = wi , L f j ; B j = β j ; Ui = wi , u .
The solution to (2.17) is

[B] = [G]−1 [U ] (2.18)

2.3.1 Selection of Basis Functions and Weight Functions

The critical factors determining the accuracy and convergence of MoM are the weight
functions and the basis functions. NEC employs different basis and weight functions
unlike Galerkin’s method. The weighing functions (wi ) are chosen as a set of delta
functions in accordance with the collocation method, where the integral equation is
enforced at specific sample points on the geometry, i.e.

wi (
r ) = δ(
r − ri )

where {ri } are the set of points on the conducting surface. In NEC, conducting wires
are sub-divided into short linear segments where a sample point is defined at the
centre of each segment. On the other hand, surfaces are represented by a set of flat
patches with a sample point defined at the patch centre.
The linear current density on wires and surface current density on surfaces are
the unknown quantities in NEC which are to be represented by a set of known basis
functions. NEC uses sub-domain basis functions to simplify the calculation of inner
product integrals and also to ensure that the matrix G is well conditioned.
In the case of wires, the current on each segment is approximated as a sum of
three terms, viz., a sine, a cosine, and a constant. The total current on ith segment
can thus be expressed as,

Ii (s) = Ai sin k(s − si ) + Bi cos k(s − si ) + Ci , |s − si | < i

(2.19)
2
2.3 Numerical Solution of Integral Equations 11

where si is the value of s at the centre of segment i, i is the length of segment i,

and k is the free space wave number.
The additional equations required to evaluate all the unknown coefficients in (2.19)
for each and every segment will be provided by the local boundary conditions. The
continuity of current and charge is one such condition to be enforced at a junction
of two segments with uniform radius. In the case of a junction with more than
two segments of different radii, the continuity of current may be generalized to
Kirchhoff’s current law.
The final solution requires the computation of the electric field at each segment
due to these currents. NEC uses three approximations of the integral equation kernel
for the computation of fields. These include (i) thin-wire approximation for most
of the cases, (ii) extended thin-wire approximation for thick wires, and (iii) current
element approximation for the case where segments are separated far apart by a
considerable distance. In all the three methods used for the evaluation of field at
a particular segment due to the induced current on a different segment, the field
computation is performed on the surface of the observation segment.
NEC models surfaces using small flat patches and the unknown surface current
on each patch are expressed as a summation of a set of pulse functions except where
a wire connection is present. When a wire connects to a surface, a special local
boundary condition has to be enforced as explained in (Burke and Poggio 1981).
The unknown surface current density over M p patches can be written as,

−
→  Mp
 
J s (
r) = (J1i t 1i + J2i t 2i )vi (
r) (2.20)
i=1

   
where t 1i = t 1 (
ri ), t 2i = t 2 (
ri ), ri is the position of the centre of patch i, vi (
r) = 1
for r on patch i and 0 otherwise. Thus, every patch has two unknown quantities and
the equation is enforced for each vector component at the sample point in a given
patch.

2.3.2 Solution of Matrix Equation

The matrix equation obtained by appropriately substituting the basis functions and
weight functions in (2.16) has to be finally solved for the unknown currents. The
equation can be alternately expressed as,

[G] [I ] = [E] (2.21)

where [G] is the interaction matrix, [I ] is the column vector with unknown currents,
and [E] is the excitation vector. NEC employs Gauss elimination method for solving
(2.21). Here, the matrix G is first factorized into the product of upper and lower
12 2 Theoretical Background for the Computation …

triangular matrices, i.e.

[G] = [L] [U ] (2.22)

The matrix equation then becomes,

[L] [U ] [I ] = [E] (2.23)

Let,

[L] [F] = [E] (2.24)

where,

[U ] [I ] = [F] (2.25)

Equation (2.24) is first solved for F by forward substitution. Then (2.25) is solved
for I by backward substitution. The order of the matrix Eq. (2.21) for a structure
having N w wires segments and N p patches is given by,

N = N p + Nw

The final matrix Eq. (2.21) can be written in terms of sub-matrices as,
    
P Q IW EW
= (2.26)
R S IP HP

where I w is the column vector of segment basis function amplitudes, I p is the column
vector of patch-current amplitudes, E w is the LHS of (2.8) evaluated at segment
centres, and H p is the LHS of (2.9) and (2.10) evaluated at patch centres.
The sub-matrix element Pij denotes the electric field at the centre of segment i
due to the basis function centred on segment j. On the other hand, the sub-matrix
element S ij denotes the tangential magnetic-field component at patch u due to the
surface current pulse on patch v where
 
i −1
u = 1 + Int ; Int[ ]indicates truncation
2
 
j −1
v = 1 + Int
2

The elements in sub-matrix Q denote the electric fields due to surface current
pulses. On the other hand, the elements in R denote magnetic fields due to segment
basis functions.
2.4 Strategy for Parallelization 13

2.4 Strategy for Parallelization

The structure of NEC can be broadly divided into two parts, namely the input section
and the calculation section. The input section reads in the geometrical data pertaining
to the structure under consideration as well as the program control commands. The
calculation section on the other hand computes the elements of the interaction matrix
[G] which holds the geometrical and electromagnetic information about the model.
This section also deals with the solution of (2.21) by Gauss elimination method. This
final step is carried out in two routines. The first one performs the computationally
intensive task of LU decomposition of [G]. Later, it is used along with the excitation
vector [E] in another routine to obtain the solution. The computation of the constituent
elements of [G] and evaluation of the final solution to the matrix equation accounts
for more than 90% of the computation time (Burke and Poggio 1981). These parts
of the code need to be parallelized for the same reason.
The parallelized version of NEC used in this book is in accordance with (Rubin-
stein et al. 2003) where the subroutines responsible for LU decomposition (FACTR)
and solution of matrix equation (SOLVE) have been revised to work in a parallel envi-
ronment (Rubinstein 2004). This version distributes the interaction matrix among the
available processors and the matrix elements are locally and individually calculated
by their respective processors. The final system of linear equations is also solved
using a parallel version of the Gauss-Doolittle algorithm. In this way, the combined
processing power and memory of several processors working together as a cluster can
be efficiently utilized. The code can automatically detect the number of processors
and hence exploit the available memory and computing resources. The parallelized
code (Rubinstein et al. 2003) is designed to run on distributed memory parallel super-
computers and it can work on any system/cluster supporting conventional message
passing parallel environments such as Message Passing Interface (MPI) and Parallel
Virtual Machine (PVM). Since it is completely based on the original NEC code, the
input files need not be changed and the output format remains the same too.
The efficient distribution of the information stored in the matrices among the
available processors is very much critical for the successful parallelization of LU
decomposition and the parallel solution of the system of equations. For this reason,
the interaction matrix and the excitation vector have to be divided into smaller sub-
matrices which will be local to the respective processors. In (Rubinstein et al. 2003),
two-dimensional block-cyclic distribution has been used to uniformly distribute the
computation effort among all available processors and it is discussed in the following
paragraphs.
To illustrate two-dimensional block-cyclic distribution, consider a parallel
computer with four processors numbered as 1, 2, 3, and 4, as shown in Fig. 2.1. A
matrix is to be distributed among them. In a conventional way, as shown in Fig. 2.1a,
the matrix is distributed column wise to each processor. For ease of understanding,
the contents of processor 1 are shaded blue. Here, once all the elements that corre-
spond to processor 1 are assigned, it will become idle for the rest of the time. The
14 2 Theoretical Background for the Computation …

Fig. 2.1 a 1D column

distribution, b 2D
block-cyclic distribution

1 2 3 4 (a)

1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4
1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4
(b)
1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4
1 2 1 2 1 2 1 2
3 4 3 4 3 4 3 4

same situation repeats for the rest of the processors too. On the other hand, two-
dimensional block-cyclic distribution as shown in Fig. 2.1b guarantees that, on an
average, all the processors remain active from the start to the end of the algorithm.
A square processor grid is a desirable attribute for better efficiency as it maintains
equitable distribution of the computational load. But, in a practical scenario, the
number of processors varies from one machine to another. Therefore, the squarest
possible processor grid on the running platform is automatically calculated at the
beginning of the algorithm. ScaLAPACK library has been used in (Rubinstein et al.
2003) for implementing the block-cyclic decomposition in NEC. A brief description
of ScaLAPACK and MPI is given in Appendix A.
It is noted that the formation of the interaction matrix on a single processor
and then the distribution among other processors is not an efficient solution for
memory problem. The interaction matrix computation is therefore parallelized in
such a manner so as to produce an already distributed matrix as it is being computed.
Each processor on the grid will have adequate information required to determine
the elements of [G] once the input section is completed. An additional function
is employed to notify each processor whether the element (i, j) being calculated
corresponds to its local sub-matrix. In that case, another sub-routine will compute
the local values of (i, j) for the same element so as to be assigned to the corresponding
local matrix. Once the matrix generation is completed, each processor will hold its
local version of the interaction matrix and there would not be any further need for
2.4 Strategy for Parallelization 15

communication among the processors. Once the matrix equation is solved, data is
compiled by a block-cyclic composition routine and all the processors become aware
of the solution.
An approximation of the memory requirements in bytes at a local node for the
parallelized model of NEC is given by (Rubinstein et al. 2003),
 
16 × Nseg
2

Memory =
P
where P is the number of processors being used during the execution of the program
and Nseg is the number of segments. On the other hand, the original version of NEC
which is not parallelized requires (16 × Nseg
2
) bytes of memory. Figure 2.2 presents
the difference in memory requirements at a particular node between the original
version of NEC and its parallelized counterpart.
Considering a structure meshed with a particular number of segments, the vari-
ation in the run time of NEC with respect to the number of processors has been
studied in (Rubinstein et al. 2003), keeping all other parameters fixed. When the
number of processors was increased from one to four, it has been found that the
total run time reduced by almost 98%, thereby emphasizing the efficiency of paral-
lelization. However, the number of processors has to be judiciously chosen based on
the electrical dimensions of the structure under consideration, as employing large
number of processors may also increase the run time due to communication delay.

700

600 Original NEC

Memory requirements (MB)

Parallel NEC (P=2)

Parallel NEC (P=4)
500 Parallel NEC (P=6)
Parallel NEC (P=8)
400

300

200

100

0
1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000

Number of segments

Fig. 2.2 Variation in memory requirements of NEC at a particular node with respect to number of
segments
Another Random Scribd Document
with Unrelated Content
cette miséricordieuse servante des pauvres ? Il faut que vous la
lisiez.
Pauline rentra, l’âme rafraîchie d’une paix où elle sentit l’avant-
goût de plus parfaites délices ; de cette heure, sa conversion devint
une chose vraiment décidée.
Toutefois, s’il avait fallu pour la conduire au premier seuil de la
vie surnaturelle, une année de tourments, la mort de Julien, et les
inestimables supplications d’âmes qu’elle ignorait, il lui restait plus
d’un doute à résoudre et d’un dégoût à vaincre.
Jusqu’alors elle avait jugé le monde « bien pensant » d’après les
Rude et son oncle. Aussi le croyait-elle supérieur à celui dont elle
était. Mais Armance lui apprit sur des gens du voisinage, réputés
dévots, quelques anecdotes qui l’indignèrent.
Il y avait au bout de la rue, dans une maison décrépite qu’elle
louait presque en entier, une vieille fille riche et sordide, Mlle Crépin.
Pauline la voyait passer tous les matins, allant à la messe de sept
heures, ratatinée sous une pèlerine noire, coiffée d’une capote de
forme archaïque, et marchant en zigzag, comme si elle cherchait,
entre les fentes des pavés, des louis d’or perdus. Mlle Crépin, qui
passait pour millionnaire, accroissait ses revenus par des
spéculations habilement conseillées ; elle participait à la fureur
d’agiotage dont était possédée cette petite ville de rentiers oisifs ; et
l’on racontait qu’en un seul mois la hausse des cuivres lui avait valu
trente mille francs de gain. Elle se mêlait d’œuvres charitables, mais
appliquait au bien des pauvres les principes qu’elle suivait pour le
sien propre ; elle plaçait l’argent recueilli à leur intention, et, même si
elle les savait dans les plus affreuses nécessités, elle les rationnait
en aumônes, ne laissait fuir de « leur capital » que des bribes
dérisoires.
« Lorsqu’on n’a pas, on sait se priver », tel était un de ses
axiomes ; elle revenait souvent du marché avec trois navets dans
son cabas en se lamentant de ce que « la vie devenait impossible » ;
elle passait l’hiver sans feu, se chauffait les mains sur le couveau où
cuisait son potage. Quand elle n’était pas à l’église, elle comptait ses
coupons ou s’occupait de faire rentrer ses loyers ; et, une fois, en
grimpant à une soupente pour sommer d’en déguerpir le locataire
qui l’habitait, elle avait failli se rompre le cou.
Quoique Pauline connût son renom, l’idée audacieuse lui vint de
sonner à sa porte et de mendier une contribution au secours que les
Rouleau attendaient.
« Étant notre voisine, peut-être n’osera-t-elle pas me refuser ; et,
si elle me reçoit mal, j’embourserai l’affront pour les pauvres. »
Un dimanche donc, après vêpres, elle se présenta chez cette
personne de dévotion. Mlle Crépin entre-bâilla son huis, laissa voir
son menton aigu, son nez sec chaussé de lunettes, ses yeux
clignotants. Dès que Pauline se nomma, elle prit un air froid et
cérémonieux : elle l’introduisit dans la salle à manger où des images
pieuses ornaient les murs, et la pria de s’asseoir, sans insister. Mais,
aux premiers mots que Pauline prononça sur les Rouleau, la vieille
demoiselle se redressa dans sa petite taille ; ses lèvres, minces
comme un fil, se pincèrent.
— On m’a déjà parlé de cette famille, dit-elle d’un ton rogue ; ce
ne sont pas des gens intéressants ; l’homme et la femme ont vécu
ensemble non mariés.
— Je crois être sûre, répliqua Pauline, qu’ils sont mariés à l’heure
actuelle.
— Oui, je le sais, le ménage est régularisé — ce mot, dans une
bouche soi-disant chrétienne, sonnait déplaisamment — ; mais nous
n’avons pas cru devoir les mettre sur nos listes. Avec les ressources
modestes dont nous disposons, il nous faut imiter les vierges sages
qui gardèrent l’huile de leur lampe… Ne pouvez-vous pas vous
adresser au bureau de bienfaisance ? Ces messieurs vous
accueilleront.
Sa façon d’articuler : « ces messieurs », s’accompagna d’un coup
d’œil qui voulait dire : « Ils sont de votre bord ; nous, foin de vos
gens ! »
— Mademoiselle, conclut Pauline en faisant deux pas vers le
seuil, je suis venue à vous, parce qu’il s’agit d’une détresse urgente.
Vous m’opposez la prudence des vierges sages ; permettez-moi de
vous rappeler l’imprudence du bon Samaritain qui n’attendit pas,
pour verser de l’huile dans les plaies du moribond, de savoir sur
quelles listes il était.
Elle salua, et partit révoltée, moins de sa ladrerie que des motifs
dont elle croyait la couvrir. Comment ! cette hypocrite pharisienne
était considérée par les prêtres, admise à communier, et on
respectait en elle une des clefs de voûte de la paroisse ! Les
diatribes de son père contre la platitude cléricale lui remontaient aux
lèvres, et, comme Victorien, elle généralisait jusqu’à l’injustice :
« Est-ce l’Église vraie du Christ qui engraisse en ses pâturages
de pareilles brebis ? Cueille-t-on des raisins sur des épines et des
figues sur des ronces ? »
Le lundi, jour de marché, en sortant de bonne heure avec
Armance, elle traversa la place de la cathédrale où les merciers
tendaient leurs bannes.
— Mademoiselle, demanda la servante, veut-elle que j’entre à
l’église pour dire un bout de prière ?
— Oui ; je vous accompagne.
Une messe, à l’autel de la sainte Vierge, justement commençait ;
Pauline s’était initiée à l’ordonnance et aux phases du sacrifice ; elle
s’agenouilla dans une pensée de vénération émue, se remémorant
la parole : « Chaque fois que vous ferez cette chose, vous
annoncerez la mort du Sauveur » ; et, sans croire d’une foi pleine à
la Présence réelle, son esprit suivait attentivement la succession des
rites.
Mais le prêtre qui célébrait briffa l’Introït, la Collecte et l’Évangile
avec une vélocité qui la déconcerta.
« Dit-il sa messe pour lui seul ou pour les fidèles qui sont là ? »
Il mettait en ses génuflexions une nonchalance d’habitude
presque irrévérencieuse ; pendant le Canon, il traçait des signes de
croix sur la patène, éleva l’hostie, puis le calice, ôtait et remettait la
pale, tournait les feuillets du missel, se frappait la poitrine, communia
comme pressé d’en finir, et, en quinze minutes, la messe fut
expédiée. Pauline eut une déception, un froid lui tomba sur le cœur :
pour l’homme qui venait de réitérer la Cène, rien ne vivait donc sous
les mots et les gestes où il s’identifiait pourtant à Jésus-Christ ?
L’accoutumance émoussait-elle à ce point la ferveur ? Et alors, était-
ce la peine de pratiquer un culte dont les liturgies, à la longue, se
vidaient de toute émotion ?
Lorsqu’elle revit, le samedi d’après, l’abbé Charmoy, elle ne lui
dissimula rien de ses désenchantements. Il parut contrarié, mais en
prit occasion pour l’éclairer sur ses faiblesses qu’il pénétrait.
— Vous êtes trop impressionnable, la blâma-t-il tranquillement.
En principe, ce n’est point tout à fait un mal ; si vous sentiez peu,
vous vous seriez endurcie dans l’abstraction, et je ne connais guère
d’état plus triste, plus irrémissible. Il faut, néanmoins, apprendre à
gouverner vos sentiments, vous faire, comme disent les Provençaux,
une tête bien cerclée. Quand vous rencontrerez de mauvais
chrétiens, des prêtres négligents… ou même scandaleux, ne vous
pressez pas de conclure que l’Église, dont ils sont, est coupable de
leur indignité. D’abord, nous sommes plus tentés que les autres,
c’est incontestable. Interrogez votre jeune expérience ; vous aviez
plus de sécurité, de fausse sécurité, avant le jour où le premier appel
d’En Haut vous troubla. On ne mérite pas la grâce sans souffrance,
et il est si commode de s’engourdir, au lieu de s’évertuer ! Le démon
de la paresse glisse dans nos veines à notre insu ; nos plus belles
résolutions font souvent comme ces petits ermitages que sainte
Thérèse, enfant, bâtissait en posant les unes sur les autres des
pierres qui tombaient presque aussitôt. Cet abbé, dont la messe
vous afflige, il ne se doute pas, je crois, de son inconvenance. Il
oublie qu’on nous juge sévèrement, plus sévèrement que d’autres,
et avec raison parce que jamais la médiocrité, en nous, n’est licite.
Mais, vous, soyez plus humble ; chaque fois que l’esprit de critique
vous tourmente, même si vos griefs sont justes, appliquez-vous à
trouver dans l’œil du voisin une paille et, dans le vôtre, une poutre.
Pauline n’acceptait pas, sans résister, le langage de l’abbé
Charmoy ; elle en recevait pourtant la notion précieuse de
l’indulgence catholique et acquérait, à son contact, ce discernement
des « valeurs » que la raison indépendante oblitère, neuf fois sur dix,
par un manque d’équilibre. Le bon sens du prêtre n’était pas
simplement le sien ; sa pensée filtrait dix-neuf siècles de certitude
expérimentale et de tradition.
Une circonstance inopinée devait bientôt faire sentir à Pauline
combien la vie de l’Église s’incorporait à sa vie.
L’oncle Hippolyte, en octobre, se mit au lit ; ce vieillard, jusque-là
ferme comme un roc, déclina soudain de telle sorte que sa fin parut
prochaine. Il garda quelques semaines encore l’illusion de se
remettre, et, le 1er novembre, le ciel étant clair, il dit à sa nièce :
— Un beau jour de Toussaint… Comme en 1840. Si je pouvais
sortir demain… Il faut que j’aille à la banque toucher de l’argent… Tu
me prépareras mes bottines et mon manteau…
Son grand souci restait de manger le plus possible ; une heure
après son repas, il soutenait qu’il n’avait pas dîné et qu’on voulait le
laisser mourir de faim. Cependant, la paralysie gagnait ses organes ;
il dormait parfois des journées entières, avec une respiration si faible
qu’il ne semblait plus devoir se réveiller. Ou bien des hallucinations
obsédaient son cerveau dont les artères s’atrophiaient. Il parlait seul,
d’une voix sourde et absorbée, dans un délire sans fièvre. Il se
croyait invité à des ripailles et répétait durant des heures les
mouvements d’un homme qui mâche ou qui boit.
— Je ne crois pas qu’il aille bien loin, fit un soir M. Ardel, peiné
de perdre son oncle et davantage de voir la mort assise sur le toit de
sa maison.
Une autre anxiété préoccupait Pauline : « Mon oncle va-t-il mourir
sans sacrements ? » L’importance involontaire que prenait pour elle
un acte religieux l’avertit à quel point la foi devenait « l’os de ses os
et la moelle de ses moelles ». Malgré tout, elle n’osait en parler au
malade et justifiait sa timidité par des motifs contestables :
« Mon oncle a eu, en somme, une conduite probe. S’il s’est
racorni dans des enfantillages d’égoïste et d’avare, il a cru faire son
devoir en gagnant bien sa vie. Il n’a jamais eu beaucoup d’idées, et,
même, les gens qui en ont, pour lui, sont « des fléaux ». Dieu lui
pardonnera, parce qu’il aura beaucoup ignoré. Il s’en va plein de
jours, après une vieillesse somnolente et calme. A vrai dire, il aura
toujours été un dormant ; son entrée dans l’autre monde sera la
réelle naissance d’une âme qui n’a pas vécu. Dois-je l’éveiller avant
la lumière ? »
Mais, un jour que son délire avait cessé, comme Pauline, pour
l’égayer, parlait de la belle saison où il redescendrait au jardin :
— La belle saison, fit-il, je ne la verrai pas. Je suis au bout de
mon rouleau.
— Puissions-nous vous garder longtemps encore ! répondit-elle,
entraînée par une décision subite. Seulement, à votre âge, mon
oncle, des surprises sont possibles. Verriez-vous sans déplaisir un
prêtre ?
Il pâlit à cette question, comme s’il eût entendu son arrêt ; car, au
fond, il espérait vivre, et se disait mourant, dans l’espoir qu’elle le
rassurerait contre ses angoisses. Il croyait Pauline une parfaite
païenne, plus païenne que lui qui conservait pour les principes de
son enfance un respect latent ; si elle lui proposait un prêtre, c’était
donc qu’elle le savait bien fini.
— Quand le moment sera venu, répondit-il après un silence…
Mais amène-m’en un vieux… qui me comprenne…
Huit jours plus tard, en buvant une tasse de lait, il s’aperçut qu’il
ne pouvait avaler ; les muscles de sa gorge se paralysaient.
— Maintenant, bégaya-t-il d’une langue déjà embarrassée ; fais
ce que tu m’as dit.
L’abbé Charmoy vint le voir, le confessa, et lui apporta, le
lendemain, l’Extrême-Onction ; vers le coucher du soleil, il entra en
agonie. Pauline retrouva, près de son oncle moribond, certaines des
impressions qui l’avaient accablée au chevet de Julien : un râle
sifflait dans la poitrine du vieillard ; sa tête, renversée en arrière,
oscillait de gauche à droite et de droite à gauche, ses mains se
crispaient sous le drap. Elle lut tout haut les prières qu’elle avait
entendu lire à l’abbé Jacques ; et il lui semblait qu’elle les lisait pour
Julien. A son âme radieuse convenait cette anticipation des
triomphes célestes :
« Que la multitude splendide des Anges accoure au-devant de
toi ; que le sénat des Apôtres juges vienne aussi, et l’armée, vêtue
de blanc, des Martyrs ; que la troupe rutilante des Confesseurs
t’environne, portant des lis ; que le chœur des Vierges te reçoive ;
que les Patriarches te serrent dans leur embrassement au sein
d’une bienheureuse quiétude ; que le visage de Jésus-Christ
t’apparaisse doux et te fasse fête… »
Qu’importait-il vraiment de traverser le couloir sinistre de la mort,
si, à l’issue, doivent se déployer les portes éternelles ?
Pauline discernait quelle lourde erreur aveuglait Victorien,
lorsqu’il induisait des passagères défaites de l’Église la faillite de sa
mission. « L’Église visible, se disait-elle, s’appuie sur l’invisible
assemblée des Saints ; nous n’apercevons de ce concile immense
que les rangs infimes, et n’entendons de ses voix qu’un écho
assourdi… »
Peu de jours après l’oncle Hippolyte, le vieil archevêque décéda.
Son successeur, Mgr Chênedru, fit, trois mois plus tard, son entrée
solennelle ; Pauline et Edmée voulurent assister à cette cérémonie.
Les bourdons l’annoncèrent de leur formidable mugissement. Au
ventre de la tour leurs volées s’élargirent, et les chocs des deux
battants retombaient ensemble, comme titubant d’ivresse. Un
Hosanna pontifical se propageait avec les cercles ondulatoires de
leurs vibrations sur les collines et la campagne que le soleil de Mars
ranimait.
A l’intérieur, bruissait un vaste peuple ; on sentait dans l’attente
de la foule sourdre une allégresse ; et les pierres des arceaux
s’égayaient sous les oriflammes appendues. Précédé de ses
prêtres, Mgr Chênedru pénétra en sa cathédrale où son premier acte
fut de vénérer le reliquaire de la vraie Croix.
Edmée et Pauline se tenaient au milieu de la grande nef, sur le
passage du cortège ; elles reconnurent l’abbé Charmoy et l’abbé
Jacques ; celui-ci nageait dans une exultation ; il dirigea vers Pauline
un regard de victorieuse espérance. L’avènement du nouvel
archevêque signifiait pour lui la résurrection du diocèse ; et les plus
inertes paroisses tressailliraient sous la rafale d’enthousiasme qu’en
arrivant l’homme de Dieu suscitait déjà.
Pauline et son amie virent s’avancer près d’elles Mgr Chênedru ;
il éleva sa main et les bénit. Edmée s’inclina vivement pour baiser
l’émeraude de son anneau ; Pauline eut un léger recul ; n’étant pas
chrétienne, elle se trouvait indigne de cette faveur ; la recevoir eût
été un acquiescement de croyante, et elle n’allait pas encore jusque-
là.
Mais l’aspect de l’archevêque la conquit sur-le-champ. Bien qu’il
fût replet et d’une stature moyenne, il s’imposait par une puissance
d’autorité sans raideur et rayonnante ; il avait l’air d’être né pour tenir
la crosse et porter la mitre. Un visage opime, un œil qui étincelait,
des lèvres fines aux coins souriants d’où la parole semblait prête à
jaillir alors même qu’il se taisait, et surtout une bonté chaude qui
s’élançait au-devant des cœurs dans la persuasion de les atteindre,
tout faisait de lui une force en marche, douce et impérieuse ; quand
on le regardait, on ne pouvait plus douter que l’Église triomphante
existe.
Pendant qu’il lisait en chaire, d’une voix perçante, sa lettre
pastorale, Pauline sentit tomber ses dernières hésitations : d’une
telle bouche, comme de celle d’un apôtre qui aurait reçu du Christ
même sa doctrine, la vérité descendait.
Après le Salut, parmi la foule qui se pressait aux portes, elle
rencontra son oncle jubilant.
— Eh bien ! que vous semble de notre archevêque ?
— Je l’aime, répondit avant elle Edmée avec ferveur.
— Moi aussi, déclara Pauline, il m’a transportée.
— Je vous présenterai à lui, dit en les quittant l’abbé qu’un de
ses confrères entraîna vers la sortie.
La semaine suivante, dès que se fut écoulé le flot des réceptions,
il prit jour avec Mgr Chênedru et conduisit Pauline à l’archevêché.
Elle ne se vit pas introduite sans émotion dans le cabinet où
l’archevêque les attendait, debout derrière son bureau, contre sa
bibliothèque. Cette robe violette et l’idée d’omnipotence qu’elle
attachait à la dignité épiscopale lui imposaient une gêne :
« Que doit-il penser de moi, une infidèle ? »
Mais il embrassa, comme un père, l’abbé Ardel et fit à Pauline
elle-même un si affable accueil qu’elle reprit toute son aisance. Elle
exposa franchement l’état de son âme, l’indifférence orgueilleuse
d’où elle était partie, ses préventions d’ignorante contre les prêtres,
les doutes qui l’avaient retardée sur le chemin de la foi.
— J’étais encore incertaine, Monseigneur, continua-t-elle, quand
vous êtes venu parmi nous ; depuis que je vous ai vu et entendu, ma
décision est nette. J’ai médité tous les articles du Credo ; il n’en est
aucun auquel je ne puisse me soumettre, même l’enfer, bien qu’il n’y
soit pas nommé explicitement. Les peines éternelles me paraissaient
une chose monstrueuse, lorsque j’avais les plus fortes chances de
les mériter. Maintenant que j’espère ma rédemption, les rigueurs de
la justice divine ne me révoltent plus…
L’archevêque sourit de ce mot naïf et pénétrant.
— Réfléchissez et priez, ma chère enfant, approuva-t-il avec cet
accent du Béarn qui ajoutait une saveur à la bonhomie de son parler.
Je n’ai pas besoin de vous dire la grande joie que vous apporterez à
Notre-Seigneur et à moi, le jour où je pourrai vous donner le saint
baptême avant de vous confirmer.
— Dans combien de temps, demanda Pauline à son oncle
aussitôt qu’ils eurent pris congé de l’archevêque, pourrai-je être
baptisée ?
— A la Pentecôte, je pense. Avant peu, je reverrai Monseigneur
et l’abbé Charmoy.
Elle se sépara de lui et entra dans l’église qui allait être sa
paroisse, à Saint-Pierre-le-Rond. Fruste au dehors comme un vieux
sanctuaire de campagne, Saint-Pierre, sur une petite place toujours
déserte, est enclos entre les murs de logis silencieux. Pauline
affectionna, dès sa première visite, la nef étroite et longue, avec les
fenêtres du chœur découpées en trèfle. Un recueillement obscur
l’habitait où la pensée, mieux qu’ailleurs, pouvait « prendre son vol
sans bruit vers Dieu ». Elle se mit à genoux, dans le bas-côté, près
des fonts baptismaux, et là elle songea au mystère du Sacrement
qu’elle recevrait, à ce sceau du baptême qui, une fois imprimé sur un
front, ne peut plus s’en effacer.
Pourquoi cette efficacité surnaturelle de l’eau ? L’abbé Charmoy
lui avait un jour enseigné que toutes les eaux terrestres possédaient
une vertu de sainteté, depuis l’heure où Jésus sanctifia celle du
Jourdain en se courbant sous le baptême de Jean. L’eau purifie,
l’eau féconde, l’eau submerge ; elle se souvenait du verset d’un
psaume : Lavabis me et super nivem dealbabor, « tu me laveras et
je serai plus blanche que la neige » ; cette image lui plaisait en ce
qu’elle lui dévoilait le sens prophétique de sa course dans la neige
avec Julien. Maintenant il lui fallait s’ensevelir dans le Christ et
renaître avec lui, de même que les néophytes s’immergeaient trois
fois dans la piscine du baptistère primitif, de même que la novice du
Carmel, avant de prendre place au festin nuptial, s’anéantit sous le
drap mortuaire. Elle devait mourir à ses impiétés d’antan, à ses
vanités enfantines, à toutes les sensualités, et alors elle serait pure
comme la neige, comme l’eau d’une source où nul n’a jamais bu.
L’abbé Charmoy n’avait pas hésité à le lui dire : « Vous tomberiez
morte à l’instant de votre baptême, le Paradis vous recevrait aussi
sainte qu’un enfant qui n’a pas encore péché. »
Mais la vertu de l’eau n’opérait qu’unie à la parole, au nom du
Père, du Fils et du Saint-Esprit. Les Trois Personnes devant qui les
Anges trouvent à balbutier un seul mot : Sanctus, Sanctus, Sanctus,
viendraient donc en elle illuminer son âme et la diviniser !
« Qu’ai-je fait, ô Seigneur, pour mériter vos dons ? Je vous ai
dédaigné, nié, crucifié. Et, à présent, est-ce que je vous aime ? Est-
ce que je désire vous faire aimer ? »
Elle s’abîmait dans la repentance de ses rébellions, de son
aridité, de sa tiédeur, quand, tout d’un coup, une voix mauvaise, au
fond d’elle, articula :
 — Si, pourtant, tu t’abusais, si tu faisais la demande et la
réponse, si, là-bas, dans ce tabernacle, il n’y avait rien…
Mais elle repoussa l’idée affreuse, la piétina ; elle eût préféré
mille fois ne plus vivre que de ne plus croire ; et elle sortit, comblée
de la certitude que ses doutes étaient tués à tout jamais.
L’archevêque décida que le baptême, la confirmation et la
communion seraient donnés par lui à Pauline dans son oratoire, la
veille de la Pentecôte. Elle avait choisi pour parrain M. Rude, et
Mme Rude, en dépit de sa répugnance à se mêler d’aucune fête,
même liturgique, accepta d’être sa marraine. N’était-ce pas Julien
lui-même qui l’exigeait ? Et ce baptême serait l’accomplissement des
fiançailles mystiques de Julien avec sa bien-aimée.
Les Rude estimaient prudent de n’en point parler d’avance autour
d’eux. La nouvelle s’ébruita on ne sut comment, et les gens hostiles
à Victorien clabaudèrent par toute la ville.
— Hein ! répétait Galibert, cet Ardel qui se vantait de n’avoir pas
fait baptiser sa fille, ils l’ont retourné comme un gant !
Le plus venimeux furent Mlle Total, qu’il avait cessé de voir, et
Flug, avec qui il s’était brouillé, après avoir qualifié de « stupides »
ses paradoxes délirants. Un journal publia un entrefilet plaisantin sur
« les dragées du baptême », et ajouta que « seul, le poupon
manquerait ».
Informé de ces ragots par un obligeant collègue, Victorien
haussa les épaules, et répliqua très fort :
— Oui, ma fille va être baptisée, et c’est moi qui l’ai voulu. Leur
clique me dégoûte si bien qu’ils me donnent envie d’aller à la messe,
pour les faire enrager !
Depuis sa conversion, Pauline témoignait à son père une
tendresse de plus en plus prévenante ; elle évitait à son amour-
propre endolori les moindres blessures. Mais elle ne lui avait dit mot
de l’événement qui se préparait pour elle ; ce mutisme le peinait ; il
finit par se décider à le rompre lui-même :
— C’est pour bientôt la cérémonie ?
 — Pour samedi, répondit-elle avec un battement de cils et en
rougissant.
Et elle insinua d’une façon câline :
— Tu n’y viendras pas ?
— Non, je serais un trouble-fête, une fausse note. L’archevêque
et moi, vis-à-vis l’un de l’autre, nous serions mal à notre aise…
Elle n’insista point, ayant peur de moins bien prier, si la présence
de son père incrédule pesait sur sa ferveur.
Le matin du grand jour, elle se réveilla, comme une mariée qui va
mettre sa robe de noces, dans une attente extraordinaire. Elle avait
jeûné la veille, et son esprit se mouvait, presque dégagé de son
corps, avec une alacrité lumineuse. Dès six heures, elle sortit,
devant rejoindre, en la chapelle de son couvent, l’abbé Charmoy qui
la confesserait.
Il avait plu avant l’aurore ; sur le mail, où personne ne passait, un
vent d’est léger, le « matinal », comme disent les paysans de
Bourgogne, agitait dans les feuilles mouillées des platanes la
lumière aussi fraîche que la rosée. La nappe verte et claire des
frondaisons d’un acacia remuait dans le vivier du ciel ; les rossignols
se répondaient à travers les jardins ; une buée fumait sur des
massifs de fleurs ; la tour de la cathédrale était rose au soleil
montant.
« Tout à l’heure, se disait Pauline, je serai joyeuse comme ces
atomes de rayons qui dansent et qui scintillent. »
Mais la perspective de sa confession couvrait encore d’une
ombre le bonheur dont elle palpitait. Quoique ses entretiens avec
l’abbé Charmoy eussent, d’avance, allégé, pour elle, l’humiliation
des aveux, elle entra, presque tremblante, à l’intérieur du
confessionnal. L’exiguité noire et nue du recoin où ses yeux ne
distinguaient qu’une image de Jésus en croix et la grille fermée
d’une planche l’inquiétait comme un accusé qui attend, dans une
cellule austère, le moment de comparaître devant un juge infaillible.
Elle entendit l’abbé Charmoy enfiler son surplis, mettre son étole et
s’asseoir : était-ce le même prêtre dont elle connaissait le visage
bénin ? Mais, dès qu’il eut ouvert la grille et parlé, elle respira. A
chacune des fautes qu’elle énumérait scrupuleusement, il
prononçait, pour l’encourager, un : Bien, paisible. Son exhortation fut
une parole, moins de reproche que d’espoir grave. Pauline s’étonna
de la pénitence facile qu’il lui infligea ; trois psaumes à lire pour dix-
neuf ans d’infidélité ! Une critique qu’elle fit taire s’ébaucha en elle, à
l’idée d’une indulgence si exorbitante !
Elle lut aussitôt les trois psaumes ; car c’était ceux précisément
qu’on récite dans la liturgie du baptême, et, de tout son cœur, elle
s’appropria ces versets :
« Seigneur, notre Dieu, comme votre nom est admirable sur la
terre ! Votre magnificence est élevée au-dessus des cieux… Qu’est
l’homme, pour que vous vous en souveniez, et le fils de l’homme
pour que vous le visitiez ? Vous l’avez établi un peu au-dessous des
anges, vous l’avez couronné d’honneur et de gloire, vous l’avez
constitué sur les œuvres de vos mains…
« Comme le cerf désire les sources des eaux, ainsi le désir de
mon âme va vers vous, ô Dieu !… Quand viendrai-je et quand
paraîtrai-je devant la face de Dieu ? Mes larmes ont été, jour et nuit,
mon pain, tandis qu’on me disait : Où est ton Dieu ?… »
Elle partit en se chantant comme une mélodie les mots
extatiques : Quare tristis es, anima mea ?… Pourquoi étais-tu triste,
ô mon âme, et pourquoi me troublais-tu ? Espère en Dieu, puisque
tu le confesseras.
Les ailes de sa joie la portaient ; elle aurait couru sur des
charbons ardents avec l’illusion de marcher sur des roses. L’espace
se faisait bleu comme le vitrail du Paradis, dans la cathédrale ; elle
pensait, les yeux dirigés vers le soleil, à la vision de la Sibylle qui
aperçut, autour de l’astre, un cercle d’or, et au milieu du cercle une
Vierge merveilleuse, portant contre sa poitrine un enfant.
Armance et Antoinette, qu’elle avait invitées toutes deux à son
baptême, l’attendaient devant la porte de l’archevêché. Bientôt, le
parrain et la marraine arrivèrent avec Edmée et Marthe ; le grand
voile noir de Mme Rude et d’Edmée semblait cacher derrière elles le
fantôme de Julien. L’abbé Jacques et l’abbé Charmoy les suivirent
de près ; le secrétaire de l’archevêque, un jeune prêtre suave et
modeste, les fit tous monter dans l’oratoire, une chambre peu vaste
transformée en chapelle, et qui faisait songer à ces réduits où les
prêtres réfractaires, sous la Terreur, célébraient la messe. Mgr
Chênedru, en pluvial violet, entra presque aussitôt ; il s’agenouilla et
se recueillit ; on sentait dans son oraison muette qu’il soulevait vers
le Très-Haut les misères et l’imploration de tout un peuple ; en
baptisant Pauline il restituait au Christ une France qui ne peut cesser
d’être à Lui.
Il se tourna vers l’assistance, et s’adressant à la néophyte,
montra le prodige des largesses que Dieu, en un seul moment, allait
faire pleuvoir sur elle à pleines mains, la veille du jour où les langues
de feu étaient descendues, où les sept dons du Paraclet emplirent
les apôtres. Il évoqua les voies singulières par où elle avait été
conduite ; des allusions chaleureuses et pleines de tact à l’abbé
Ardel, à l’influence tacite des Rude, à Julien, à l’abbé Charmoy,
touchèrent d’un trait si juste le cœur de chacun que Mme Rude et
Edmée rabattirent leur voile devant leur figure, afin de pleurer
librement.
Mais, ajouta Mgr Chênedru, ce n’était point pour elle seule qu’elle
devait être chrétienne ; il fallait que sa naissance à la grâce fût un
signe et un exemple, et qu’autour d’elle la lampe ardente de sa piété
resplendît…
Ensuite, le baptême commença. L’archevêque, s’étant assis,
énonça, selon les formules rituelles, les mêmes questions que les
évêques des premiers siècles posaient, dans les catacombes, aux
jeunes chrétiennes de Rome.
Pauline y répondait en latin, et, chaque fois qu’elle réitérait le
simple mot : Credo, la conviction de sa foi s’implantait plus avant
dans son être, par cela seul qu’elle l’affirmait.
Puis, il se leva, l’exorcisa en soufflant sur elle ; et elle s’humilia
sans effort sous l’idée que sa personne avait pu être un temple de
l’Esprit immonde. Son âme, à cette heure, était souple, fondue
d’amour, telle que l’or liquide et rouge, quand on le verse dans le
creuset.
Il lui fit avec le pouce une croix sur le front et dit en même
temps :
— Signe ton front, pour que tu reçoives la Croix du Seigneur.
Et il continua :
— Signe tes oreilles, pour que tu entendes les divins préceptes.
Tes yeux, pour que tu voies la clarté de Dieu. Ton nez, pour que tu
sentes l’odeur de suavité du Christ. Ta bouche, pour que tu dises les
paroles de vie. Ta poitrine, pour que tu croies en Dieu. Tes épaules,
pour que tu prennes sur toi le joug de sa servitude…
Le Christ prenait possession de sa servante, l’investissait tout
entière, la voulant sienne « dans les siècles des siècles ».
L’archevêque exorcisa et bénit le sel qu’il mit sur la langue de
« l’Élue », afin que ce principe de force et de sagesse pénétrât dans
sa chair et y demeurât éternellement. Le parrain et la marraine
marquèrent, à leur tour, avec le pouce, le front de Pauline d’un signe
de croix. De la main du père et de la mère qui, par Julien, avaient
mis en elle les premiers rudiments de sa croyance, ce geste, trois
fois recommencé, équivalait à une attestation de leur paternité
acquise dans la douleur ; et ce fut, pour eux tous, une des minutes
les plus solennelles de la cérémonie.
Pauline suivait sur son Rituale romanum le sens intime des
oraisons, en apparence impersonnelles, mais exactement faites à
son intention. Dans un des exorcismes le célébrant disait :
— Tentateur maudit, ne viole jamais ce signe de la Croix sainte
que nous mettons sur son front… Va-t’en, tremblant et gémissant.
C’est Jésus-Christ qui te le commande, lui qui marcha sur la mer, et
tendit sa droite à Pierre qui sombrait.
C’était là une des images où elle se reconnaissait le plus
familièrement : l’élan de Pierre marchant sur les vagues à la
rencontre du Maître qu’il avait d’abord pris pour un fantôme, son cri
d’angoisse : « Seigneur, sauve-moi ! » et la main toute-puissante
tendue à sa faiblesse : « Homme de peu de foi, pourquoi as-tu
douté ? »
De même, plus loin, l’invocation au Dieu qui a ouvert les yeux de
l’aveugle-né lui remémora le mot de Julien, si vrai dans sa sévérité !
L’archevêque lui imposa sa main sur la tête ; avec lui et les
assistants elle prononça le Credo et le Pater ; il trempa son pouce
dans l’huile sainte, accomplit des onctions sur la poitrine et entre les
épaules de celle que le baptême allait sanctifier. Car l’instant était
venu pour Pauline de recevoir l’eau de la vie éternelle ; sa personne
était soustraite au Prince de ce monde ; elle pouvait devenir le
tabernacle de l’Esprit-Saint.
Une fois encore l’archevêque, en latin, lui demanda :
— Crois-tu au Dieu omnipotent, créateur du ciel et de la terre ?
— J’y crois.
— Crois-tu en Jésus-Christ, son fils unique, notre Seigneur, qui
est né et qui a souffert ?
— J’y crois.
— Crois-tu en l’Esprit-Saint, en la sainte Église catholique, en la
rémission des péchés ?…
— J’y crois.
— Veux-tu être baptisée ?
— Je le veux.
Alors elle s’inclina, il lui versa trois fois l’eau sainte sur sa tête ;
puis il lui mit un cierge entre les doigts, comme à une Vierge prête à
suivre le cortège de l’Époux.
Tout à l’heure, il demandait au Christ pour elle, en l’une des
oraisons, « de ne pas la laisser avoir faim longtemps, jusqu’à ce
qu’elle fût rassasiée de la nourriture céleste ». Cette nourriture, elle
l’attendait avidement. Lorsque l’archevêque l’eut confirmée, il ôta sa
mitre, revêtit une chasuble et dit la Messe, que lui servirent l’abbé
Charmoy et l’abbé Jacques.
 Pauline ne venait pas en vain de recevoir l’Esprit de sagesse et
d’intelligence. Tandis que la Messe se développait, elle entrait — ce
qu’elle n’aurait su faire auparavant — dans la sublimité du mystère
célébré devant elle et avec elle, puisque les chrétiens présents
officiaient, selon leur part de ferveur, en même temps que le prêtre
et l’invisible Officiant qui s’immolait.
Toute signée de la croix, elle la retrouvait multipliée sur la pierre
de l’autel, sur la chasuble, sur les instruments du sacrifice, dans les
gestes du célébrant. Mgr Chênedru articulait d’un ton haut les
prières du rite ; à la Consécration, il baissa la voix, mais proféra
lentes et distinctes les syllabes miraculeuses. Pauline sentit
réellement s’opérer la divine Présence, elle se vit couverte du sang
brûlant de la Victime ; elle aurait été confondue de tristesse en
pensant qu’elle-même avait ouvert ces veines et transpercé cette
chair, si l’attente de la communion ne l’eût saturée d’un bonheur
qu’ensuite elle s’étonna d’avoir pu porter. Ah ! comment des
hommes pouvaient-ils croire vivre, en ignorant de telles extases !
Elle ne se laissa point aller pourtant à une adoration passive. Elle
pria pour son malheureux père :
« S’il ne se convertit, ô mon Dieu, disait-elle, c’est que je ne
saurai pas vous aimer… »
Elle pria pour sa mère défunte, pour toutes les âmes perdues,
pour la déplorable paroisse de son oncle, pour le diocèse dénué de
prêtres, pour la France à ressusciter. Elle pria pour les pauvres sans
consolateur, pour les morts dont nul ne se souvient, pour les juifs et
les hérétiques, pour les immenses peuples qui seront idolâtres
jusqu’à la fin des temps…
Puis elle revint à ceux qu’elle aimait, aux Rouleau, aux deux
servantes, aux bons Rude, à l’abbé Charmoy, au saint archevêque,
à Julien qui lui méritait sa félicité. Elle s’unissait à lui dans le Christ,
comme jamais un amour terrestre ne les aurait unis. Dans la salle du
festin où l’Époux les conviait tous deux, elle entrait avec sa robe
blanche, immaculée, sa robe baptismale qu’elle ne quitterait plus.
 IX

Pauline était, depuis deux ans, chrétienne. La joie de son

baptême continuait, approfondie par l’intimité des Sacrements et la
richesse de méditations ardentes. Des peines cependant
l’obscurcissaient par intervalles : d’abord, elle avait honte d’elle-
même, quand elle évaluait son peu de charité, sa médiocre ferveur
de pénitence. Une chose l’humiliait surtout : elle ne pouvait prier,
même un temps court, sans distraction. Quelquefois la claire vue de
ses insuffisances la décourageait ; elle eût volontiers renoncé à
l’effort, se croyant vouée à trop d’imperfection. Puis elle rebondissait,
opposait ce qu’elle était à ce qu’elle avait cessé d’être, et s’exaltait
d’une gratitude inexprimable, lorsqu’elle mesurait son changement.
Mais il lui pesait de ressonger à ses années vaines : comme tout
cela était loin maintenant ! Une seule amertume les prolongeait,
l’incrédulité persistante de M. Ardel.
Il avait néanmoins changé, lui aussi. Le matin du baptême, au
retour de sa fille, dans ce beau visage une transfiguration l’avait
frappé ; un autre sang paraissait couler en ses joues, et la
transparence heureuse de ses prunelles renvoyait une lumière
séraphique. Il ne songea plus à nier que les vieux rites de l’Église
continssent encore une efficacité vitale. Mais il s’attendait à voir
Pauline, enflée par l’orgueil de sa conversion, s’éloigner de lui ; au
rebours, elle resta simple, affectueuse, soumise à ses désirs. Elle
rappelait une des figures de l’incomparable tapisserie du Trésor,
l’Esther couronnée par Assuérus, modeste dans sa gloire, comme si
elle devait en être toujours indigne. L’arome de paix qui sortait d’elle
agit peu à peu sur l’aigreur de Victorien ; il supportait plus
légèrement les déboires de sa carrière ; ses méfiances
s’atténuaient ; la sympathie plus équitable qu’il accordait aux
croyances de sa fille modifiait l’ensemble de son attitude critique.
Seulement, endurci à saisir les faits sous l’angle sec de l’intelligence,
il ne concluait pas le moins du monde qu’elle eût raison de croire, ni
qu’il dût la suivre.
— A mon âge, lui redit-il une fois, comme il l’avait dit à son frère,
on ne change guère son pli.
Il voyait souvent les Rude dans une amitié de plus en plus
étroite, et le peintre venait de lui annoncer un cadeau dont il était
charmé : le portrait de Pauline. Toutes les semaines, elle passait
donc, chez M. Rude, une après-midi. Elle s’asseyait, prenait un livre
captivant, et c’était en cette attitude de liseuse que l’artiste la fixait.
Un jeudi d’avril, M. Rude lui dit, d’assez mauvaise humeur :
— Nous n’aurons pas une longue séance aujourd’hui ; une visite
nous dérangera, Gabriel Authelin avec sa mère.
Gabriel Authelin remplaçait, dans la chaire de philosophie, Flug
que les suites d’une extravagance avaient contraint de s’en aller. A
propos d’une dissertation sur le mot de Montaigne : « Tout ciel m’est
un », Flug avait exposé que, pour le philosophe, la notion de patrie
demeurait inexistante : « Il me serait égal, avait-il déclaré, d’être
Allemand aussi bien que Français. » Là-dessus, deux de ses élèves,
se levant, avaient quitté la classe ; les autres, sauf un seul, s’étaient
empressés d’en faire autant ; des familles s’étaient plaintes, ses
chefs l’avaient blâmé, de sorte qu’il jugea prudent de porter ailleurs
sa métaphysique. Venant après lui, Authelin semblait justifier la
théorie platonicienne sur le rythme des contraires ; dogmatiste et
catholique, il était le neveu de cet abbé Authelin qui assista Mme
Rovère [1] dans sa maladie et sa prodigieuse guérison. Il avait connu
Daniel Rovère, le doux martyr, mort à Tarragone où les Chartreux
l’avaient recueilli. Sa philosophie, imbibée du mysticisme de Blanc
de Saint-Bonnet et d’Hello, y ajoutait un sens de la vie concrète,
d’autant plus surprenant qu’il était aveugle.
[1] V. l’Immolé.

C’était à l’âge de quatre ans, quand il vivait avec sa mère déjà

veuve et ses trois frères, à la campagne, près de Lyon, sur les
hauteurs du mont Cindre, qu’au moment d’un orage un coup de
foudre l’avait terrassé et avait brûlé ses yeux. Mais cette privation de
la vue stimula ses facultés natives ; ses autres sens s’étaient
emparés du monde extérieur avec une finesse suraiguë. A quinze
ans, il parlait sept langues ; sa mémoire, comme sa dialectique, se
faisait un jeu des connaissances les plus compliquées. Et il souffrait
peu de n’y plus voir ; car il se conduisait seul au dehors, distinguait si
son chemin était à droite ou à gauche, s’il longeait une place ou une
rue. Il voyait par les oreilles et le toucher ; la canne dont il s’aidait lui
communiquait sur les objets voisins des données précises. Ses
doigts lisaient aussi aisément que l’eussent fait ses yeux ; et, quand
il aimait un livre, sa mère patiente le copiait à son usage d’aveugle.
Sa vie méditative s’accroissait de tout ce que ses regards pouvaient
perdre ; il disait que sa « chambre obscure » ressemblait à certaines
chapelles de la cathédrale Saint-Jean où les ténèbres, en plein midi,
restent opaques, pour que l’on y puisse mieux faire oraison.
L’abbé Ardel, qui avait rencontré à Lyon Mme Authelin, lui inspira
le désir de connaître les Rude ; c’est pourquoi, ce jeudi, elle devait
leur conduire son fils.
Pauline et Edmée attendaient curieusement cette visite. Gabriel
entra, suivant sa mère, une femme de noble mine, plus grande que
lui, lente et mesurée dans sa démarche, par l’habitude qu’elle avait
de se mettre au pas de l’aveugle. Il tâtait, du bout de sa canne,
d’une façon discrète, le plancher. Il atteignit un fauteuil et s’assit
sans embarras. On se fût à peine douté, en l’apercevant, qu’il n’y
voyait rien. Il tenait ses paupières baissées, à la façon d’un
somnambule ; mais son front bombé, poli comme un marbre, ne
laissait point voir ce plissement douloureux, si habituel chez les
aveugles. Ses cheveux étaient longs, bruns comme sa barbe ; il
Welcome to our website – the ideal destination for book lovers and
knowledge seekers. With a mission to inspire endlessly, we offer a
vast collection of books, ranging from classic literary works to
specialized publications, self-development books, and children's
literature. Each book is a new journey of discovery, expanding
knowledge and enriching the soul of the reade

Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.

Let us accompany you on the journey of exploring knowledge and

personal growth!

textbookfull.com