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 ) THROWDEF_RTE_IAE
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 ) THROWDEF_RTE_IAE
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 ) THROWDEF_RTE_IAE
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 ) THROWDEF_RTE_IAE
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 ) THROWDEF_RTE_IAE
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 ) THROWDEF_RTE_IAE
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 	THROWDEF_RTE_IAE
105 {
106 	if( fCoup < 0.0 || fYield < 0.0 || CHK_Freq || nSettle >= nMat )
107 		THROW_IAE;
108 
109     double fRet = GetDuration( GetNullDate( xOpt ),  nSettle, nMat, fCoup, fYield, nFreq, getDateMode( xOpt, rOB ) );
110     RETURN_FINITE( fRet );
111 }
112 
113 
getEffect(double fNominal,sal_Int32 nPeriods)114 double SAL_CALL AnalysisAddIn::getEffect( double fNominal, sal_Int32 nPeriods ) THROWDEF_RTE_IAE
115 {
116 	if( nPeriods < 1 || fNominal <= 0.0 )
117 		THROW_IAE;
118 
119 	double	fPeriods = nPeriods;
120 
121     double fRet = pow( 1.0 + fNominal / fPeriods, fPeriods ) - 1.0;
122     RETURN_FINITE( fRet );
123 }
124 
125 
getCumprinc(double fRate,sal_Int32 nNumPeriods,double fVal,sal_Int32 nStartPer,sal_Int32 nEndPer,sal_Int32 nPayType)126 double SAL_CALL AnalysisAddIn::getCumprinc( double fRate, sal_Int32 nNumPeriods, double fVal,
127 	sal_Int32 nStartPer, sal_Int32 nEndPer, sal_Int32 nPayType ) THROWDEF_RTE_IAE
128 {
129 	double fRmz, fKapZ;
130 
131 	if( nStartPer < 1 || nEndPer < nStartPer || fRate <= 0.0 || nEndPer > nNumPeriods  || nNumPeriods <= 0 ||
132 		fVal <= 0.0 || ( nPayType != 0 && nPayType != 1 ) )
133 		THROW_IAE;
134 
135 	fRmz = GetRmz( fRate, nNumPeriods, fVal, 0.0, nPayType );
136 
137 	fKapZ = 0.0;
138 
139     sal_uInt32  nStart = sal_uInt32( nStartPer );
140     sal_uInt32  nEnd = sal_uInt32( nEndPer );
141 
142     if( nStart == 1 )
143 	{
144 		if( nPayType <= 0 )
145 			fKapZ = fRmz + fVal * fRate;
146 		else
147 			fKapZ = fRmz;
148 
149         nStart++;
150 	}
151 
152     for( sal_uInt32 i = nStart ; i <= nEnd ; i++ )
153 	{
154 		if( nPayType > 0 )
155 			fKapZ += fRmz - ( GetZw( fRate, double( i - 2 ), fRmz, fVal, 1 ) - fRmz ) * fRate;
156 		else
157 			fKapZ += fRmz - GetZw( fRate, double( i - 1 ), fRmz, fVal, 0 ) * fRate;
158 	}
159 
160     RETURN_FINITE( fKapZ );
161 }
162 
163 
getCumipmt(double fRate,sal_Int32 nNumPeriods,double fVal,sal_Int32 nStartPer,sal_Int32 nEndPer,sal_Int32 nPayType)164 double SAL_CALL AnalysisAddIn::getCumipmt( double fRate, sal_Int32 nNumPeriods, double fVal,
165 	sal_Int32 nStartPer, sal_Int32 nEndPer, sal_Int32 nPayType ) THROWDEF_RTE_IAE
166 {
167 	double fRmz, fZinsZ;
168 
169 	if( nStartPer < 1 || nEndPer < nStartPer || fRate <= 0.0 || nEndPer > nNumPeriods  || nNumPeriods <= 0 ||
170 		fVal <= 0.0 || ( nPayType != 0 && nPayType != 1 ) )
171 		THROW_IAE;
172 
173 	fRmz = GetRmz( fRate, nNumPeriods, fVal, 0.0, nPayType );
174 
175 	fZinsZ = 0.0;
176 
177     sal_uInt32  nStart = sal_uInt32( nStartPer );
178     sal_uInt32  nEnd = sal_uInt32( nEndPer );
179 
180     if( nStart == 1 )
181 	{
182 		if( nPayType <= 0 )
183 			fZinsZ = -fVal;
184 
185         nStart++;
186 	}
187 
188     for( sal_uInt32 i = nStart ; i <= nEnd ; i++ )
189 	{
190 		if( nPayType > 0 )
191 			fZinsZ += GetZw( fRate, double( i - 2 ), fRmz, fVal, 1 ) - fRmz;
192 		else
193 			fZinsZ += GetZw( fRate, double( i - 1 ), fRmz, fVal, 0 );
194 	}
195 
196 	fZinsZ *= fRate;
197 
198     RETURN_FINITE( fZinsZ );
199 }
200 
201 
getPrice(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fRate,double fYield,double fRedemp,sal_Int32 nFreq,const ANY & rOB)202 double SAL_CALL AnalysisAddIn::getPrice( constREFXPS& xOpt,
203 	sal_Int32 nSettle, sal_Int32 nMat, double fRate, double fYield, double fRedemp, sal_Int32 nFreq,
204 	const ANY& rOB ) THROWDEF_RTE_IAE
205 {
206 	if( fYield < 0.0 || fRate < 0.0 || fRedemp <= 0 || CHK_Freq || nSettle >= nMat )
207 		THROW_IAE;
208 
209     double fRet = getPrice_( GetNullDate( xOpt ), nSettle, nMat, fRate, fYield, fRedemp, nFreq, getDateMode( xOpt, rOB ) );
210     RETURN_FINITE( fRet );
211 }
212 
213 
getPricedisc(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fDisc,double fRedemp,const ANY & rOB)214 double SAL_CALL AnalysisAddIn::getPricedisc( constREFXPS& xOpt,
215 	sal_Int32 nSettle, sal_Int32 nMat, double fDisc, double fRedemp, const ANY& rOB ) THROWDEF_RTE_IAE
216 {
217 	if( fDisc <= 0.0 || fRedemp <= 0 || nSettle >= nMat )
218 		THROW_IAE;
219 
220     double fRet = fRedemp * ( 1.0 - fDisc * GetYearDiff( GetNullDate( xOpt ), nSettle, nMat, getDateMode( xOpt, rOB ) ) );
221     RETURN_FINITE( fRet );
222 }
223 
224 
getPricemat(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nIssue,double fRate,double fYield,const ANY & rOB)225 double SAL_CALL AnalysisAddIn::getPricemat( constREFXPS& xOpt,
226 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, double fRate, double fYield, const ANY& rOB )
227 	THROWDEF_RTE_IAE
228 {
229 	if( fRate < 0.0 || fYield < 0.0 || nSettle >= nMat )
230 		THROW_IAE;
231 
232 	sal_Int32	nNullDate = GetNullDate( xOpt );
233     sal_Int32   nBase = getDateMode( xOpt, rOB );
234 
235 	double		fIssMat = GetYearFrac( nNullDate, nIssue, nMat, nBase );
236 	double		fIssSet = GetYearFrac( nNullDate, nIssue, nSettle, nBase );
237 	double		fSetMat = GetYearFrac( nNullDate, nSettle, nMat, nBase );
238 
239 	double		fRet = 1.0 + fIssMat * fRate;
240 	fRet /= 1.0 + fSetMat * fYield;
241 	fRet -= fIssSet * fRate;
242 	fRet *= 100.0;
243 
244     RETURN_FINITE( fRet );
245 }
246 
247 
getMduration(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fCoup,double fYield,sal_Int32 nFreq,const ANY & rOB)248 double SAL_CALL AnalysisAddIn::getMduration( constREFXPS& xOpt,
249 	sal_Int32 nSettle, sal_Int32 nMat, double fCoup, double fYield, sal_Int32 nFreq, const ANY& rOB )
250 	THROWDEF_RTE_IAE
251 {
252 	if( fCoup < 0.0 || fYield < 0.0 || CHK_Freq )
253 		THROW_IAE;
254 
255     double      fRet = GetDuration( GetNullDate( xOpt ),  nSettle, nMat, fCoup, fYield, nFreq, getDateMode( xOpt, rOB ) );
256 	fRet /= 1.0 + ( fYield / double( nFreq ) );
257     RETURN_FINITE( fRet );
258 }
259 
260 
getNominal(double fRate,sal_Int32 nPeriods)261 double SAL_CALL AnalysisAddIn::getNominal( double fRate, sal_Int32 nPeriods ) THROWDEF_RTE_IAE
262 {
263 	if( fRate <= 0.0 || nPeriods < 0 )
264 		THROW_IAE;
265 
266 	double	fPeriods = nPeriods;
267     double fRet = ( pow( fRate + 1.0, 1.0 / fPeriods ) - 1.0 ) * fPeriods;
268     RETURN_FINITE( fRet );
269 }
270 
271 
getDollarfr(double fDollarDec,sal_Int32 nFrac)272 double SAL_CALL AnalysisAddIn::getDollarfr( double fDollarDec, sal_Int32 nFrac ) THROWDEF_RTE_IAE
273 {
274 	if( nFrac <= 0 )
275 		THROW_IAE;
276 
277 	double	fInt;
278 	double	fFrac = nFrac;
279 
280 	double	fRet = modf( fDollarDec, &fInt );
281 
282 	fRet *= fFrac;
283 
284 	fRet *= pow( 10.0, -ceil( log10( fFrac ) ) );
285 
286 	fRet += fInt;
287 
288     RETURN_FINITE( fRet );
289 }
290 
291 
getDollarde(double fDollarFrac,sal_Int32 nFrac)292 double SAL_CALL AnalysisAddIn::getDollarde( double fDollarFrac, sal_Int32 nFrac ) THROWDEF_RTE_IAE
293 {
294 	if( nFrac <= 0 )
295 		THROW_IAE;
296 
297 	double	fInt;
298 	double	fFrac = nFrac;
299 
300 	double	fRet = modf( fDollarFrac, &fInt );
301 
302 	fRet /= fFrac;
303 
304 	fRet *= pow( 10.0, ceil( log10( fFrac ) ) );
305 
306 	fRet += fInt;
307 
308     RETURN_FINITE( fRet );
309 }
310 
311 
getYield(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fCoup,double fPrice,double fRedemp,sal_Int32 nFreq,const ANY & rOB)312 double SAL_CALL AnalysisAddIn::getYield( constREFXPS& xOpt,
313 	sal_Int32 nSettle, sal_Int32 nMat, double fCoup, double fPrice, double fRedemp, sal_Int32 nFreq, const ANY& rOB )
314 	THROWDEF_RTE_IAE
315 {
316 	if( fCoup < 0.0 || fPrice <= 0.0 || fRedemp <= 0.0 || CHK_Freq || nSettle >= nMat )
317 		THROW_IAE;
318 
319     double fRet = getYield_( GetNullDate( xOpt ), nSettle, nMat, fCoup, fPrice, fRedemp, nFreq, getDateMode( xOpt, rOB ) );
320     RETURN_FINITE( fRet );
321 }
322 
323 
getYielddisc(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fPrice,double fRedemp,const ANY & rOB)324 double SAL_CALL AnalysisAddIn::getYielddisc( constREFXPS& xOpt,
325 	sal_Int32 nSettle, sal_Int32 nMat, double fPrice, double fRedemp, const ANY& rOB ) THROWDEF_RTE_IAE
326 {
327 	if( fPrice <= 0.0 || fRedemp <= 0.0 || nSettle >= nMat )
328 		THROW_IAE;
329 
330 	sal_Int32	nNullDate = GetNullDate( xOpt );
331 
332 #if 0
333 	double		fRet = 1.0 - fPrice / fRedemp;
334     fRet /= GetYearFrac( nNullDate, nSettle, nMat, getDateMode( xOpt, rOB ) );
335     fRet /= 0.99795;  // don't know what this constant means in original
336 #endif
337 
338     double fRet = ( fRedemp / fPrice ) - 1.0;
339     fRet /= GetYearFrac( nNullDate, nSettle, nMat, getDateMode( xOpt, rOB ) );
340 
341     RETURN_FINITE( fRet );
342 }
343 
344 
getYieldmat(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nIssue,double fRate,double fPrice,const ANY & rOB)345 double SAL_CALL AnalysisAddIn::getYieldmat( constREFXPS& xOpt,
346 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, double fRate, double fPrice, const ANY& rOB )
347 	THROWDEF_RTE_IAE
348 {
349 	if( fRate < 0.0 || fRate <= 0.0 || nSettle >= nMat )
350 		THROW_IAE;
351 
352     double fRet = GetYieldmat( GetNullDate( xOpt ),  nSettle, nMat, nIssue, fRate, fPrice, getDateMode( xOpt, rOB ) );
353     RETURN_FINITE( fRet );
354 }
355 
356 
getTbilleq(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fDisc)357 double SAL_CALL AnalysisAddIn::getTbilleq( constREFXPS& xOpt,
358 	sal_Int32 nSettle, sal_Int32 nMat, double fDisc ) THROWDEF_RTE_IAE
359 {
360 	nMat++;
361 
362 	sal_Int32	nDiff = GetDiffDate360( xOpt, nSettle, nMat, sal_True );
363 
364 	if( fDisc <= 0.0 || nSettle >= nMat || nDiff > 360 )
365 		THROW_IAE;
366 
367     double fRet = ( 365 * fDisc ) / ( 360 - ( fDisc * double( nDiff ) ) );
368     RETURN_FINITE( fRet );
369 }
370 
371 
getTbillprice(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fDisc)372 double SAL_CALL AnalysisAddIn::getTbillprice( constREFXPS& xOpt,
373 	sal_Int32 nSettle, sal_Int32 nMat, double fDisc ) THROWDEF_RTE_IAE
374 {
375 	if( fDisc <= 0.0 || nSettle > nMat )
376 		THROW_IAE;
377 
378 	nMat++;
379 
380 	double	fFraction = GetYearFrac( xOpt, nSettle, nMat, 0 );	// method: USA 30/360
381 
382 	double	fDummy;
383 	if( modf( fFraction, &fDummy ) == 0.0 )
384 		THROW_IAE;
385 
386     double fRet = 100.0 * ( 1.0 - fDisc * fFraction );
387     RETURN_FINITE( fRet );
388 }
389 
390 
getTbillyield(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fPrice)391 double SAL_CALL AnalysisAddIn::getTbillyield( constREFXPS& xOpt, sal_Int32 nSettle, sal_Int32 nMat, double fPrice )
392 	THROWDEF_RTE_IAE
393 {
394 	sal_Int32	nDiff = GetDiffDate360( xOpt, nSettle, nMat, sal_True );
395 	nDiff++;
396 
397 	if( fPrice <= 0.0 || nSettle >= nMat || nDiff > 360 )
398 		THROW_IAE;
399 
400 	double		fRet = 100.0;
401 	fRet /= fPrice;
402 	fRet--;
403 	fRet /= double( nDiff );
404 	fRet *= 360.0;
405 
406     RETURN_FINITE( fRet );
407 }
408 
409 
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)410 double SAL_CALL AnalysisAddIn::getOddfprice( constREFXPS& xOpt,
411 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, sal_Int32 nFirstCoup,
412 	double fRate, double fYield, double fRedemp, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
413 {
414 	if( fRate < 0 || fYield < 0 || CHK_Freq || nMat <= nFirstCoup || nFirstCoup <= nSettle || nSettle <= nIssue )
415 		THROW_IAE;
416 
417     double fRet = GetOddfprice( GetNullDate( xOpt ), nSettle, nMat, nIssue, nFirstCoup, fRate, fYield, fRedemp, nFreq, getDateMode( xOpt, rOB ) );
418     RETURN_FINITE( fRet );
419 }
420 
421 
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)422 double SAL_CALL AnalysisAddIn::getOddfyield( constREFXPS& xOpt,
423 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nIssue, sal_Int32 nFirstCoup,
424 	double fRate, double fPrice, double fRedemp, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
425 {
426 	if( fRate < 0 || fPrice <= 0 || CHK_Freq || nMat <= nFirstCoup || nFirstCoup <= nSettle || nSettle <= nIssue )
427 		THROW_IAE;
428 
429     double fRet = GetOddfyield( GetNullDate( xOpt ), nSettle, nMat, nIssue, nFirstCoup, fRate, fPrice, fRedemp, nFreq,
430                         getDateMode( xOpt, rOB ) );
431     RETURN_FINITE( fRet );
432 }
433 
434 
getOddlprice(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nLastInterest,double fRate,double fYield,double fRedemp,sal_Int32 nFreq,const ANY & rOB)435 double SAL_CALL AnalysisAddIn::getOddlprice( constREFXPS& xOpt,
436 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nLastInterest,
437 	double fRate, double fYield, double fRedemp, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
438 {
439 	if( fRate < 0 || fYield < 0 || CHK_Freq || nMat <= nSettle || nSettle <= nLastInterest )
440 		THROW_IAE;
441 
442     double fRet = GetOddlprice( GetNullDate( xOpt ), nSettle, nMat, nLastInterest, fRate, fYield, fRedemp, nFreq,
443                         getDateMode( xOpt, rOB ) );
444     RETURN_FINITE( fRet );
445 }
446 
447 
getOddlyield(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nLastInterest,double fRate,double fPrice,double fRedemp,sal_Int32 nFreq,const ANY & rOB)448 double SAL_CALL AnalysisAddIn::getOddlyield( constREFXPS& xOpt,
449 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nLastInterest,
450 	double fRate, double fPrice, double fRedemp, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
451 {
452 	if( fRate < 0 || fPrice <= 0 || CHK_Freq || nMat <= nSettle || nSettle <= nLastInterest )
453 		THROW_IAE;
454 
455     double fRet = GetOddlyield( GetNullDate( xOpt ), nSettle, nMat, nLastInterest, fRate, fPrice, fRedemp, nFreq,
456                         getDateMode( xOpt, rOB ) );
457     RETURN_FINITE( fRet );
458 }
459 
460 
461 // ============================================================================
462 // XIRR helper functions
463 
464 #define V_(i) (*rValues.Get(i))
465 #define D_(i) (*rDates.Get(i))
466 
467 /** 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)468 double lcl_sca_XirrResult( const ScaDoubleList& rValues, const ScaDoubleList& rDates, double fRate )
469 {
470     /*  V_0 ... V_n = input values.
471         D_0 ... D_n = input dates.
472         R           = input interest rate.
473 
474         r   := R+1
475         E_i := (D_i-D_0) / 365
476 
477                     n    V_i                n    V_i
478         f(R)  =  SUM   -------  =  V_0 + SUM   ------- .
479                    i=0  r^E_i              i=1  r^E_i
480     */
481     double D_0 = D_(0);
482     double r = fRate + 1.0;
483     double fResult = V_(0);
484     for( sal_uInt32 i = 1, nCount = rValues.Count(); i < nCount; ++i )
485         fResult += V_(i) / pow( r, (D_(i) - D_0) / 365.0 );
486     return fResult;
487 }
488 
489 /** Calculates the first derivation of lcl_sca_XirrResult(). */
lcl_sca_XirrResult_Deriv1(const ScaDoubleList & rValues,const ScaDoubleList & rDates,double fRate)490 double lcl_sca_XirrResult_Deriv1( const ScaDoubleList& rValues, const ScaDoubleList& rDates, double fRate )
491 {
492     /*  V_0 ... V_n = input values.
493         D_0 ... D_n = input dates.
494         R           = input interest rate.
495 
496         r   := R+1
497         E_i := (D_i-D_0) / 365
498 
499                              n    V_i
500         f'(R)  =  [ V_0 + SUM   ------- ]'
501                             i=1  r^E_i
502 
503                          n           V_i                 n    E_i V_i
504                =  0 + SUM   -E_i ----------- r'  =  - SUM   ----------- .
505                         i=1       r^(E_i+1)             i=1  r^(E_i+1)
506     */
507     double D_0 = D_(0);
508     double r = fRate + 1.0;
509     double fResult = 0.0;
510     for( sal_uInt32 i = 1, nCount = rValues.Count(); i < nCount; ++i )
511     {
512         double E_i = (D_(i) - D_0) / 365.0;
513         fResult -= E_i * V_(i) / pow( r, E_i + 1.0 );
514     }
515     return fResult;
516 }
517 
518 #undef V_
519 #undef D_
520 
521 
522 // ----------------------------------------------------------------------------
523 // XIRR calculation
524 
getXirr(constREFXPS & xOpt,const SEQSEQ (double)& rValues,const SEQSEQ (sal_Int32)& rDates,const ANY & rGuessRate)525 double SAL_CALL AnalysisAddIn::getXirr(
526     constREFXPS& xOpt, const SEQSEQ( double )& rValues, const SEQSEQ( sal_Int32 )& rDates, const ANY& rGuessRate ) THROWDEF_RTE_IAE
527 {
528     ScaDoubleList aValues, aDates;
529     aValues.Append( rValues );
530     aDates.Append( rDates );
531 
532     if( (aValues.Count() < 2) || (aValues.Count() != aDates.Count()) )
533         THROW_IAE;
534 
535     // result interest rate, initialized with passed guessed rate, or 10%
536     double fResultRate = aAnyConv.getDouble( xOpt, rGuessRate, 0.1 );
537     if( fResultRate <= -1 )
538         THROW_IAE;
539 
540     // maximum epsilon for end of iteration
541     static const double fMaxEps = 1e-10;
542     // maximum number of iterations
543     static const sal_Int32 nMaxIter = 50;
544 
545     // Newton's method - try to find a fResultRate, so that lcl_sca_XirrResult() returns 0.
546     double fNewRate, fRateEps, fResultValue;
547     sal_Int32 nIter = 0;
548     bool bContLoop;
549     do
550     {
551         fResultValue = lcl_sca_XirrResult( aValues, aDates, fResultRate );
552         fNewRate = fResultRate - fResultValue / lcl_sca_XirrResult_Deriv1( aValues, aDates, fResultRate );
553         fRateEps = fabs( fNewRate - fResultRate );
554         fResultRate = fNewRate;
555         bContLoop = (fRateEps > fMaxEps) && (fabs( fResultValue ) > fMaxEps);
556     }
557     while( bContLoop && (++nIter < nMaxIter) );
558 
559     if( bContLoop )
560         THROW_IAE;
561     RETURN_FINITE( fResultRate );
562 }
563 
564 
565 // ============================================================================
566 
getXnpv(double fRate,const SEQSEQ (double)& rValues,const SEQSEQ (sal_Int32)& rDates)567 double SAL_CALL AnalysisAddIn::getXnpv(
568 	double fRate, const SEQSEQ( double )& rValues, const SEQSEQ( sal_Int32 )& rDates ) THROWDEF_RTE_IAE
569 {
570     ScaDoubleList aValList;
571     ScaDoubleList aDateList;
572 
573 	aValList.Append( rValues );
574 	aDateList.Append( rDates );
575 
576 	sal_Int32			nNum = aValList.Count();
577 
578 	if( nNum != sal_Int32( aDateList.Count() ) || nNum < 2 )
579 		THROW_IAE;
580 
581 	double				fRet = 0.0;
582 	double				fNull = *aDateList.Get( 0 );
583 	fRate++;
584 
585 	for( sal_Int32 i = 0 ; i < nNum ; i++ )
586 		fRet += *aValList.Get( i ) / ( pow( fRate, ( *aDateList.Get( i ) - fNull ) / 365.0 ) );
587 
588     RETURN_FINITE( fRet );
589 }
590 
591 
getIntrate(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,double fInvest,double fRedemp,const ANY & rOB)592 double SAL_CALL AnalysisAddIn::getIntrate( constREFXPS& xOpt,
593 	sal_Int32 nSettle, sal_Int32 nMat, double fInvest, double fRedemp, const ANY& rOB ) THROWDEF_RTE_IAE
594 {
595 	if( fInvest <= 0.0 || fRedemp <= 0.0 || nSettle >= nMat )
596 		THROW_IAE;
597 
598     double fRet = ( ( fRedemp / fInvest ) - 1.0 ) / GetYearDiff( GetNullDate( xOpt ), nSettle, nMat, getDateMode( xOpt, rOB ) );
599     RETURN_FINITE( fRet );
600 }
601 
602 
getCoupncd(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)603 double SAL_CALL AnalysisAddIn::getCoupncd( constREFXPS& xOpt,
604 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
605 {
606     double fRet = GetCoupncd( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
607     RETURN_FINITE( fRet );
608 }
609 
610 
getCoupdays(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)611 double SAL_CALL AnalysisAddIn::getCoupdays( constREFXPS& xOpt,
612 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
613 {
614     double fRet = GetCoupdays( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
615     RETURN_FINITE( fRet );
616 }
617 
618 
getCoupdaysnc(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)619 double SAL_CALL AnalysisAddIn::getCoupdaysnc( constREFXPS& xOpt,
620 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
621 {
622     double fRet = GetCoupdaysnc( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
623     RETURN_FINITE( fRet );
624 }
625 
626 
getCoupdaybs(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)627 double SAL_CALL AnalysisAddIn::getCoupdaybs( constREFXPS& xOpt,
628 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
629 {
630     double fRet = GetCoupdaybs( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
631     RETURN_FINITE( fRet );
632 }
633 
634 
getCouppcd(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)635 double SAL_CALL AnalysisAddIn::getCouppcd( constREFXPS& xOpt,
636 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
637 {
638     double fRet = GetCouppcd( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
639     RETURN_FINITE( fRet );
640 }
641 
642 
getCoupnum(constREFXPS & xOpt,sal_Int32 nSettle,sal_Int32 nMat,sal_Int32 nFreq,const ANY & rOB)643 double SAL_CALL AnalysisAddIn::getCoupnum( constREFXPS& xOpt,
644 	sal_Int32 nSettle, sal_Int32 nMat, sal_Int32 nFreq, const ANY& rOB ) THROWDEF_RTE_IAE
645 {
646     double fRet = GetCoupnum( GetNullDate( xOpt ), nSettle, nMat, nFreq, getDateMode( xOpt, rOB ) );
647     RETURN_FINITE( fRet );
648 }
649 
650 
getFvschedule(double fPrinc,const SEQSEQ (double)& rSchedule)651 double SAL_CALL AnalysisAddIn::getFvschedule( double fPrinc, const SEQSEQ( double )& rSchedule ) THROWDEF_RTE_IAE
652 {
653     ScaDoubleList aSchedList;
654 
655 	aSchedList.Append( rSchedule );
656 
657 	for( const double* p = aSchedList.First() ; p ; p = aSchedList.Next() )
658 		fPrinc *= 1.0 + *p;
659 
660     RETURN_FINITE( fPrinc );
661 }
662 
663 
664