![]() |
pl-nk v0.4.5
Plonk|Plink|Plank are a set of cross-platform C/C++ frameworks for audio software development
|
00001 /* 00002 ------------------------------------------------------------------------------- 00003 This file is part of the Plink, Plonk, Plank libraries 00004 by Martin Robinson 00005 00006 http://code.google.com/p/pl-nk/ 00007 00008 Copyright University of the West of England, Bristol 2011-14 00009 All rights reserved. 00010 00011 Redistribution and use in source and binary forms, with or without 00012 modification, are permitted provided that the following conditions are met: 00013 00014 * Redistributions of source code must retain the above copyright 00015 notice, this list of conditions and the following disclaimer. 00016 * Redistributions in binary form must reproduce the above copyright 00017 notice, this list of conditions and the following disclaimer in the 00018 documentation and/or other materials provided with the distribution. 00019 * Neither the name of University of the West of England, Bristol nor 00020 the names of its contributors may be used to endorse or promote products 00021 derived from this software without specific prior written permission. 00022 00023 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND 00024 ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED 00025 WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE 00026 DISCLAIMED. IN NO EVENT SHALL UNIVERSITY OF THE WEST OF ENGLAND, BRISTOL BE 00027 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 00028 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE 00029 GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 00030 HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 00031 LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT 00032 OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00033 00034 This software makes use of third party libraries. For more information see: 00035 doc/license.txt included in the distribution. 00036 ------------------------------------------------------------------------------- 00037 */ 00038 00039 #ifndef PLONK_ZMULCHANNEL_H 00040 #define PLONK_ZMULCHANNEL_H 00041 00042 #include "../channel/plonk_ChannelInternalCore.h" 00043 #include "../plonk_GraphForwardDeclarations.h" 00044 00046 template<class SampleType> 00047 class ZMulChannelInternal 00048 : public ChannelInternal<SampleType, ChannelInternalCore::Data> 00049 { 00050 public: 00051 typedef ChannelInternalCore::Data Data; 00052 typedef ChannelBase<SampleType> ChannelType; 00053 typedef ObjectArray<ChannelType> ChannelArrayType; 00054 typedef ZMulChannelInternal<SampleType> ZMulInternal; 00055 typedef ChannelInternal<SampleType,Data> Internal; 00056 typedef ChannelInternalBase<SampleType> InternalBase; 00057 typedef UnitBase<SampleType> UnitType; 00058 typedef InputDictionary Inputs; 00059 typedef NumericalArray<SampleType> Buffer; 00060 00061 ZMulChannelInternal (Inputs const& inputs, 00062 Data const& data, 00063 BlockSize const& blockSize, 00064 SampleRate const& sampleRate) throw() 00065 : Internal (inputs, data, blockSize, sampleRate) 00066 { 00067 } 00068 00069 Text getName() const throw() 00070 { 00071 return "ZMul"; 00072 } 00073 00074 IntArray getInputKeys() const throw() 00075 { 00076 const IntArray keys (IOKey::LeftOperand, IOKey::RightOperand); 00077 return keys; 00078 } 00079 00080 InternalBase* getChannel (const int index) throw() 00081 { 00082 const Inputs channelInputs = this->getInputs().getChannel (index); 00083 return new ZMulInternal (channelInputs, 00084 this->getState(), 00085 this->getBlockSize(), 00086 this->getSampleRate()); 00087 } 00088 00089 void initChannel (const int channel) throw() 00090 { 00091 const UnitType& leftOperand = this->getInputAsUnit (IOKey::LeftOperand); 00092 this->setBlockSize (leftOperand.getBlockSize (channel)); 00093 this->setSampleRate (leftOperand.getSampleRate (channel)); 00094 this->setOverlap (leftOperand.getOverlap (channel)); 00095 00096 #ifdef PLONK_DEBUG 00097 const UnitType& rightOperand = this->getInputAsUnit (IOKey::RightOperand); 00098 plonk_assert (leftOperand.getBlockSize (channel) == rightOperand.getBlockSize (channel)); 00099 plonk_assert (leftOperand.getSampleRate (channel) == rightOperand.getSampleRate (channel)); 00100 plonk_assert (leftOperand.getOverlap (channel) == rightOperand.getOverlap (channel)); 00101 #endif 00102 00103 this->initValue (SampleType (0)); // todo: should really calculate this... 00104 } 00105 00106 void process (ProcessInfo& info, const int channel) throw() 00107 { 00108 UnitType& leftUnit (this->getInputAsUnit (IOKey::LeftOperand)); 00109 UnitType& rightUnit (this->getInputAsUnit (IOKey::RightOperand)); 00110 00111 const Buffer& leftBuffer (leftUnit.process (info, channel)); 00112 const Buffer& rightBuffer (rightUnit.process (info, channel)); 00113 00114 const int outputBufferLength = this->getOutputBuffer().length(); 00115 const int outputBufferLengthHalved = (unsigned int)outputBufferLength / 2; 00116 SampleType* const realOutputSamples = this->getOutputSamples(); 00117 SampleType* const imagOutputSamples = realOutputSamples + outputBufferLengthHalved; 00118 00119 plonk_assert ((outputBufferLength == leftBuffer.length()) && 00120 (outputBufferLength == rightBuffer.length())); 00121 00122 const SampleType* const leftRealSamples = leftBuffer.getArray(); 00123 const SampleType* const rightRealSamples = rightBuffer.getArray(); 00124 const SampleType* const leftImagSamples = leftRealSamples + outputBufferLengthHalved; 00125 const SampleType* const rightImagSamples = rightRealSamples + outputBufferLengthHalved; 00126 00127 const SampleType leftDC = leftRealSamples[0]; 00128 const SampleType rightDC = rightRealSamples[0]; 00129 const SampleType leftNyquist = leftImagSamples[0]; 00130 const SampleType rightNyquist = rightImagSamples[0]; 00131 00132 NumericalArrayComplex<SampleType>::zmul (realOutputSamples, imagOutputSamples, 00133 leftRealSamples, leftImagSamples, 00134 rightRealSamples, rightImagSamples, 00135 outputBufferLengthHalved); 00136 00137 realOutputSamples[0] = leftDC * rightDC; 00138 imagOutputSamples[0] = leftNyquist * rightNyquist; 00139 } 00140 }; 00141 00142 00143 00144 //------------------------------------------------------------------------------ 00145 00157 template<class SampleType> 00158 class ZMulUnit 00159 { 00160 public: 00161 typedef ZMulChannelInternal<SampleType> ZMulInternal; 00162 typedef typename ZMulInternal::Data Data; 00163 typedef UnitBase<SampleType> UnitType; 00164 typedef InputDictionary Inputs; 00165 00166 static inline UnitInfos getInfo() throw() 00167 { 00168 return UnitInfo ("ZMul", "Complex multiplication.", 00169 00170 // output 00171 ChannelCount::VariableChannelCount, 00172 IOKey::FFTPacked, Measure::FFTPacked, IOInfo::NoDefault, IOLimit::None, 00173 IOKey::End, 00174 00175 // inputs 00176 IOKey::LeftOperand, Measure::FFTPacked, IOInfo::NoDefault, IOLimit::None, 00177 IOKey::RightOperand, Measure::FFTPacked, IOInfo::NoDefault, IOLimit::None, 00178 IOKey::End); 00179 } 00180 00182 static inline UnitType ar (UnitType const& leftOperand, UnitType const& rightOperand) throw() 00183 { 00184 Inputs inputs; 00185 inputs.put (IOKey::LeftOperand, leftOperand); 00186 inputs.put (IOKey::RightOperand, rightOperand); 00187 00188 Data data = { -1.0, -1.0 }; 00189 00190 return UnitType::template createFromInputs<ZMulInternal> (inputs, 00191 data, 00192 BlockSize::noPreference(), 00193 SampleRate::noPreference()); 00194 } 00195 00196 }; 00197 00198 typedef ZMulUnit<PLONK_TYPE_DEFAULT> ZMul; 00199 00200 #endif // PLONK_ZMULCHANNEL_H 00201 00202