xref: /trunk/main/scaddins/source/analysis/financial.cxx (revision 91144cd0085a7583d2099b982122deb2184ab956)
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