/* TA-LIB Copyright (c) 1999-2007, Mario Fortier * All rights reserved. * * Redistribution and use in source and binary forms, with or * without modification, are permitted provided that the following * conditions are met: * * - Redistributions of source code must retain the above copyright * notice, this list of conditions and the following disclaimer. * * - Redistributions in binary form must reproduce the above copyright * notice, this list of conditions and the following disclaimer in * the documentation and/or other materials provided with the * distribution. * * - Neither name of author nor the names of its contributors * may be used to endorse or promote products derived from this * software without specific prior written permission. * * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE * REGENTS OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, * WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE * OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, * EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. */ /* List of contributors: * * Initial Name/description * ------------------------------------------------------------------- * MF Mario Fortier * * * Change history: * * MMDDYY BY Description * ------------------------------------------------------------------- * 120802 MF Template creation. * 052603 MF Adapt code to compile with .NET Managed C++ * */ /**** START GENCODE SECTION 1 - DO NOT DELETE THIS LINE ****/ /* All code within this section is automatically * generated by gen_code. Any modification will be lost * next time gen_code is run. */ /* Generated */ /* Generated */ #if defined( _MANAGED ) /* Generated */ #include "TA-Lib-Core.h" /* Generated */ #define TA_INTERNAL_ERROR(Id) (RetCode::InternalError) /* Generated */ namespace TicTacTec { namespace TA { namespace Library { /* Generated */ #elif defined( _JAVA ) /* Generated */ #include "ta_defs.h" /* Generated */ #include "ta_java_defs.h" /* Generated */ #define TA_INTERNAL_ERROR(Id) (RetCode.InternalError) /* Generated */ #else /* Generated */ #include /* Generated */ #include /* Generated */ #include "ta_func.h" /* Generated */ #endif /* Generated */ /* Generated */ #ifndef TA_UTILITY_H /* Generated */ #include "ta_utility.h" /* Generated */ #endif /* Generated */ /* Generated */ #ifndef TA_MEMORY_H /* Generated */ #include "ta_memory.h" /* Generated */ #endif /* Generated */ /* Generated */ #define TA_PREFIX(x) TA_##x /* Generated */ #define INPUT_TYPE double /* Generated */ /* Generated */ #if defined( _MANAGED ) /* Generated */ int Core::HtTrendModeLookback( void ) /* Generated */ /* Generated */ #elif defined( _JAVA ) /* Generated */ public int htTrendModeLookback( ) /* Generated */ /* Generated */ #else /* Generated */ int TA_HT_TRENDMODE_Lookback( void ) /* Generated */ /* Generated */ #endif /**** END GENCODE SECTION 1 - DO NOT DELETE THIS LINE ****/ { /* insert local variable here */ /**** START GENCODE SECTION 2 - DO NOT DELETE THIS LINE ****/ /* Generated */ /* No parameters to validate. */ /**** END GENCODE SECTION 2 - DO NOT DELETE THIS LINE ****/ /* insert lookback code here. */ /* 31 input are skip * +32 output are skip to account for misc lookback * --- * 63 Total Lookback * * 31 is for being compatible with Tradestation. * See TA_MAMA_Lookback for an explanation of the "32". */ return 63 + TA_GLOBALS_UNSTABLE_PERIOD(TA_FUNC_UNST_HT_TRENDMODE,HtTrendMode); } /**** START GENCODE SECTION 3 - DO NOT DELETE THIS LINE ****/ /* * TA_HT_TRENDMODE - Hilbert Transform - Trend vs Cycle Mode * * Input = double * Output = int * */ /* Generated */ /* Generated */ #if defined( _MANAGED ) && defined( USE_SUBARRAY ) /* Generated */ enum class Core::RetCode Core::HtTrendMode( int startIdx, /* Generated */ int endIdx, /* Generated */ SubArray^ inReal, /* Generated */ [Out]int% outBegIdx, /* Generated */ [Out]int% outNBElement, /* Generated */ cli::array^ outInteger ) /* Generated */ #elif defined( _MANAGED ) /* Generated */ enum class Core::RetCode Core::HtTrendMode( int startIdx, /* Generated */ int endIdx, /* Generated */ cli::array^ inReal, /* Generated */ [Out]int% outBegIdx, /* Generated */ [Out]int% outNBElement, /* Generated */ cli::array^ outInteger ) /* Generated */ #elif defined( _JAVA ) /* Generated */ public RetCode htTrendMode( int startIdx, /* Generated */ int endIdx, /* Generated */ double inReal[], /* Generated */ MInteger outBegIdx, /* Generated */ MInteger outNBElement, /* Generated */ int outInteger[] ) /* Generated */ #else /* Generated */ TA_RetCode TA_HT_TRENDMODE( int startIdx, /* Generated */ int endIdx, /* Generated */ const double inReal[], /* Generated */ int *outBegIdx, /* Generated */ int *outNBElement, /* Generated */ int outInteger[] ) /* Generated */ #endif /**** END GENCODE SECTION 3 - DO NOT DELETE THIS LINE ****/ { /* insert local variable here */ int outIdx, i; int lookbackTotal, today; double tempReal, tempReal2; double adjustedPrevPeriod, period; /* Variable used for the price smoother (a weighted moving average). */ int trailingWMAIdx; double periodWMASum, periodWMASub, trailingWMAValue; double smoothedValue; /* Variable to keep track of the last 3 ITrend */ double iTrend1, iTrend2, iTrend3; /* Variables used for the Hilbert Transormation */ CONSTANT_DOUBLE(a) = 0.0962; CONSTANT_DOUBLE(b) = 0.5769; double hilbertTempReal; int hilbertIdx; HILBERT_VARIABLES( detrender ); HILBERT_VARIABLES( Q1 ); HILBERT_VARIABLES( jI ); HILBERT_VARIABLES( jQ ); double Q2, I2, prevQ2, prevI2, Re, Im; double I1ForOddPrev2, I1ForOddPrev3; double I1ForEvenPrev2, I1ForEvenPrev3; double rad2Deg, deg2Rad, constDeg2RadBy360; double todayValue, smoothPeriod; /* Variable used to keep track of the previous * smooth price. In the case of this algorithm, * we will never need more than 50 values. */ #define SMOOTH_PRICE_SIZE 50 CIRCBUF_PROLOG(smoothPrice,double,SMOOTH_PRICE_SIZE); int idx; /* Variable used to calculate the dominant cycle phase */ int DCPeriodInt; double DCPhase, DCPeriod, imagPart, realPart; /* Variable used to calculate the trend mode */ int daysInTrend, trend; double prevDCPhase, trendline; double prevSine, prevLeadSine, sine, leadSine; /**** START GENCODE SECTION 4 - DO NOT DELETE THIS LINE ****/ /* Generated */ /* Generated */ #ifndef TA_FUNC_NO_RANGE_CHECK /* Generated */ /* Generated */ /* Validate the requested output range. */ /* Generated */ if( startIdx < 0 ) /* Generated */ return ENUM_VALUE(RetCode,TA_OUT_OF_RANGE_START_INDEX,OutOfRangeStartIndex); /* Generated */ if( (endIdx < 0) || (endIdx < startIdx)) /* Generated */ return ENUM_VALUE(RetCode,TA_OUT_OF_RANGE_END_INDEX,OutOfRangeEndIndex); /* Generated */ /* Generated */ #if !defined(_JAVA) /* Generated */ if( !inReal ) return ENUM_VALUE(RetCode,TA_BAD_PARAM,BadParam); /* Generated */ #endif /* !defined(_JAVA)*/ /* Generated */ #if !defined(_JAVA) /* Generated */ if( !outInteger ) /* Generated */ return ENUM_VALUE(RetCode,TA_BAD_PARAM,BadParam); /* Generated */ /* Generated */ #endif /* !defined(_JAVA) */ /* Generated */ #endif /* TA_FUNC_NO_RANGE_CHECK */ /* Generated */ /**** END GENCODE SECTION 4 - DO NOT DELETE THIS LINE ****/ /* Insert TA function code here. */ CIRCBUF_INIT_LOCAL_ONLY(smoothPrice,double); iTrend1 = iTrend2 = iTrend3 = 0.0; daysInTrend = 0; prevDCPhase = DCPhase = 0.0; prevSine = sine = 0.0; prevLeadSine = leadSine = 0.0; /* The following could be replaced by constant eventually. */ tempReal = std_atan(1); rad2Deg = 45.0/tempReal; deg2Rad = 1.0/rad2Deg; constDeg2RadBy360 = tempReal*8.0; /* Identify the minimum number of price bar needed * to calculate at least one output. */ lookbackTotal = 63 + TA_GLOBALS_UNSTABLE_PERIOD(TA_FUNC_UNST_HT_TRENDMODE,HtTrendMode); /* Move up the start index if there is not * enough initial data. */ if( startIdx < lookbackTotal ) startIdx = lookbackTotal; /* Make sure there is still something to evaluate. */ if( startIdx > endIdx ) { VALUE_HANDLE_DEREF_TO_ZERO(outBegIdx); VALUE_HANDLE_DEREF_TO_ZERO(outNBElement); return ENUM_VALUE(RetCode,TA_SUCCESS,Success); } VALUE_HANDLE_DEREF(outBegIdx) = startIdx; /* Initialize the price smoother, which is simply a weighted * moving average of the price. * To understand this algorithm, I strongly suggest to understand * first how TA_WMA is done. */ trailingWMAIdx = startIdx - lookbackTotal; today = trailingWMAIdx; /* Initialization is same as WMA, except loop is unrolled * for speed optimization. */ tempReal = inReal[today++]; periodWMASub = tempReal; periodWMASum = tempReal; tempReal = inReal[today++]; periodWMASub += tempReal; periodWMASum += tempReal*2.0; tempReal = inReal[today++]; periodWMASub += tempReal; periodWMASum += tempReal*3.0; trailingWMAValue = 0.0; /* Subsequent WMA value are evaluated by using * the DO_PRICE_WMA macro. */ #define DO_PRICE_WMA(varNewPrice,varToStoreSmoothedValue) { \ periodWMASub += varNewPrice; \ periodWMASub -= trailingWMAValue; \ periodWMASum += varNewPrice*4.0; \ trailingWMAValue = inReal[trailingWMAIdx++]; \ varToStoreSmoothedValue = periodWMASum*0.1; \ periodWMASum -= periodWMASub; \ } i = 34; do { tempReal = inReal[today++]; DO_PRICE_WMA(tempReal,smoothedValue); } while( --i != 0); /* Initialize the circular buffers used by the hilbert * transform logic. * A buffer is used for odd day and another for even days. * This minimize the number of memory access and floating point * operations needed (note also that by using static circular buffer, * no large dynamic memory allocation is needed for storing * intermediate calculation!). */ hilbertIdx = 0; INIT_HILBERT_VARIABLES(detrender); INIT_HILBERT_VARIABLES(Q1); INIT_HILBERT_VARIABLES(jI); INIT_HILBERT_VARIABLES(jQ); period = 0.0; outIdx = 0; prevI2 = prevQ2 = 0.0; Re = Im = 0.0; I1ForOddPrev3 = I1ForEvenPrev3 = 0.0; I1ForOddPrev2 = I1ForEvenPrev2 = 0.0; smoothPeriod = 0.0; for( i=0; i < SMOOTH_PRICE_SIZE; i++ ) smoothPrice[i] = 0.0; /* The code is speed optimized and is most likely very * hard to follow if you do not already know well the * original algorithm. * To understadn better, it is strongly suggested to look * first at the Excel implementation in "test_MAMA.xls" included * in this package. */ DCPhase = 0.0; while( today <= endIdx ) { adjustedPrevPeriod = (0.075*period)+0.54; todayValue = inReal[today]; DO_PRICE_WMA(todayValue,smoothedValue); /* Remember the smoothedValue into the smoothPrice * circular buffer. */ smoothPrice[smoothPrice_Idx] = smoothedValue; if( (today%2) == 0 ) { /* Do the Hilbert Transforms for even price bar */ DO_HILBERT_EVEN(detrender,smoothedValue); DO_HILBERT_EVEN(Q1,detrender); DO_HILBERT_EVEN(jI,I1ForEvenPrev3); DO_HILBERT_EVEN(jQ,Q1); if( ++hilbertIdx == 3 ) hilbertIdx = 0; Q2 = (0.2*(Q1 + jI)) + (0.8*prevQ2); I2 = (0.2*(I1ForEvenPrev3 - jQ)) + (0.8*prevI2); /* The variable I1 is the detrender delayed for * 3 price bars. * * Save the current detrender value for being * used by the "odd" logic later. */ I1ForOddPrev3 = I1ForOddPrev2; I1ForOddPrev2 = detrender; } else { /* Do the Hilbert Transforms for odd price bar */ DO_HILBERT_ODD(detrender,smoothedValue); DO_HILBERT_ODD(Q1,detrender); DO_HILBERT_ODD(jI,I1ForOddPrev3); DO_HILBERT_ODD(jQ,Q1); Q2 = (0.2*(Q1 + jI)) + (0.8*prevQ2); I2 = (0.2*(I1ForOddPrev3 - jQ)) + (0.8*prevI2); /* The varaiable I1 is the detrender delayed for * 3 price bars. * * Save the current detrender value for being * used by the "even" logic later. */ I1ForEvenPrev3 = I1ForEvenPrev2; I1ForEvenPrev2 = detrender; } /* Adjust the period for next price bar */ Re = (0.2*((I2*prevI2)+(Q2*prevQ2)))+(0.8*Re); Im = (0.2*((I2*prevQ2)-(Q2*prevI2)))+(0.8*Im); prevQ2 = Q2; prevI2 = I2; tempReal = period; if( (Im != 0.0) && (Re != 0.0) ) period = 360.0 / (std_atan(Im/Re)*rad2Deg); tempReal2 = 1.5*tempReal; if( period > tempReal2) period = tempReal2; tempReal2 = 0.67*tempReal; if( period < tempReal2 ) period = tempReal2; if( period < 6 ) period = 6; else if( period > 50 ) period = 50; period = (0.2*period) + (0.8 * tempReal); smoothPeriod = (0.33*period)+(0.67*smoothPeriod); /* Compute Dominant Cycle Phase */ prevDCPhase = DCPhase; DCPeriod = smoothPeriod+0.5; DCPeriodInt = (int)DCPeriod; realPart = 0.0; imagPart = 0.0; /* idx is used to iterate for up to 50 of the last * value of smoothPrice. */ idx = smoothPrice_Idx; for( i=0; i < DCPeriodInt; i++ ) { tempReal = ((double)i*constDeg2RadBy360)/(double)DCPeriodInt; tempReal2 = smoothPrice[idx]; realPart += std_sin(tempReal)*tempReal2; imagPart += std_cos(tempReal)*tempReal2; if( idx == 0 ) idx = SMOOTH_PRICE_SIZE-1; else idx--; } tempReal = std_fabs(imagPart); if( tempReal > 0.0 ) DCPhase = std_atan(realPart/imagPart)*rad2Deg; else if( tempReal <= 0.01 ) { if( realPart < 0.0 ) DCPhase -= 90.0; else if( realPart > 0.0 ) DCPhase += 90.0; } DCPhase += 90.0; /* Compensate for one bar lag of the weighted moving average */ DCPhase += 360.0 / smoothPeriod; if( imagPart < 0.0 ) DCPhase += 180.0; if( DCPhase > 315.0 ) DCPhase -= 360.0; prevSine = sine; prevLeadSine = leadSine; sine = std_sin(DCPhase*deg2Rad); leadSine = std_sin((DCPhase+45)*deg2Rad); /* Compute Trendline */ DCPeriod = smoothPeriod+0.5; DCPeriodInt = (int)DCPeriod; /* idx is used to iterate for up to 50 of the last * value of smoothPrice. */ idx = today; tempReal = 0.0; for( i=0; i < DCPeriodInt; i++ ) tempReal += inReal[idx--]; if( DCPeriodInt > 0 ) tempReal = tempReal/(double)DCPeriodInt; trendline = (4.0*tempReal + 3.0*iTrend1 + 2.0*iTrend2 + iTrend3) / 10.0; iTrend3 = iTrend2; iTrend2 = iTrend1; iTrend1 = tempReal; /* Compute the trend Mode , and assume trend by default */ trend = 1; /* Measure days in trend from last crossing of the SineWave Indicator lines */ if( ((sine > leadSine) && (prevSine <= prevLeadSine)) || ((sine < leadSine) && (prevSine >= prevLeadSine)) ) { daysInTrend = 0; trend = 0; } daysInTrend++; if( daysInTrend < (0.5*smoothPeriod) ) trend = 0; tempReal = DCPhase - prevDCPhase; if( (smoothPeriod != 0.0) && ((tempReal > (0.67*360.0/smoothPeriod)) && (tempReal < (1.5*360.0/smoothPeriod))) ) { trend = 0; } tempReal = smoothPrice[smoothPrice_Idx]; if( (trendline != 0.0) && (std_fabs( (tempReal - trendline)/trendline ) >= 0.015) ) trend = 1; if( today >= startIdx ) { outInteger[outIdx++] = trend; } /* Ooof... let's do the next price bar now! */ CIRCBUF_NEXT(smoothPrice); today++; } VALUE_HANDLE_DEREF(outNBElement) = outIdx; return ENUM_VALUE(RetCode,TA_SUCCESS,Success); } /**** START GENCODE SECTION 5 - DO NOT DELETE THIS LINE ****/ /* Generated */ /* Generated */ #define USE_SINGLE_PRECISION_INPUT /* Generated */ #if !defined( _MANAGED ) && !defined( _JAVA ) /* Generated */ #undef TA_PREFIX /* Generated */ #define TA_PREFIX(x) TA_S_##x /* Generated */ #endif /* Generated */ #undef INPUT_TYPE /* Generated */ #define INPUT_TYPE float /* Generated */ #if defined( _MANAGED ) /* Generated */ enum class Core::RetCode Core::HtTrendMode( int startIdx, /* Generated */ int endIdx, /* Generated */ cli::array^ inReal, /* Generated */ [Out]int% outBegIdx, /* Generated */ [Out]int% outNBElement, /* Generated */ cli::array^ outInteger ) /* Generated */ #elif defined( _JAVA ) /* Generated */ public RetCode htTrendMode( int startIdx, /* Generated */ int endIdx, /* Generated */ float inReal[], /* Generated */ MInteger outBegIdx, /* Generated */ MInteger outNBElement, /* Generated */ int outInteger[] ) /* Generated */ #else /* Generated */ TA_RetCode TA_S_HT_TRENDMODE( int startIdx, /* Generated */ int endIdx, /* Generated */ const float inReal[], /* Generated */ int *outBegIdx, /* Generated */ int *outNBElement, /* Generated */ int outInteger[] ) /* Generated */ #endif /* Generated */ { /* Generated */ int outIdx, i; /* Generated */ int lookbackTotal, today; /* Generated */ double tempReal, tempReal2; /* Generated */ double adjustedPrevPeriod, period; /* Generated */ int trailingWMAIdx; /* Generated */ double periodWMASum, periodWMASub, trailingWMAValue; /* Generated */ double smoothedValue; /* Generated */ double iTrend1, iTrend2, iTrend3; /* Generated */ CONSTANT_DOUBLE(a) = 0.0962; /* Generated */ CONSTANT_DOUBLE(b) = 0.5769; /* Generated */ double hilbertTempReal; /* Generated */ int hilbertIdx; /* Generated */ HILBERT_VARIABLES( detrender ); /* Generated */ HILBERT_VARIABLES( Q1 ); /* Generated */ HILBERT_VARIABLES( jI ); /* Generated */ HILBERT_VARIABLES( jQ ); /* Generated */ double Q2, I2, prevQ2, prevI2, Re, Im; /* Generated */ double I1ForOddPrev2, I1ForOddPrev3; /* Generated */ double I1ForEvenPrev2, I1ForEvenPrev3; /* Generated */ double rad2Deg, deg2Rad, constDeg2RadBy360; /* Generated */ double todayValue, smoothPeriod; /* Generated */ #define SMOOTH_PRICE_SIZE 50 /* Generated */ CIRCBUF_PROLOG(smoothPrice,double,SMOOTH_PRICE_SIZE); /* Generated */ int idx; /* Generated */ int DCPeriodInt; /* Generated */ double DCPhase, DCPeriod, imagPart, realPart; /* Generated */ int daysInTrend, trend; /* Generated */ double prevDCPhase, trendline; /* Generated */ double prevSine, prevLeadSine, sine, leadSine; /* Generated */ #ifndef TA_FUNC_NO_RANGE_CHECK /* Generated */ if( startIdx < 0 ) /* Generated */ return ENUM_VALUE(RetCode,TA_OUT_OF_RANGE_START_INDEX,OutOfRangeStartIndex); /* Generated */ if( (endIdx < 0) || (endIdx < startIdx)) /* Generated */ return ENUM_VALUE(RetCode,TA_OUT_OF_RANGE_END_INDEX,OutOfRangeEndIndex); /* Generated */ #if !defined(_JAVA) /* Generated */ if( !inReal ) return ENUM_VALUE(RetCode,TA_BAD_PARAM,BadParam); /* Generated */ #endif /* Generated */ #if !defined(_JAVA) /* Generated */ if( !outInteger ) /* Generated */ return ENUM_VALUE(RetCode,TA_BAD_PARAM,BadParam); /* Generated */ #endif /* Generated */ #endif /* Generated */ CIRCBUF_INIT_LOCAL_ONLY(smoothPrice,double); /* Generated */ iTrend1 = iTrend2 = iTrend3 = 0.0; /* Generated */ daysInTrend = 0; /* Generated */ prevDCPhase = DCPhase = 0.0; /* Generated */ prevSine = sine = 0.0; /* Generated */ prevLeadSine = leadSine = 0.0; /* Generated */ tempReal = std_atan(1); /* Generated */ rad2Deg = 45.0/tempReal; /* Generated */ deg2Rad = 1.0/rad2Deg; /* Generated */ constDeg2RadBy360 = tempReal*8.0; /* Generated */ lookbackTotal = 63 + TA_GLOBALS_UNSTABLE_PERIOD(TA_FUNC_UNST_HT_TRENDMODE,HtTrendMode); /* Generated */ if( startIdx < lookbackTotal ) /* Generated */ startIdx = lookbackTotal; /* Generated */ if( startIdx > endIdx ) /* Generated */ { /* Generated */ VALUE_HANDLE_DEREF_TO_ZERO(outBegIdx); /* Generated */ VALUE_HANDLE_DEREF_TO_ZERO(outNBElement); /* Generated */ return ENUM_VALUE(RetCode,TA_SUCCESS,Success); /* Generated */ } /* Generated */ VALUE_HANDLE_DEREF(outBegIdx) = startIdx; /* Generated */ trailingWMAIdx = startIdx - lookbackTotal; /* Generated */ today = trailingWMAIdx; /* Generated */ tempReal = inReal[today++]; /* Generated */ periodWMASub = tempReal; /* Generated */ periodWMASum = tempReal; /* Generated */ tempReal = inReal[today++]; /* Generated */ periodWMASub += tempReal; /* Generated */ periodWMASum += tempReal*2.0; /* Generated */ tempReal = inReal[today++]; /* Generated */ periodWMASub += tempReal; /* Generated */ periodWMASum += tempReal*3.0; /* Generated */ trailingWMAValue = 0.0; /* Generated */ #define DO_PRICE_WMA(varNewPrice,varToStoreSmoothedValue) { \ /* Generated */ periodWMASub += varNewPrice; \ /* Generated */ periodWMASub -= trailingWMAValue; \ /* Generated */ periodWMASum += varNewPrice*4.0; \ /* Generated */ trailingWMAValue = inReal[trailingWMAIdx++]; \ /* Generated */ varToStoreSmoothedValue = periodWMASum*0.1; \ /* Generated */ periodWMASum -= periodWMASub; \ /* Generated */ } /* Generated */ i = 34; /* Generated */ do /* Generated */ { /* Generated */ tempReal = inReal[today++]; /* Generated */ DO_PRICE_WMA(tempReal,smoothedValue); /* Generated */ } while( --i != 0); /* Generated */ hilbertIdx = 0; /* Generated */ INIT_HILBERT_VARIABLES(detrender); /* Generated */ INIT_HILBERT_VARIABLES(Q1); /* Generated */ INIT_HILBERT_VARIABLES(jI); /* Generated */ INIT_HILBERT_VARIABLES(jQ); /* Generated */ period = 0.0; /* Generated */ outIdx = 0; /* Generated */ prevI2 = prevQ2 = 0.0; /* Generated */ Re = Im = 0.0; /* Generated */ I1ForOddPrev3 = I1ForEvenPrev3 = 0.0; /* Generated */ I1ForOddPrev2 = I1ForEvenPrev2 = 0.0; /* Generated */ smoothPeriod = 0.0; /* Generated */ for( i=0; i < SMOOTH_PRICE_SIZE; i++ ) /* Generated */ smoothPrice[i] = 0.0; /* Generated */ DCPhase = 0.0; /* Generated */ while( today <= endIdx ) /* Generated */ { /* Generated */ adjustedPrevPeriod = (0.075*period)+0.54; /* Generated */ todayValue = inReal[today]; /* Generated */ DO_PRICE_WMA(todayValue,smoothedValue); /* Generated */ smoothPrice[smoothPrice_Idx] = smoothedValue; /* Generated */ if( (today%2) == 0 ) /* Generated */ { /* Generated */ DO_HILBERT_EVEN(detrender,smoothedValue); /* Generated */ DO_HILBERT_EVEN(Q1,detrender); /* Generated */ DO_HILBERT_EVEN(jI,I1ForEvenPrev3); /* Generated */ DO_HILBERT_EVEN(jQ,Q1); /* Generated */ if( ++hilbertIdx == 3 ) /* Generated */ hilbertIdx = 0; /* Generated */ Q2 = (0.2*(Q1 + jI)) + (0.8*prevQ2); /* Generated */ I2 = (0.2*(I1ForEvenPrev3 - jQ)) + (0.8*prevI2); /* Generated */ I1ForOddPrev3 = I1ForOddPrev2; /* Generated */ I1ForOddPrev2 = detrender; /* Generated */ } /* Generated */ else /* Generated */ { /* Generated */ DO_HILBERT_ODD(detrender,smoothedValue); /* Generated */ DO_HILBERT_ODD(Q1,detrender); /* Generated */ DO_HILBERT_ODD(jI,I1ForOddPrev3); /* Generated */ DO_HILBERT_ODD(jQ,Q1); /* Generated */ Q2 = (0.2*(Q1 + jI)) + (0.8*prevQ2); /* Generated */ I2 = (0.2*(I1ForOddPrev3 - jQ)) + (0.8*prevI2); /* Generated */ I1ForEvenPrev3 = I1ForEvenPrev2; /* Generated */ I1ForEvenPrev2 = detrender; /* Generated */ } /* Generated */ Re = (0.2*((I2*prevI2)+(Q2*prevQ2)))+(0.8*Re); /* Generated */ Im = (0.2*((I2*prevQ2)-(Q2*prevI2)))+(0.8*Im); /* Generated */ prevQ2 = Q2; /* Generated */ prevI2 = I2; /* Generated */ tempReal = period; /* Generated */ if( (Im != 0.0) && (Re != 0.0) ) /* Generated */ period = 360.0 / (std_atan(Im/Re)*rad2Deg); /* Generated */ tempReal2 = 1.5*tempReal; /* Generated */ if( period > tempReal2) /* Generated */ period = tempReal2; /* Generated */ tempReal2 = 0.67*tempReal; /* Generated */ if( period < tempReal2 ) /* Generated */ period = tempReal2; /* Generated */ if( period < 6 ) /* Generated */ period = 6; /* Generated */ else if( period > 50 ) /* Generated */ period = 50; /* Generated */ period = (0.2*period) + (0.8 * tempReal); /* Generated */ smoothPeriod = (0.33*period)+(0.67*smoothPeriod); /* Generated */ prevDCPhase = DCPhase; /* Generated */ DCPeriod = smoothPeriod+0.5; /* Generated */ DCPeriodInt = (int)DCPeriod; /* Generated */ realPart = 0.0; /* Generated */ imagPart = 0.0; /* Generated */ idx = smoothPrice_Idx; /* Generated */ for( i=0; i < DCPeriodInt; i++ ) /* Generated */ { /* Generated */ tempReal = ((double)i*constDeg2RadBy360)/(double)DCPeriodInt; /* Generated */ tempReal2 = smoothPrice[idx]; /* Generated */ realPart += std_sin(tempReal)*tempReal2; /* Generated */ imagPart += std_cos(tempReal)*tempReal2; /* Generated */ if( idx == 0 ) /* Generated */ idx = SMOOTH_PRICE_SIZE-1; /* Generated */ else /* Generated */ idx--; /* Generated */ } /* Generated */ tempReal = std_fabs(imagPart); /* Generated */ if( tempReal > 0.0 ) /* Generated */ DCPhase = std_atan(realPart/imagPart)*rad2Deg; /* Generated */ else if( tempReal <= 0.01 ) /* Generated */ { /* Generated */ if( realPart < 0.0 ) /* Generated */ DCPhase -= 90.0; /* Generated */ else if( realPart > 0.0 ) /* Generated */ DCPhase += 90.0; /* Generated */ } /* Generated */ DCPhase += 90.0; /* Generated */ DCPhase += 360.0 / smoothPeriod; /* Generated */ if( imagPart < 0.0 ) /* Generated */ DCPhase += 180.0; /* Generated */ if( DCPhase > 315.0 ) /* Generated */ DCPhase -= 360.0; /* Generated */ prevSine = sine; /* Generated */ prevLeadSine = leadSine; /* Generated */ sine = std_sin(DCPhase*deg2Rad); /* Generated */ leadSine = std_sin((DCPhase+45)*deg2Rad); /* Generated */ DCPeriod = smoothPeriod+0.5; /* Generated */ DCPeriodInt = (int)DCPeriod; /* Generated */ idx = today; /* Generated */ tempReal = 0.0; /* Generated */ for( i=0; i < DCPeriodInt; i++ ) /* Generated */ tempReal += inReal[idx--]; /* Generated */ if( DCPeriodInt > 0 ) /* Generated */ tempReal = tempReal/(double)DCPeriodInt; /* Generated */ trendline = (4.0*tempReal + 3.0*iTrend1 + 2.0*iTrend2 + iTrend3) / 10.0; /* Generated */ iTrend3 = iTrend2; /* Generated */ iTrend2 = iTrend1; /* Generated */ iTrend1 = tempReal; /* Generated */ trend = 1; /* Generated */ if( ((sine > leadSine) && (prevSine <= prevLeadSine)) || /* Generated */ ((sine < leadSine) && (prevSine >= prevLeadSine)) ) /* Generated */ { /* Generated */ daysInTrend = 0; /* Generated */ trend = 0; /* Generated */ } /* Generated */ daysInTrend++; /* Generated */ if( daysInTrend < (0.5*smoothPeriod) ) /* Generated */ trend = 0; /* Generated */ tempReal = DCPhase - prevDCPhase; /* Generated */ if( (smoothPeriod != 0.0) && /* Generated */ ((tempReal > (0.67*360.0/smoothPeriod)) && (tempReal < (1.5*360.0/smoothPeriod))) ) /* Generated */ { /* Generated */ trend = 0; /* Generated */ } /* Generated */ tempReal = smoothPrice[smoothPrice_Idx]; /* Generated */ if( (trendline != 0.0) && (std_fabs( (tempReal - trendline)/trendline ) >= 0.015) ) /* Generated */ trend = 1; /* Generated */ if( today >= startIdx ) /* Generated */ { /* Generated */ outInteger[outIdx++] = trend; /* Generated */ } /* Generated */ CIRCBUF_NEXT(smoothPrice); /* Generated */ today++; /* Generated */ } /* Generated */ VALUE_HANDLE_DEREF(outNBElement) = outIdx; /* Generated */ return ENUM_VALUE(RetCode,TA_SUCCESS,Success); /* Generated */ } /* Generated */ /* Generated */ #if defined( _MANAGED ) /* Generated */ }}} // Close namespace TicTacTec.TA.Lib /* Generated */ #endif /**** END GENCODE SECTION 5 - DO NOT DELETE THIS LINE ****/