1 /**************************************************************
2 *
3 * Licensed to the Apache Software Foundation (ASF) under one
4 * or more contributor license agreements. See the NOTICE file
5 * distributed with this work for additional information
6 * regarding copyright ownership. The ASF licenses this file
7 * to you under the Apache License, Version 2.0 (the
8 * "License"); you may not use this file except in compliance
9 * with the License. You may obtain a copy of the License at
10 *
11 * http://www.apache.org/licenses/LICENSE-2.0
12 *
13 * Unless required by applicable law or agreed to in writing,
14 * software distributed under the License is distributed on an
15 * "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
16 * KIND, either express or implied. See the License for the
17 * specific language governing permissions and limitations
18 * under the License.
19 *
20 *************************************************************/
21
22
23
24 #include "analysis.hxx"
25 #include "analysishelper.hxx"
26 #include <rtl/math.hxx>
27
28
29
getAmordegrc(constREFXPS & xOpt,double fCost,sal_Int32 nDate,sal_Int32 nFirstPer,double fRestVal,double fPer,double fRate,const ANY & rOB)30 double SAL_CALL AnalysisAddIn::getAmordegrc( constREFXPS& xOpt,
31 double fCost, sal_Int32 nDate, sal_Int32 nFirstPer, double fRestVal,
32 double fPer, double fRate, const ANY& rOB )
33 {
34 if( nDate > nFirstPer || fRate <= 0.0 || fRestVal > fCost )
35 THROW_IAE;
36
37 double fRet = GetAmordegrc( GetNullDate( xOpt ), fCost, nDate, nFirstPer, fRestVal, fPer, fRate, getDateMode( xOpt, rOB ) );
38 RETURN_FINITE( fRet );
39 }
40
41
getAmorlinc(constREFXPS & xOpt,double fCost,sal_Int32 nDate,sal_Int32 nFirstPer,double fRestVal,double fPer,double fRate,const ANY & rOB)42 double SAL_CALL AnalysisAddIn::getAmorlinc( constREFXPS& xOpt,
43 double fCost, sal_Int32 nDate, sal_Int32 nFirstPer, double fRestVal,
44 double fPer, double fRate, const ANY& rOB )
45 {
46 if( nDate > nFirstPer || fRate <= 0.0 || fRestVal > fCost )
47 THROW_IAE;
48
49 double fRet = GetAmorlinc( GetNullDate( xOpt ), fCost, nDate, nFirstPer, fRestVal, fPer, fRate, getDateMode( xOpt, rOB ) );
50 RETURN_FINITE( fRet );
51 }
52
53
getAccrint(constREFXPS & xOpt,sal_Int32 nIssue,sal_Int32,sal_Int32 nSettle,double fRate,const ANY & rVal,sal_Int32 nFreq,const ANY & rOB)54 double SAL_CALL AnalysisAddIn::getAccrint( constREFXPS& xOpt,
55 sal_Int32 nIssue, sal_Int32 /*nFirstInter*/, sal_Int32 nSettle, double fRate,
56 const ANY &rVal, sal_Int32 nFreq, const ANY& rOB )
57 {
58 double fVal = aAnyConv.getDouble( xOpt, rVal, 1000.0 );
59
60 if( fRate <= 0.0 || fVal <= 0.0 || CHK_Freq || nIssue >= nSettle )
61 THROW_IAE;
62
63 double fRet = fVal * fRate * GetYearDiff( GetNullDate( xOpt ), nIssue, nSettle, getDateMode( xOpt, rOB ) );
64 RETURN_FINITE( fRet );
65 }
66
67
getAccrintm(constREFXPS & xOpt,sal_Int32 nIssue,sal_Int32 nSettle,double fRate,const ANY & rVal,const ANY & rOB)68 double SAL_CALL AnalysisAddIn::getAccrintm( constREFXPS& xOpt,
69 sal_Int32 nIssue, sal_Int32 nSettle, double fRate, const ANY& rVal, const ANY& rOB )
70 {
71 double fVal = aAnyConv.getDouble( xOpt, rVal, 1000.0 );
72
73 if( fRate <= 0.0 || fVal <= 0.0 || nIssue >= nSettle )
74 THROW_IAE;
75
76 double fRet = fVal * fRate * GetYearDiff( GetNullDate( xOpt ), nIssue, nSettle, getDateMode( xOpt, rOB ) );
77 RETURN_FINITE( fRet );
78 }
79
80
getReceived(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fInvest,double fDisc,const ANY & rOB)81 double SAL_CALL AnalysisAddIn::getReceived( constREFXPS& xOpt,
82 sal_Int32 nSettle, sal_Int32 nMat, double fInvest, double fDisc, const ANY& rOB )
83 {
84 if( fInvest <= 0.0 || fDisc <= 0.0 )
85 THROW_IAE;
86
87 double fRet = fInvest / ( 1.0 - ( fDisc * GetYearDiff( GetNullDate( xOpt ), nSettle, nMat, getDateMode( xOpt, rOB ) ) ) );
88 RETURN_FINITE( fRet );
89 }
90
91
getDisc(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fPrice,double fRedemp,const ANY & rOB)92 double SAL_CALL AnalysisAddIn::getDisc( constREFXPS& xOpt,
93 sal_Int32 nSettle, sal_Int32 nMat, double fPrice, double fRedemp, const ANY& rOB )
94 {
95 if( fPrice <= 0.0 || fRedemp <= 0.0 || nSettle >= nMat )
96 THROW_IAE;
97 double fRet = ( 1.0 - fPrice / fRedemp ) / GetYearFrac( xOpt, nSettle, nMat, getDateMode( xOpt, rOB ) );
98 RETURN_FINITE( fRet );
99 }
100
101
getDuration(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fCoup,double fYield,sal_Int32 nFreq,const ANY & rOB)102 double SAL_CALL AnalysisAddIn::getDuration( constREFXPS& xOpt,
103 sal_Int32 nSettle, sal_Int32 nMat, double fCoup, double fYield, sal_Int32 nFreq, const ANY& rOB )
104 {
105 if( fCoup < 0.0 || fYield < 0.0 || CHK_Freq || nSettle >= nMat )
106 THROW_IAE;
107
108 double fRet = GetDuration( GetNullDate( xOpt ), nSettle, nMat, fCoup, fYield, nFreq, getDateMode( xOpt, rOB ) );
109 RETURN_FINITE( fRet );
110 }
111
112
getEffect(double fNominal,sal_Int32 nPeriods)113 double SAL_CALL AnalysisAddIn::getEffect( double fNominal, sal_Int32 nPeriods )
114 {
115 if( nPeriods < 1 || fNominal <= 0.0 )
116 THROW_IAE;
117
118 double fPeriods = nPeriods;
119
120 double fRet = pow( 1.0 + fNominal / fPeriods, fPeriods ) - 1.0;
121 RETURN_FINITE( fRet );
122 }
123
124
getCumprinc(double fRate,sal_Int32 nNumPeriods,double fVal,sal_Int32 nStartPer,sal_Int32 nEndPer,sal_Int32 nPayType)125 double SAL_CALL AnalysisAddIn::getCumprinc( double fRate, sal_Int32 nNumPeriods, double fVal,
126 sal_Int32 nStartPer, sal_Int32 nEndPer, sal_Int32 nPayType )
127 {
128 double fRmz, fKapZ;
129
130 if( nStartPer < 1 || nEndPer < nStartPer || fRate <= 0.0 || nEndPer > nNumPeriods || nNumPeriods <= 0 ||
131 fVal <= 0.0 || ( nPayType != 0 && nPayType != 1 ) )
132 THROW_IAE;
133
134 fRmz = GetRmz( fRate, nNumPeriods, fVal, 0.0, nPayType );
135
136 fKapZ = 0.0;
137
138 sal_uInt32 nStart = sal_uInt32( nStartPer );
139 sal_uInt32 nEnd = sal_uInt32( nEndPer );
140
141 if( nStart == 1 )
142 {
143 if( nPayType <= 0 )
144 fKapZ = fRmz + fVal * fRate;
145 else
146 fKapZ = fRmz;
147
148 nStart++;
149 }
150
151 for( sal_uInt32 i = nStart ; i <= nEnd ; i++ )
152 {
153 if( nPayType > 0 )
154 fKapZ += fRmz - ( GetZw( fRate, double( i - 2 ), fRmz, fVal, 1 ) - fRmz ) * fRate;
155 else
156 fKapZ += fRmz - GetZw( fRate, double( i - 1 ), fRmz, fVal, 0 ) * fRate;
157 }
158
159 RETURN_FINITE( fKapZ );
160 }
161
162
getCumipmt(double fRate,sal_Int32 nNumPeriods,double fVal,sal_Int32 nStartPer,sal_Int32 nEndPer,sal_Int32 nPayType)163 double SAL_CALL AnalysisAddIn::getCumipmt( double fRate, sal_Int32 nNumPeriods, double fVal,
164 sal_Int32 nStartPer, sal_Int32 nEndPer, sal_Int32 nPayType )
165 {
166 double fRmz, fZinsZ;
167
168 if( nStartPer < 1 || nEndPer < nStartPer || fRate <= 0.0 || nEndPer > nNumPeriods || nNumPeriods <= 0 ||
169 fVal <= 0.0 || ( nPayType != 0 && nPayType != 1 ) )
170 THROW_IAE;
171
172 fRmz = GetRmz( fRate, nNumPeriods, fVal, 0.0, nPayType );
173
174 fZinsZ = 0.0;
175
176 sal_uInt32 nStart = sal_uInt32( nStartPer );
177 sal_uInt32 nEnd = sal_uInt32( nEndPer );
178
179 if( nStart == 1 )
180 {
181 if( nPayType <= 0 )
182 fZinsZ = -fVal;
183
184 nStart++;
185 }
186
187 for( sal_uInt32 i = nStart ; i <= nEnd ; i++ )
188 {
189 if( nPayType > 0 )
190 fZinsZ += GetZw( fRate, double( i - 2 ), fRmz, fVal, 1 ) - fRmz;
191 else
192 fZinsZ += GetZw( fRate, double( i - 1 ), fRmz, fVal, 0 );
193 }
194
195 fZinsZ *= fRate;
196
197 RETURN_FINITE( fZinsZ );
198 }
199
200
getPrice(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fRate,double fYield,double fRedemp,sal_Int32 nFreq,const ANY & rOB)201 double SAL_CALL AnalysisAddIn::getPrice( constREFXPS& xOpt,
202 sal_Int32 nSettle, sal_Int32 nMat, double fRate, double fYield, double fRedemp, sal_Int32 nFreq,
203 const ANY& rOB )
204 {
205 if( fYield < 0.0 || fRate < 0.0 || fRedemp <= 0 || CHK_Freq || nSettle >= nMat )
206 THROW_IAE;
207
208 double fRet = getPrice_( GetNullDate( xOpt ), nSettle, nMat, fRate, fYield, fRedemp, nFreq, getDateMode( xOpt, rOB ) );
209 RETURN_FINITE( fRet );
210 }
211
212
getPricedisc(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fDisc,double fRedemp,const ANY & rOB)213 double SAL_CALL AnalysisAddIn::getPricedisc( constREFXPS& xOpt,
214 sal_Int32 nSettle, sal_Int32 nMat, double fDisc, double fRedemp, const ANY& rOB )
215 {
216 if( fDisc <= 0.0 || fRedemp <= 0 || nSettle >= nMat )
217 THROW_IAE;
218
219 double fRet = fRedemp * ( 1.0 - fDisc * GetYearDiff( GetNullDate( xOpt ), nSettle, nMat, getDateMode( xOpt, rOB ) ) );
220 RETURN_FINITE( fRet );
221 }
222
223
getPricemat(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nIssue,double fRate,double fYield,const ANY & rOB)224 double SAL_CALL AnalysisAddIn::getPricemat( constREFXPS& xOpt,
225 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, double fRate, double fYield, const ANY& rOB )
226 {
227 if( fRate < 0.0 || fYield < 0.0 || nSettle >= nMat )
228 THROW_IAE;
229
230 sal_Int32 nNullDate = GetNullDate( xOpt );
231 sal_Int32 nBase = getDateMode( xOpt, rOB );
232
233 double fIssMat = GetYearFrac( nNullDate, nIssue, nMat, nBase );
234 double fIssSet = GetYearFrac( nNullDate, nIssue, nSettle, nBase );
235 double fSetMat = GetYearFrac( nNullDate, nSettle, nMat, nBase );
236
237 double fRet = 1.0 + fIssMat * fRate;
238 fRet /= 1.0 + fSetMat * fYield;
239 fRet -= fIssSet * fRate;
240 fRet *= 100.0;
241
242 RETURN_FINITE( fRet );
243 }
244
245
getMduration(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fCoup,double fYield,sal_Int32 nFreq,const ANY & rOB)246 double SAL_CALL AnalysisAddIn::getMduration( constREFXPS& xOpt,
247 sal_Int32 nSettle, sal_Int32 nMat, double fCoup, double fYield, sal_Int32 nFreq, const ANY& rOB )
248 {
249 if( fCoup < 0.0 || fYield < 0.0 || CHK_Freq )
250 THROW_IAE;
251
252 double fRet = GetDuration( GetNullDate( xOpt ), nSettle, nMat, fCoup, fYield, nFreq, getDateMode( xOpt, rOB ) );
253 fRet /= 1.0 + ( fYield / double( nFreq ) );
254 RETURN_FINITE( fRet );
255 }
256
257
getNominal(double fRate,sal_Int32 nPeriods)258 double SAL_CALL AnalysisAddIn::getNominal( double fRate, sal_Int32 nPeriods )
259 {
260 if( fRate <= 0.0 || nPeriods < 0 )
261 THROW_IAE;
262
263 double fPeriods = nPeriods;
264 double fRet = ( pow( fRate + 1.0, 1.0 / fPeriods ) - 1.0 ) * fPeriods;
265 RETURN_FINITE( fRet );
266 }
267
268
getDollarfr(double fDollarDec,sal_Int32 nFrac)269 double SAL_CALL AnalysisAddIn::getDollarfr( double fDollarDec, sal_Int32 nFrac )
270 {
271 if( nFrac <= 0 )
272 THROW_IAE;
273
274 double fInt;
275 double fFrac = nFrac;
276
277 double fRet = modf( fDollarDec, &fInt );
278
279 fRet *= fFrac;
280
281 fRet *= pow( 10.0, -ceil( log10( fFrac ) ) );
282
283 fRet += fInt;
284
285 RETURN_FINITE( fRet );
286 }
287
288
getDollarde(double fDollarFrac,sal_Int32 nFrac)289 double SAL_CALL AnalysisAddIn::getDollarde( double fDollarFrac, sal_Int32 nFrac )
290 {
291 if( nFrac <= 0 )
292 THROW_IAE;
293
294 double fInt;
295 double fFrac = nFrac;
296
297 double fRet = modf( fDollarFrac, &fInt );
298
299 fRet /= fFrac;
300
301 fRet *= pow( 10.0, ceil( log10( fFrac ) ) );
302
303 fRet += fInt;
304
305 RETURN_FINITE( fRet );
306 }
307
308
getYield(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fCoup,double fPrice,double fRedemp,sal_Int32 nFreq,const ANY & rOB)309 double SAL_CALL AnalysisAddIn::getYield( constREFXPS& xOpt,
310 sal_Int32 nSettle, sal_Int32 nMat, double fCoup, double fPrice, double fRedemp, sal_Int32 nFreq, const ANY& rOB )
311 {
312 if( fCoup < 0.0 || fPrice <= 0.0 || fRedemp <= 0.0 || CHK_Freq || nSettle >= nMat )
313 THROW_IAE;
314
315 double fRet = getYield_( GetNullDate( xOpt ), nSettle, nMat, fCoup, fPrice, fRedemp, nFreq, getDateMode( xOpt, rOB ) );
316 RETURN_FINITE( fRet );
317 }
318
319
getYielddisc(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fPrice,double fRedemp,const ANY & rOB)320 double SAL_CALL AnalysisAddIn::getYielddisc( constREFXPS& xOpt,
321 sal_Int32 nSettle, sal_Int32 nMat, double fPrice, double fRedemp, const ANY& rOB )
322 {
323 if( fPrice <= 0.0 || fRedemp <= 0.0 || nSettle >= nMat )
324 THROW_IAE;
325
326 sal_Int32 nNullDate = GetNullDate( xOpt );
327
328 #if 0
329 double fRet = 1.0 - fPrice / fRedemp;
330 fRet /= GetYearFrac( nNullDate, nSettle, nMat, getDateMode( xOpt, rOB ) );
331 fRet /= 0.99795; // don't know what this constant means in original
332 #endif
333
334 double fRet = ( fRedemp / fPrice ) - 1.0;
335 fRet /= GetYearFrac( nNullDate, nSettle, nMat, getDateMode( xOpt, rOB ) );
336
337 RETURN_FINITE( fRet );
338 }
339
340
getYieldmat(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nIssue,double fRate,double fPrice,const ANY & rOB)341 double SAL_CALL AnalysisAddIn::getYieldmat( constREFXPS& xOpt,
342 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, double fRate, double fPrice, const ANY& rOB )
343 {
344 if( fRate < 0.0 || fRate <= 0.0 || nSettle >= nMat )
345 THROW_IAE;
346
347 double fRet = GetYieldmat( GetNullDate( xOpt ), nSettle, nMat, nIssue, fRate, fPrice, getDateMode( xOpt, rOB ) );
348 RETURN_FINITE( fRet );
349 }
350
351
getTbilleq(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fDisc)352 double SAL_CALL AnalysisAddIn::getTbilleq( constREFXPS& xOpt,
353 sal_Int32 nSettle, sal_Int32 nMat, double fDisc )
354 {
355 nMat++;
356
357 sal_Int32 nDiff = GetDiffDate360( xOpt, nSettle, nMat, sal_True );
358
359 if( fDisc <= 0.0 || nSettle >= nMat || nDiff > 360 )
360 THROW_IAE;
361
362 double fRet = ( 365 * fDisc ) / ( 360 - ( fDisc * double( nDiff ) ) );
363 RETURN_FINITE( fRet );
364 }
365
366
getTbillprice(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fDisc)367 double SAL_CALL AnalysisAddIn::getTbillprice( constREFXPS& xOpt,
368 sal_Int32 nSettle, sal_Int32 nMat, double fDisc )
369 {
370 if( fDisc <= 0.0 || nSettle > nMat )
371 THROW_IAE;
372
373 nMat++;
374
375 double fFraction = GetYearFrac( xOpt, nSettle, nMat, 0 ); // method: USA 30/360
376
377 double fDummy;
378 if( modf( fFraction, &fDummy ) == 0.0 )
379 THROW_IAE;
380
381 double fRet = 100.0 * ( 1.0 - fDisc * fFraction );
382 RETURN_FINITE( fRet );
383 }
384
385
getTbillyield(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fPrice)386 double SAL_CALL AnalysisAddIn::getTbillyield( constREFXPS& xOpt, sal_Int32 nSettle, sal_Int32 nMat, double fPrice )
387 {
388 sal_Int32 nDiff = GetDiffDate360( xOpt, nSettle, nMat, sal_True );
389 nDiff++;
390
391 if( fPrice <= 0.0 || nSettle >= nMat || nDiff > 360 )
392 THROW_IAE;
393
394 double fRet = 100.0;
395 fRet /= fPrice;
396 fRet--;
397 fRet /= double( nDiff );
398 fRet *= 360.0;
399
400 RETURN_FINITE( fRet );
401 }
402
403
getOddfprice(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nIssue,sal_Int32 nFirstCoup,double fRate,double fYield,double fRedemp,sal_Int32 nFreq,const ANY & rOB)404 double SAL_CALL AnalysisAddIn::getOddfprice( constREFXPS& xOpt,
405 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, sal_Int32 nFirstCoup,
406 double fRate, double fYield, double fRedemp, sal_Int32 nFreq, const ANY& rOB )
407 {
408 if( fRate < 0 || fYield < 0 || CHK_Freq || nMat <= nFirstCoup || nFirstCoup <= nSettle || nSettle <= nIssue )
409 THROW_IAE;
410
411 double fRet = GetOddfprice( GetNullDate( xOpt ), nSettle, nMat, nIssue, nFirstCoup, fRate, fYield, fRedemp, nFreq, getDateMode( xOpt, rOB ) );
412 RETURN_FINITE( fRet );
413 }
414
415
getOddfyield(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nIssue,sal_Int32 nFirstCoup,double fRate,double fPrice,double fRedemp,sal_Int32 nFreq,const ANY & rOB)416 double SAL_CALL AnalysisAddIn::getOddfyield( constREFXPS& xOpt,
417 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, sal_Int32 nFirstCoup,
418 double fRate, double fPrice, double fRedemp, sal_Int32 nFreq, const ANY& rOB )
419 {
420 if( fRate < 0 || fPrice <= 0 || CHK_Freq || nMat <= nFirstCoup || nFirstCoup <= nSettle || nSettle <= nIssue )
421 THROW_IAE;
422
423 double fRet = GetOddfyield( GetNullDate( xOpt ), nSettle, nMat, nIssue, nFirstCoup, fRate, fPrice, fRedemp, nFreq,
424 getDateMode( xOpt, rOB ) );
425 RETURN_FINITE( fRet );
426 }
427
428
getOddlprice(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nLastInterest,double fRate,double fYield,double fRedemp,sal_Int32 nFreq,const ANY & rOB)429 double SAL_CALL AnalysisAddIn::getOddlprice( constREFXPS& xOpt,
430 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nLastInterest,
431 double fRate, double fYield, double fRedemp, sal_Int32 nFreq, const ANY& rOB )
432 {
433 if( fRate < 0 || fYield < 0 || CHK_Freq || nMat <= nSettle || nSettle <= nLastInterest )
434 THROW_IAE;
435
436 double fRet = GetOddlprice( GetNullDate( xOpt ), nSettle, nMat, nLastInterest, fRate, fYield, fRedemp, nFreq,
437 getDateMode( xOpt, rOB ) );
438 RETURN_FINITE( fRet );
439 }
440
441
getOddlyield(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nLastInterest,double fRate,double fPrice,double fRedemp,sal_Int32 nFreq,const ANY & rOB)442 double SAL_CALL AnalysisAddIn::getOddlyield( constREFXPS& xOpt,
443 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nLastInterest,
444 double fRate, double fPrice, double fRedemp, sal_Int32 nFreq, const ANY& rOB )
445 {
446 if( fRate < 0 || fPrice <= 0 || CHK_Freq || nMat <= nSettle || nSettle <= nLastInterest )
447 THROW_IAE;
448
449 double fRet = GetOddlyield( GetNullDate( xOpt ), nSettle, nMat, nLastInterest, fRate, fPrice, fRedemp, nFreq,
450 getDateMode( xOpt, rOB ) );
451 RETURN_FINITE( fRet );
452 }
453
454
455 // ============================================================================
456 // XIRR helper functions
457
458 #define V_(i) (*rValues.Get(i))
459 #define D_(i) (*rDates.Get(i))
460
461 /** Calculates the resulting amount for the passed interest rate and the given XIRR parameters. */
lcl_sca_XirrResult(const ScaDoubleList & rValues,const ScaDoubleList & rDates,double fRate)462 double lcl_sca_XirrResult( const ScaDoubleList& rValues, const ScaDoubleList& rDates, double fRate )
463 {
464 /* V_0 ... V_n = input values.
465 D_0 ... D_n = input dates.
466 R = input interest rate.
467
468 r := R+1
469 E_i := (D_i-D_0) / 365
470
471 n V_i n V_i
472 f(R) = SUM ------- = V_0 + SUM ------- .
473 i=0 r^E_i i=1 r^E_i
474 */
475 double D_0 = D_(0);
476 double r = fRate + 1.0;
477 double fResult = V_(0);
478 for( sal_uInt32 i = 1, nCount = rValues.Count(); i < nCount; ++i )
479 fResult += V_(i) / pow( r, (D_(i) - D_0) / 365.0 );
480 return fResult;
481 }
482
483 /** Calculates the first derivation of lcl_sca_XirrResult(). */
lcl_sca_XirrResult_Deriv1(const ScaDoubleList & rValues,const ScaDoubleList & rDates,double fRate)484 double lcl_sca_XirrResult_Deriv1( const ScaDoubleList& rValues, const ScaDoubleList& rDates, double fRate )
485 {
486 /* V_0 ... V_n = input values.
487 D_0 ... D_n = input dates.
488 R = input interest rate.
489
490 r := R+1
491 E_i := (D_i-D_0) / 365
492
493 n V_i
494 f'(R) = [ V_0 + SUM ------- ]'
495 i=1 r^E_i
496
497 n V_i n E_i V_i
498 = 0 + SUM -E_i ----------- r' = - SUM ----------- .
499 i=1 r^(E_i+1) i=1 r^(E_i+1)
500 */
501 double D_0 = D_(0);
502 double r = fRate + 1.0;
503 double fResult = 0.0;
504 for( sal_uInt32 i = 1, nCount = rValues.Count(); i < nCount; ++i )
505 {
506 double E_i = (D_(i) - D_0) / 365.0;
507 fResult -= E_i * V_(i) / pow( r, E_i + 1.0 );
508 }
509 return fResult;
510 }
511
512 #undef V_
513 #undef D_
514
515
516 // ----------------------------------------------------------------------------
517 // XIRR calculation
518
getXirr(constREFXPS & xOpt,const SEQSEQ (double)& rValues,const SEQSEQ (sal_Int32)& rDates,const ANY & rGuessRate)519 double SAL_CALL AnalysisAddIn::getXirr(
520 constREFXPS& xOpt, const SEQSEQ( double )& rValues, const SEQSEQ( sal_Int32 )& rDates, const ANY& rGuessRate )
521 {
522 ScaDoubleList aValues, aDates;
523 aValues.Append( rValues );
524 aDates.Append( rDates );
525
526 if( (aValues.Count() < 2) || (aValues.Count() != aDates.Count()) )
527 THROW_IAE;
528
529 // result interest rate, initialized with passed guessed rate, or 10%
530 double fResultRate = aAnyConv.getDouble( xOpt, rGuessRate, 0.1 );
531 if( fResultRate <= -1 )
532 THROW_IAE;
533
534 // maximum epsilon for end of iteration
535 static const double fMaxEps = 1e-10;
536 // maximum number of iterations
537 static const sal_Int32 nMaxIter = 50;
538
539 // Newton's method - try to find a fResultRate, so that lcl_sca_XirrResult() returns 0.
540 double fNewRate, fRateEps, fResultValue;
541 sal_Int32 nIter = 0;
542 sal_Int32 nIterScan = 0;
543 bool bContLoop = false;
544 bool bResultRateScanEnd = false;
545
546 // First the inner while-loop will be executed using the default Value fResultRate
547 // or the user guessed fResultRate if those do not deliver a solution for the
548 // Newton's method then the range from -0.99 to +0.99 will be scanned with a
549 // step size of 0.01 to find fResultRate's value which can deliver a solution
550 do
551 {
552 if (nIterScan >=1)
553 fResultRate = -0.99 + (nIterScan -1)* 0.01;
554 do
555 {
556 fResultValue = lcl_sca_XirrResult( aValues, aDates, fResultRate );
557 fNewRate = fResultRate - fResultValue / lcl_sca_XirrResult_Deriv1( aValues, aDates, fResultRate );
558 fRateEps = fabs( fNewRate - fResultRate );
559 fResultRate = fNewRate;
560 bContLoop = (fRateEps > fMaxEps) && (fabs( fResultValue ) > fMaxEps);
561 }
562 while( bContLoop && (++nIter < nMaxIter) );
563 nIter = 0;
564 if ( ::rtl::math::isNan(fResultRate) || ::rtl::math::isInf(fResultRate)
565 ||::rtl::math::isNan(fResultValue) || ::rtl::math::isInf(fResultValue))
566 bContLoop = true;
567
568 ++nIterScan;
569 bResultRateScanEnd = (nIterScan >= 200);
570 }
571 while(bContLoop && !bResultRateScanEnd);
572
573 if( bContLoop )
574 THROW_IAE;
575 RETURN_FINITE( fResultRate );
576 }
577
578
579 // ============================================================================
580
getXnpv(double fRate,const SEQSEQ (double)& rValues,const SEQSEQ (sal_Int32)& rDates)581 double SAL_CALL AnalysisAddIn::getXnpv(
582 double fRate, const SEQSEQ( double )& rValues, const SEQSEQ( sal_Int32 )& rDates )
583 {
584 ScaDoubleList aValList;
585 ScaDoubleList aDateList;
586
587 aValList.Append( rValues );
588 aDateList.Append( rDates );
589
590 sal_Int32 nNum = aValList.Count();
591
592 if( nNum != sal_Int32( aDateList.Count() ) || nNum < 2 )
593 THROW_IAE;
594
595 double fRet = 0.0;
596 double fNull = *aDateList.Get( 0 );
597 fRate++;
598
599 for( sal_Int32 i = 0 ; i < nNum ; i++ )
600 fRet += *aValList.Get( i ) / ( pow( fRate, ( *aDateList.Get( i ) - fNull ) / 365.0 ) );
601
602 RETURN_FINITE( fRet );
603 }
604
605
getIntrate(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fInvest,double fRedemp,const ANY & rOB)606 double SAL_CALL AnalysisAddIn::getIntrate( constREFXPS& xOpt,
607 sal_Int32 nSettle, sal_Int32 nMat, double fInvest, double fRedemp, const ANY& rOB )
608 {
609 if( fInvest <= 0.0 || fRedemp <= 0.0 || nSettle >= nMat )
610 THROW_IAE;
611
612 double fRet = ( ( fRedemp / fInvest ) - 1.0 ) / GetYearDiff( GetNullDate( xOpt ), nSettle, nMat, getDateMode( xOpt, rOB ) );
613 RETURN_FINITE( fRet );
614 }
615
616
getCoupncd(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)617 double SAL_CALL AnalysisAddIn::getCoupncd( constREFXPS& xOpt,
618 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB )
619 {
620 double fRet = GetCoupncd( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
621 RETURN_FINITE( fRet );
622 }
623
624
getCoupdays(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)625 double SAL_CALL AnalysisAddIn::getCoupdays( constREFXPS& xOpt,
626 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB )
627 {
628 double fRet = GetCoupdays( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
629 RETURN_FINITE( fRet );
630 }
631
632
getCoupdaysnc(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)633 double SAL_CALL AnalysisAddIn::getCoupdaysnc( constREFXPS& xOpt,
634 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB )
635 {
636 double fRet = GetCoupdaysnc( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
637 RETURN_FINITE( fRet );
638 }
639
640
getCoupdaybs(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)641 double SAL_CALL AnalysisAddIn::getCoupdaybs( constREFXPS& xOpt,
642 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB )
643 {
644 double fRet = GetCoupdaybs( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
645 RETURN_FINITE( fRet );
646 }
647
648
getCouppcd(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)649 double SAL_CALL AnalysisAddIn::getCouppcd( constREFXPS& xOpt,
650 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB )
651 {
652 double fRet = GetCouppcd( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
653 RETURN_FINITE( fRet );
654 }
655
656
getCoupnum(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)657 double SAL_CALL AnalysisAddIn::getCoupnum( constREFXPS& xOpt,
658 sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB )
659 {
660 double fRet = GetCoupnum( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
661 RETURN_FINITE( fRet );
662 }
663
664
getFvschedule(double fPrinc,const SEQSEQ (double)& rSchedule)665 double SAL_CALL AnalysisAddIn::getFvschedule( double fPrinc, const SEQSEQ( double )& rSchedule )
666 {
667 ScaDoubleList aSchedList;
668
669 aSchedList.Append( rSchedule );
670
671 for( const double* p = aSchedList.First() ; p ; p = aSchedList.Next() )
672 fPrinc *= 1.0 + *p;
673
674 RETURN_FINITE( fPrinc );
675 }
676