MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Numerical evaluation is Mathematica bottleneck?!

  • To: mathgroup at
  • Subject: [mg69045] Re: Numerical evaluation is Mathematica bottleneck?!
  • From: Jens-Peer Kuska <kuska at>
  • Date: Mon, 28 Aug 2006 02:54:47 -0400 (EDT)
  • References: <ecrasb$p3o$>
  • Sender: owner-wri-mathgroup at


the C/C++ version may be faster -- or not, it depend on
your programming skills and the numerical methods you are
But as long as you don't tell more about the problem and
how you have solved it in Mathematica, there is not defined answer.


gregorc wrote:
> Dear mathgroup members,
> I am solving a non-linear semiconductor PDE problem with finite 
> difference discretization. The Compiled algebraic equation and its 
> Jacobian  take as input central point and two adjacent neighbors and 
> they are Threaded over the whole mesh. A LinearSolve and some damping 
> schemes are used to obtain next step towards the solution. Typical mesh 
> is about 300 points and the solution is usually found in 7 iterations, 
> which takes about 7 seconds of computational time.This is already too 
> slow, but since I plan to solve the problem in 2D this is much too 
> slow!!! The bottleneck is the numerical evaluation of the compiled 
> equation and its Jacobian. I am thinking of using MathLink and write the 
> numerical evaluation routine in C++, which would take a list of data and 
> return evaluated function and the Jacobian.
> Is this the right way? Could the link itself be a timing bottleneck? Any 
> other suggestions?
> Regards,
> Gregor

  • Prev by Date: Re: General--Difficulties in Understanding Mathematica Syntax
  • Next by Date: Re: Re: Subscript vs SubscriptBox
  • Previous by thread: Numerical evaluation is Mathematica bottleneck?!
  • Next by thread: Re: Numerical evaluation is Mathematica bottleneck?!