<t>Hi,Just had a query on Douglas Rachford in 2 space dmiensions...I understand that the scheme is unconditionally stable for 3-D (no cross terms). Is it also stable under 2-D? Or do I have to resort to Peaceman Rachford for the 2-D case?Are there any benefits of using Douglas-Rachford over Peacman ...

Excellent... many thanks

<t>Hi all,Just curious... For a 2-D stochastic vol PDE (first dimension is the underlying, second the vol, e.g Heston)... do people use a set of 'typical' boundary conditions for the vol? For e.g. in 1-D BS, one can make impose a boundary condition vis the delta or the gamma... what could you use fo...

Does anyone understand why there is a 1/q in numerator for the middle line of equation 2.9??I understand why the R_{s}(t) is divided by q (compounding freq) but do not understand why this again needs to be adjusted by 1/q.What am I missing?Thanks,

Thanks folks!Let me look at this a bit more... I will definately report back.Sam

I would have thought that all the default information is incorporated in the 'yield' as you pointed out.In fact, you should be able to use the yield data to work out what the default probability is (assuming an intensiy model for deafult)...I agree with the 0.64 answer; seems rational.Regards,

<t>Hi,Rather than start a new post, I'll hijack this one since my question is very related...Say we have just 2-D (time and space)... and we had a PDE where the coefficients were non constant (e.g. the Black Scholes PDE without making the X = LN(S) substitution)... will the PDE still solve with an i...

<t>Hi,If anyone can shed some light on this then it will be much appreciated.Basically I am trying to understand the limitations of applying the ADI technique as it is written in volume 2 of wilmott. The 3-D PDE has coefficients a,b and in the general case they are functions of t and the 2 underlyin...

It is an excellent conference, although in some of the talks it is quite hard to keep up.

A.real(2) does work and it sets the real part to 2. I glanced over the header file where the class and its members/methods are declared and there isnt an overloaded version of the A.real() function which returns by reference, hence what I was initially trying to do was failing...Thanks,

I knew I was missing something!Many Thanks Lee.Regards,

<t>Hi,Has anyone here used the complex C++ class before Header <complex>? I am trying to do something VERY straightforward but there doesnt seem to be any way to do it!Basically, I define a complex, and then want to re-asign its real and imaginary parts:complex<double> A(1,1); // assign real and im ...

<t>Hi,Thank you all for yoru comments and insights.I have moved a step further and started to implement the integral numerically. I am setting the model to the B-S case for testing purposes (i.e. vol of vol =0) and I am comparing the probabilities from my numerical integration with the closed form v...

Hi Rez,Thanks for the reply.I'm with you right until the end. Where does the 1/2 come from again? How does setting y = infinity change the integral into a 1/2??Regards,

<t>Hi,Can anyone please give me an outline for the derivation of the probability function by inverting its fourier transform, i.e.where f is the characteristic function.Basically, I do not understand where the 1/2 comes from. My approach was to calculate the fourier transform of the probabity functi...