Numerical Analysis

Archive for May, 2010

Call for Papers: IAENG International Workshop on Scientific Computing and Computational Statistics (in IMECS 2006)

CFP From:
IAENG: International Association of Engineers (
Engineering Letters (

Call for Papers: The 2006 IAENG International Workshop on Scientific
Computing and Computational Statistics

(Part of The International MultiConference of Engineers and Computer
Scientists IMECS 2006)

20-22 June, 2006, Hong Kong

The IWSCCS’06 has the focus on the frontier topics in the theoretical
and applied scientific computing and computational statistics subjects
and their applications. It serves as good platforms for the scientific
computing and computational statistics community members to meet with
each other and to exchange ideas.

The IMECS 2006 multiconference has the focus on the frontier topics in
the theoretical and applied engineering and computer science subjects.
It consists of 14 workshops (see the details at IMECS website: The multiconference serves as good platforms
for the engineering community members of different disciplines to meet
with each other and to exchange ideas. The current conference committee
of the IMECS 2006 includes over 140 workshop co-chairs and committee
members of mainly research center heads, department heads, professors,
research scientists from over 20 countries, while a few of the
committee members are also experienced software development directors
and engineers.

All submitted papers will be under peer review and accepted papers will
be published in the conference proceeding (ISBN: 988-98671-3-3). The
abstracts will be indexed and available at major academic databases.
The Technology Research Databases (TRD) of CSA (Cambridge Scientific
Abstracts), DBLP and Computer Science Bibliographies have promised to
index the print proceeding in advance of its publication. And after the
publication of the proceeding, print copies will also be sent to
databases like IEE INSPEC, Engineering Index (EI) and ISI Thomson
Scientific for indexing. The accepted papers will also be considered
for publication in the special issues of the journal Engineering
Letters. Some participants may also be invited to submit extended
version of their conference papers for considering as book chapters
(soon after the conference).

IWSCCS Workshop Co-Chairs:

Chin-Chen Chang
IEEE Fellow, IEE Fellow
Chair Professor in Department of Information Engineering and Computer
Science, Feng Chia University, Taiwan

Chiang-Ho Cheng (co-chair)
Assistant Professor, Department of Mechanical and Automation
Da-Yeh University, Taiwan

Dr. Alexandar Djordjevich
Programme Leader for the BEng in Mechatronic Engineering Programme
City University of Hong Kong, Hong Kong

G. Ganesan
Associate Professor of Mathematics
Jayaprakash Narayan College of Engineering, India

Wan-Kai Pang (co-chair)
Assistant Professor, Department of Applied Mathematics
The Hong Kong Polytechnic University, Hong Kong

Josep R. Herrero (co-chair)
Assistant professor, Department of Computer Architecture
Polytechnic University of Catalonia, Spain

Chih-Jer Lin (co-chair)
Associate Professor of Mechanical and Automation Engineering
Da-Yeh University, Taiwan

Tadeusz Sawik
Professor and Head
Computer Integrated Manufacturing,
AGH University of Science & Technology, Poland

Lev V. Utkin, PhD, DSc, (co-chair)
Professor, Dept. of Computer Science
St.Petersburg Forest Technical Academy, Russia

Jiunn-Lin Wu
Assistant Professor, Dept. of Computer Science,
National Chung Hsing University, Taiwan

An-Shik Yang (co-chair)
Associate Professor, Department of Mechanical and Automation
Da-Yeh University, Taiwan

The topics of the workshop include, but not limited to, the following:

Numerical analysis
Finite difference methods
Finite element methods
High order difference approximations
Methods for integration on a uniform and non-uniform mesh
Monte Carlo methods
Numerical Linear Algebra
Factorization methods
Transformation methods and applications
Dynamical systems
Bayesian computing
And applications of scientific computing
Computer intensive inferential methods
Data exploration and data mining
Neural networks
Numerical methods
Statistical databases
Statistical systems
And statistical software


Prospective authors are invited to submit their draft paper in abstract
format (one page) or in full paper format to im… by 12
March, 2006. The submitted file can be in MS Word format, PS format, or
PDF formats.

The first page of the draft paper should include:

· Title of the paper;
· Name, affiliation and e-mail address for each author;
· A maximum of 5 keywords of the paper;

Also, the name of the workshop session that the paper is being
submitted to should be stated in the email.

Important Dates:

Proposals for special conference sessions and tutorials deadline: 30
December, 2005
Draft Manuscript / Abstract submission deadline: 12 March, 2006
Camera-Ready papers & Pre-registration due: 2 April, 2006
IMECS 2006: 20-22 June, 2006

More details about the IWFE 2006 can be found at:

It will be highly appreciated if you can circulate these calls for
papers to your colleagues.

posted by admin in Uncategorized and have No Comments

New mathematical twists

I would like to inform all of you that there is an interesting book recently published"New Mathematical Twists". This book is about a whole new way of thinking about mathematics. The solutions to the problems in the book will lead you to a hidden treasure. Who knows, the lucky winner may be one of the members of this board. On top of the hidden treasure, you will find many interesting facts about life in general. You can order a copy of the book on this website:

posted by admin in Uncategorized and have No Comments

Exponent manipulation

i am currently programming a numberical solution in matlab.
exponents in the solution result in numbers larger than the
double-percision floating point limit (~10^308).
expo = exp(2*g1*(L+x)/h1)
(L+x)/h1 is very small.
g1 is an array based on the properties of the system. the higher
elements of g1 cause the limit to be exceeded.
These exponents are then formed into a three-dimensional array.
is there any way to manipulate the numbers so that they do not exceed
the floating point limit?

I have tried
expo_1 = exp((2*g1*(L+x)/h1)/n)
expo = expo_1^n
but expo still exceeds the limit and also as n becomes larger, the
percision of the answer drops due to computing errors.


posted by admin in Uncategorized and have Comment (1)

need help: non-linear optimization with linear constraint.

I am totally new to the field of linear/nonlinear programming. I only
have a rough idea about some basic concepts.

The funtion (about 100 variables) to be minimized doesn’t have a close
form. And evaluating it requires some sort of simulation. In another
word, first derivative is not avaiable. There are some linear
constraints on those variables.

Based on my readings, the simple solutions are simplex method or
Powell’s method.  However, I haven’t found any material on how to
handle constraint in those methods.

So, any suggestion, any pointer to readings or implementation will be
helpful. Thanks.

posted by admin in Uncategorized and have Comments (2)

matrix examples

does anyone have any links of any site giving examples of any real or
complex matrices having repeat eigenvalues?


posted by admin in Uncategorized and have Comments (7)

GMRES: what to use for MSOLVE, MATVEC

Please point me to a more appropriate group if this group
isn’t the right one to ask my question.

I grabbed GMRES Fortran code (specifically, single precision
SGMRES) from This code requires 2 user supplied
functions, MSOLVE and MATVEC. I was trying to find out what
functions people typically hook these up to, but searching
for more than one hour with Google didn’t turn up any info
that conclusively answered that question. I also checked the
online book at

Since MSOLVE is solving a system of linear equations, and
MATVEC is computing a matrix-vector product, I think that
I might be able to use the LAPACK routine SGESV for MSOLVE,
and the BLAS2 routine SGEMV for MATVEC.

I have only rudimentary Fortran skills, so I would like to
find out whether I am on the correct path before charging
in that direction. Any insights are much appreciated.

– Norbert

posted by admin in Uncategorized and have Comments (2)

Model reduction for finite element models


I would like to announce my paper that is just appeared online

E. B. Rudnyi and J. G. Korvink. Model Order Reduction for Large Scale
Engineering Models Developed in ANSYS. Lecture Notes in Computer
Science, v. 3732, pp. 349-356, 2005.

If you do not have access to Springer, please contact me directly for

Best wishes,



posted by admin in Uncategorized and have No Comments


Hi !

I would like to code the inverse of erf. Any pointers ?


posted by admin in Uncategorized and have Comments (2)

Slow transport & turn-table-OUT are both UP.

$$           Slow transport & turn-table-OUT are both UP.
$$     "A clock moving OUT [ HORZ ] on a turn-table", moves UP.
$$     Any HORiZONTAL clock "SLOW transport" is *UP* transport.
$$     Where HORiZONTAL means something for a WORLD-line POiNT.
$$     Where OUT and UP means something for a WORLD-line POiNT.
$$     Emphasis Mine,    “`Brian A M Stuckless, Ph.T (Tivity).
p.s.   GR "meanings" are MENTAL; There’s no "PHYSiCAL" meaning.

- Hide quoted text — Show quoted text -

Bilge wrote: >
> Hexenmeister, Princess of Wails, sniffles some more, >
> Hey Sekret Agent Arse Bandit,

>   I still don’t see the androsissy task farce from scotland
> yard surrounding the house. Apparently, they aren’t taking
> their duty to enforce the androsissy usenet tantrum
> protection act seriously enough. When you sober up, why don’t
> you climb onto your tricycle, pedal down to ghq and tell them
> that you are one hopping mad sot and insist that scotland yard
> set up a task force to respond to calls from your shoe phone
> and apprehend the scofflaws who make fun of your tantrums? But,
> that would require sobering up and finding some cajones, so
> nevermind.

>  >"Bilge" <dubi…> wrote in message
>  >…
>  >> Chris:
>  >> >The  Lorenz transformation

>  >> >This contraction is not real it is an effect like perspective
>  >> >because we use light to measure the distance, it has no
>  >> >physical meaning.

>  >>  That is not true. It does have physical meaning and it has
>  >> nothing to do with light.

>  >> >Quite simple really, it is an example of a geometric mean
>  >> >and is called a metric as in differential geometry where you
>  >> >look at a hill, going over the hill is further than along the
>  >> >ground and if you are unaware of the slope the longer distance
>  >> >would be the only way you could know of the hump.
>  >> >Relativity is like that, we cannot see the real path of the
>  >> >light it is distorted by the space distortion caused by the
>  >> >movement of the source or observer.

>  >>  Except for the fact that what an observer actually sees from
>  >> the light rays which arrive simultaneously is not just a length
>  >> contraction.

>  >Fucking delusional cunt Bilge thinks someone has factually
>  > observed a length contraction.
>  >What a moron!

>  >Androcles >  > >  >

>  >> Spheres still look like spheres, not flattened ellipsoids,
>  >> for example.

>  >> >There is no contraction really and no upper limit on velocity.

>  >>  That is like saying that rotated objects aren’t rotated and
>  >> that \pi = 0.

>  >> [...]
>  >> >Much more fun you can get velocities higher than light with a
>  >> >rocket or accelerators where the phase velocity of the
>  >> >travelling wave (as in a rhombadron) is higher than light.
>  >> > Remember group velocityXphase velocity is c for light in a
>  >> >waveguide and the same is true of matter waves so a stationary
>  >> >electrons phase wave occupies the whole universe and the
>  >> >energy quanta of a free electron is influenced by the size of
>  >> >the universe and all the matter in it. >  >> > >  >> >  >> >

Re: The  Lorenz transformation
Re: Slow transport & turn-table-OUT are both UP.   End of POST.

posted by admin in Uncategorized and have No Comments