source: icGREP/icgrep-devel/icgrep/kernels/kernel.cpp @ 5609

Last change on this file since 5609 was 5609, checked in by cameron, 22 months ago

Ensure available/writable items for more buffer types

File size: 58.7 KB
Line 
1/*
2 *  Copyright (c) 2016 International Characters.
3 *  This software is licensed to the public under the Open Software License 3.0.
4 */
5
6#include "kernel.h"
7#include <toolchain/toolchain.h>
8#include <kernels/streamset.h>
9#include <llvm/IR/Constants.h>
10#include <llvm/IR/Function.h>
11#include <llvm/IR/Instructions.h>
12#include <llvm/IR/MDBuilder.h>
13#include <llvm/IR/Module.h>
14#include <llvm/Support/raw_ostream.h>
15#include <llvm/Bitcode/ReaderWriter.h>
16#include <llvm/Transforms/Utils/Local.h>
17#include <kernels/streamset.h>
18#include <sstream>
19#include <kernels/kernel_builder.h>
20#include <llvm/Support/Debug.h>
21
22using namespace llvm;
23using namespace parabix;
24
25namespace kernel {
26
27const std::string Kernel::DO_BLOCK_SUFFIX = "_DoBlock";
28const std::string Kernel::FINAL_BLOCK_SUFFIX = "_FinalBlock";
29const std::string Kernel::MULTI_BLOCK_SUFFIX = "_MultiBlock";
30const std::string Kernel::LOGICAL_SEGMENT_NO_SCALAR = "logicalSegNo";
31const std::string Kernel::PROCESSED_ITEM_COUNT_SUFFIX = "_processedItemCount";
32const std::string Kernel::CONSUMED_ITEM_COUNT_SUFFIX = "_consumedItemCount";
33const std::string Kernel::PRODUCED_ITEM_COUNT_SUFFIX = "_producedItemCount";
34const std::string Kernel::TERMINATION_SIGNAL = "terminationSignal";
35const std::string Kernel::BUFFER_PTR_SUFFIX = "_bufferPtr";
36const std::string Kernel::CONSUMER_SUFFIX = "_consumerLocks";
37const std::string Kernel::CYCLECOUNT_SCALAR = "CPUcycles";
38
39unsigned Kernel::addScalar(Type * const type, const std::string & name) {
40    if (LLVM_UNLIKELY(mKernelStateType != nullptr)) {
41        report_fatal_error("Cannot add field " + name + " to " + getName() + " after kernel state finalized");
42    }
43    if (LLVM_UNLIKELY(mKernelMap.count(name))) {
44        report_fatal_error(getName() + " already contains scalar field " + name);
45    }
46    const auto index = mKernelFields.size();
47    mKernelMap.emplace(name, index);
48    mKernelFields.push_back(type);
49    return index;
50}
51
52unsigned Kernel::addUnnamedScalar(Type * const type) {
53    if (LLVM_UNLIKELY(mKernelStateType != nullptr)) {
54        report_fatal_error("Cannot add unnamed field  to " + getName() + " after kernel state finalized");
55    }
56    const auto index = mKernelFields.size();
57    mKernelFields.push_back(type);
58    return index;
59}
60
61void Kernel::prepareStreamSetNameMap() {
62    for (unsigned i = 0; i < mStreamSetInputs.size(); i++) {
63        mStreamMap.emplace(mStreamSetInputs[i].name, std::make_pair(Port::Input, i));
64    }
65    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
66        mStreamMap.emplace(mStreamSetOutputs[i].name, std::make_pair(Port::Output, i));
67    }
68}
69
70void Kernel::bindPorts(const StreamSetBuffers & inputs, const StreamSetBuffers & outputs) {
71    assert (mModule == nullptr);
72    assert (mStreamSetInputBuffers.empty());
73    assert (mStreamSetOutputBuffers.empty());
74
75    if (LLVM_UNLIKELY(mStreamSetInputs.size() != inputs.size())) {
76        report_fatal_error(getName() + ": expected " + std::to_string(mStreamSetInputs.size()) +
77                           " input stream sets but was given "
78                           + std::to_string(inputs.size()));
79    }
80
81    for (unsigned i = 0; i < inputs.size(); ++i) {
82        StreamSetBuffer * const buf = inputs[i];
83        if (LLVM_UNLIKELY(buf == nullptr)) {
84            report_fatal_error(getName() + ": input stream set " + std::to_string(i)
85                               + " cannot be null");
86        }
87        buf->addConsumer(this);
88    }
89
90    if (LLVM_UNLIKELY(mStreamSetOutputs.size() != outputs.size())) {
91        report_fatal_error(getName() + ": expected " + std::to_string(mStreamSetOutputs.size())
92                           + " output stream sets but was given "
93                           + std::to_string(outputs.size()));
94    }
95
96    for (unsigned i = 0; i < outputs.size(); ++i) {
97        StreamSetBuffer * const buf = outputs[i];
98        if (LLVM_UNLIKELY(buf == nullptr)) {
99            report_fatal_error(getName() + ": output stream set " + std::to_string(i) + " cannot be null");
100        }
101        if (LLVM_LIKELY(buf->getProducer() == nullptr)) {
102            buf->setProducer(this);
103        } else {
104            report_fatal_error(getName() + ": output stream set " + std::to_string(i)
105                               + " is already produced by kernel " + buf->getProducer()->getName());
106        }
107    }
108
109    mStreamSetInputBuffers.assign(inputs.begin(), inputs.end());
110    mStreamSetOutputBuffers.assign(outputs.begin(), outputs.end());
111}
112
113Module * Kernel::makeModule(const std::unique_ptr<KernelBuilder> & idb) {
114    assert (mModule == nullptr);
115    std::stringstream cacheName;   
116    cacheName << getName() << '_' << idb->getBuilderUniqueName();
117    for (const StreamSetBuffer * b: mStreamSetInputBuffers) {
118        cacheName <<  ':' <<  b->getUniqueID();
119    }
120    for (const StreamSetBuffer * b: mStreamSetOutputBuffers) {
121        cacheName <<  ':' <<  b->getUniqueID();
122    }
123    mModule = new Module(cacheName.str(), idb->getContext());
124    prepareKernel(idb);
125    return mModule;
126}
127
128Module * Kernel::setModule(const std::unique_ptr<KernelBuilder> & idb, llvm::Module * const module) {
129    assert (mModule == nullptr);
130    mModule = module;
131    prepareKernel(idb);
132    return mModule;
133}
134
135void Kernel::prepareKernel(const std::unique_ptr<KernelBuilder> & idb) {
136    assert ("KernelBuilder does not have a valid IDISA Builder" && idb);
137    if (LLVM_UNLIKELY(mKernelStateType != nullptr)) {
138        report_fatal_error("Cannot prepare kernel after kernel state finalized");
139    }
140    const auto blockSize = idb->getBitBlockWidth();
141    if (mStride == 0) {
142        // Set the default kernel stride.
143        mStride = blockSize;
144    }
145    const auto requiredBlocks = codegen::SegmentSize + ((blockSize + mLookAheadPositions - 1) / blockSize);
146
147    for (unsigned i = 0; i < mStreamSetInputs.size(); i++) {
148        if ((mStreamSetInputBuffers[i]->getBufferBlocks() != 0) && (mStreamSetInputBuffers[i]->getBufferBlocks() < requiredBlocks)) {
149            report_fatal_error(getName() + ": " + mStreamSetInputs[i].name + " requires buffer size " + std::to_string(requiredBlocks));
150        }
151        mScalarInputs.emplace_back(mStreamSetInputBuffers[i]->getStreamSetHandle()->getType(), mStreamSetInputs[i].name + BUFFER_PTR_SUFFIX);
152        if ((i == 0) || !mStreamSetInputs[i].rate.isExact()) {
153            addScalar(idb->getSizeTy(), mStreamSetInputs[i].name + PROCESSED_ITEM_COUNT_SUFFIX);
154        }
155    }
156
157    IntegerType * const sizeTy = idb->getSizeTy();
158    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
159        mScalarInputs.emplace_back(mStreamSetOutputBuffers[i]->getStreamSetHandle()->getType(), mStreamSetOutputs[i].name + BUFFER_PTR_SUFFIX);
160        if ((mStreamSetInputs.empty() && (i == 0)) || !mStreamSetOutputs[i].rate.isExact()) {
161            addScalar(sizeTy, mStreamSetOutputs[i].name + PRODUCED_ITEM_COUNT_SUFFIX);
162        }
163    }
164    for (const auto & binding : mScalarInputs) {
165        addScalar(binding.type, binding.name);
166    }
167    for (const auto & binding : mScalarOutputs) {
168        addScalar(binding.type, binding.name);
169    }
170    if (mStreamMap.empty()) {
171        prepareStreamSetNameMap();
172    }
173    for (const auto & binding : mInternalScalars) {
174        addScalar(binding.type, binding.name);
175    }
176
177    Type * const consumerSetTy = StructType::get(sizeTy, sizeTy->getPointerTo()->getPointerTo(), nullptr)->getPointerTo();
178    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
179        addScalar(consumerSetTy, mStreamSetOutputs[i].name + CONSUMER_SUFFIX);
180    }
181
182    addScalar(sizeTy, LOGICAL_SEGMENT_NO_SCALAR);
183    addScalar(idb->getInt1Ty(), TERMINATION_SIGNAL);
184
185    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
186        addScalar(sizeTy, mStreamSetOutputs[i].name + CONSUMED_ITEM_COUNT_SUFFIX);
187    }
188
189    // We compile in a 64-bit CPU cycle counter into every kernel.   It will remain unused
190    // in normal execution, but when codegen::EnableCycleCounter is specified, pipelines
191    // will be able to add instrumentation to cached modules without recompilation.
192    addScalar(idb->getInt64Ty(), CYCLECOUNT_SCALAR);
193   
194    mKernelStateType = StructType::create(idb->getContext(), mKernelFields, getName());
195   
196    processingRateAnalysis();
197}
198   
199   
200void Kernel::processingRateAnalysis() {
201   
202    const unsigned inputSetCount = mStreamSetInputs.size();
203    const unsigned outputSetCount = mStreamSetOutputs.size();
204    const unsigned totalSetCount = inputSetCount + outputSetCount;
205   
206    mItemsPerStride.resize(totalSetCount);
207    mIsDerived.resize(totalSetCount);
208
209    mItemsPerStride[0] = mStride;
210    mIsDerived[0] = true;
211   
212    for (unsigned i = 0; i < inputSetCount; i++) {
213        // Default reference stream set is the principal input stream set.
214        auto & rate = mStreamSetInputs[i].rate;
215        if (rate.referenceStreamSet() == "") {
216            rate.setReferenceStreamSet(mStreamSetInputs[0].name);
217        }
218        Port port; unsigned ssIdx;
219        std::tie(port, ssIdx) = getStreamPort(rate.referenceStreamSet());
220        if ((port == Port::Output) || (ssIdx > i) || ((ssIdx == i) && (i > 0))) {
221            report_fatal_error(getName() + ": input set " + mStreamSetInputs[i].name + ": forward or circular rate dependency");
222        }
223        if ((rate.isExact() || rate.isMaxRatio()) && mIsDerived[ssIdx]) {
224            if ((mItemsPerStride[ssIdx] % rate.getRatioDenominator()) != 0) {
225                report_fatal_error(getName() + ": " + mStreamSetInputs[i].name + " processing rate denominator does not exactly divide items per stride.");
226            }
227            mItemsPerStride[i] = rate.calculateRatio(mItemsPerStride[ssIdx]);
228            mIsDerived[i] = rate.isExact();
229        }
230        else {
231            mIsDerived[i] = false;
232            mItemsPerStride[i] = 0;  // For unknown input rate, no items will be copied to temp buffers.
233        }
234    }
235   
236    for (unsigned i = inputSetCount; i < totalSetCount; i++) {
237        auto & rate = mStreamSetOutputs[i-inputSetCount].rate;
238        // Default reference stream set is the principal input stream set for the principal output stream set.
239        // Default reference stream set is the principal output stream set for other output stream sets.
240        if (rate.referenceStreamSet() == "") {
241            if ((mStreamSetInputs.size() > 0) && (i == inputSetCount)) {
242                rate.setReferenceStreamSet(mStreamSetInputs[0].name);
243            }
244            else {
245                rate.setReferenceStreamSet(mStreamSetOutputs[0].name);
246            }
247        }
248        Port port; unsigned ssIdx;
249        std::tie(port, ssIdx) = getStreamPort(rate.referenceStreamSet());
250        if (port == Port::Output) ssIdx += inputSetCount;
251        if ((ssIdx > i) || ((ssIdx == i) && (i > 0))) {
252            report_fatal_error(getName() + ": output set " + mStreamSetOutputs[i].name + ": forward or circular rate dependency");
253        }
254        if ((rate.isExact() || rate.isMaxRatio()) && mIsDerived[ssIdx]) {
255            if ((mItemsPerStride[ssIdx] % rate.getRatioDenominator()) != 0) {
256                report_fatal_error(getName() + ": " + mStreamSetOutputs[i-inputSetCount].name + " processing rate denominator does not exactly divide items per stride.");
257            }
258            mItemsPerStride[i] = rate.calculateRatio(mItemsPerStride[ssIdx]);
259            mIsDerived[i] = rate.isExact();
260        }
261        else {
262            mIsDerived[i] = false;
263            mItemsPerStride[i] = 0;  // For unknown output rate, no items will be copied to temp buffers.
264        }
265    }
266}
267
268
269// Default kernel signature: generate the IR and emit as byte code.
270std::string Kernel::makeSignature(const std::unique_ptr<kernel::KernelBuilder> & idb) {
271    assert ("KernelBuilder does not have a valid IDISA Builder" && idb.get());
272    if (LLVM_UNLIKELY(hasSignature())) {
273        generateKernel(idb);
274        std::string signature;
275        raw_string_ostream OS(signature);
276        WriteBitcodeToFile(getModule(), OS);
277        return signature;
278    } else {
279        return getModule()->getModuleIdentifier();
280    }
281}
282
283void Kernel::generateKernel(const std::unique_ptr<kernel::KernelBuilder> & idb) {
284    assert ("KernelBuilder does not have a valid IDISA Builder" && idb.get());
285    // If the module id cannot uniquely identify this kernel, "generateKernelSignature()" will have already
286    // generated the unoptimized IR.
287    if (!mIsGenerated) {
288        const auto m = idb->getModule();
289        const auto ip = idb->saveIP();
290        const auto saveInstance = getInstance();
291        idb->setModule(mModule);
292        addKernelDeclarations(idb);
293        callGenerateInitializeMethod(idb);
294        callGenerateDoSegmentMethod(idb);
295        callGenerateFinalizeMethod(idb);
296        setInstance(saveInstance);
297        idb->setModule(m);
298        idb->restoreIP(ip);
299        mIsGenerated = true;
300    }
301}
302
303inline void Kernel::callGenerateInitializeMethod(const std::unique_ptr<kernel::KernelBuilder> & idb) {
304    mCurrentMethod = getInitFunction(idb->getModule());
305    idb->SetInsertPoint(BasicBlock::Create(idb->getContext(), "entry", mCurrentMethod));
306    Function::arg_iterator args = mCurrentMethod->arg_begin();
307    setInstance(&*(args++));
308    idb->CreateStore(ConstantAggregateZero::get(mKernelStateType), getInstance());
309    for (const auto & binding : mScalarInputs) {
310        idb->setScalarField(binding.name, &*(args++));
311    }
312    for (const auto & binding : mStreamSetOutputs) {
313        idb->setConsumerLock(binding.name, &*(args++));
314    }
315    generateInitializeMethod(idb);
316    idb->CreateRetVoid();
317}
318
319inline void Kernel::callGenerateDoSegmentMethod(const std::unique_ptr<kernel::KernelBuilder> & idb) {
320    mCurrentMethod = getDoSegmentFunction(idb->getModule());
321    idb->SetInsertPoint(BasicBlock::Create(idb->getContext(), "entry", mCurrentMethod));
322    auto args = mCurrentMethod->arg_begin();
323    setInstance(&*(args++));
324    mIsFinal = &*(args++);
325    const auto n = mStreamSetInputs.size();
326    mAvailableItemCount.resize(n, nullptr);
327    for (unsigned i = 0; i < mStreamSetInputs.size(); i++) {
328        mAvailableItemCount[i] = &*(args++);
329    }
330    generateDoSegmentMethod(idb); // must be overridden by the KernelBuilder subtype
331    mIsFinal = nullptr;
332    mAvailableItemCount.clear();
333    idb->CreateRetVoid();
334}
335
336inline void Kernel::callGenerateFinalizeMethod(const std::unique_ptr<KernelBuilder> & idb) {
337    mCurrentMethod = getTerminateFunction(idb->getModule());
338    idb->SetInsertPoint(BasicBlock::Create(idb->getContext(), "entry", mCurrentMethod));
339    auto args = mCurrentMethod->arg_begin();
340    setInstance(&*(args++));
341    generateFinalizeMethod(idb); // may be overridden by the KernelBuilder subtype
342    const auto n = mScalarOutputs.size();
343    if (n == 0) {
344        idb->CreateRetVoid();
345    } else {
346        Value * outputs[n];
347        for (unsigned i = 0; i < n; ++i) {
348            outputs[i] = idb->getScalarField(mScalarOutputs[i].name);
349        }
350        if (n == 1) {
351            idb->CreateRet(outputs[0]);
352        } else {
353            idb->CreateAggregateRet(outputs, n);
354        }
355    }
356}
357
358unsigned Kernel::getScalarIndex(const std::string & name) const {
359    const auto f = mKernelMap.find(name);
360    if (LLVM_UNLIKELY(f == mKernelMap.end())) {
361        report_fatal_error(getName() + " does not contain scalar: " + name);
362    }
363    return f->second;
364}
365
366Value * Kernel::createInstance(const std::unique_ptr<KernelBuilder> & idb) {
367    assert ("KernelBuilder does not have a valid IDISA Builder" && idb);
368    if (LLVM_UNLIKELY(mKernelStateType == nullptr)) {
369        report_fatal_error("Cannot instantiate " + getName() + " before calling prepareKernel()");
370    }
371    setInstance(idb->CreateCacheAlignedAlloca(mKernelStateType));
372    return getInstance();
373}
374
375void Kernel::initializeInstance(const std::unique_ptr<KernelBuilder> & idb) {
376    assert ("KernelBuilder does not have a valid IDISA Builder" && idb);
377    if (LLVM_UNLIKELY(getInstance() == nullptr)) {
378        report_fatal_error("Cannot initialize " + getName() + " before calling createInstance()");
379    }
380    std::vector<Value *> args;
381    args.reserve(1 + mInitialArguments.size() + mStreamSetInputBuffers.size() + (mStreamSetOutputBuffers.size() * 2));
382    args.push_back(getInstance());
383    for (unsigned i = 0; i < mInitialArguments.size(); ++i) {
384        Value * arg = mInitialArguments[i];
385        if (LLVM_UNLIKELY(arg == nullptr)) {
386            report_fatal_error(getName() + ": initial argument " + std::to_string(i)
387                               + " cannot be null when calling createInstance()");
388        }
389        args.push_back(arg);
390    }
391    for (unsigned i = 0; i < mStreamSetInputBuffers.size(); ++i) {
392        assert (mStreamSetInputBuffers[i]);
393        Value * arg = mStreamSetInputBuffers[i]->getStreamSetHandle();
394        if (LLVM_UNLIKELY(arg == nullptr)) {
395            report_fatal_error(getName() + ": input stream set " + std::to_string(i)
396                               + " was not allocated prior to calling createInstance()");
397        }
398        args.push_back(arg);
399    }
400    assert (mStreamSetInputs.size() == mStreamSetInputBuffers.size());
401    for (unsigned i = 0; i < mStreamSetOutputBuffers.size(); ++i) {
402        assert (mStreamSetOutputBuffers[i]);
403        Value * arg = mStreamSetOutputBuffers[i]->getStreamSetHandle();
404        if (LLVM_UNLIKELY(arg == nullptr)) {
405            report_fatal_error(getName() + ": output stream set " + std::to_string(i)
406                               + " was not allocated prior to calling createInstance()");
407        }
408        args.push_back(arg);
409    }
410    assert (mStreamSetOutputs.size() == mStreamSetOutputBuffers.size());
411    IntegerType * const sizeTy = idb->getSizeTy();
412    PointerType * const sizePtrTy = sizeTy->getPointerTo();
413    PointerType * const sizePtrPtrTy = sizePtrTy->getPointerTo();
414    StructType * const consumerTy = StructType::get(sizeTy, sizePtrPtrTy, nullptr);
415    for (unsigned i = 0; i < mStreamSetOutputBuffers.size(); ++i) {
416        const auto output = mStreamSetOutputBuffers[i];
417        const auto & consumers = output->getConsumers();
418        const auto n = consumers.size();
419        AllocaInst * const outputConsumers = idb->CreateAlloca(consumerTy);
420        Value * const consumerSegNoArray = idb->CreateAlloca(ArrayType::get(sizePtrTy, n));
421        for (unsigned i = 0; i < n; ++i) {
422            Kernel * const consumer = consumers[i];
423            assert ("all instances must be created prior to initialization of any instance" && consumer->getInstance());
424            idb->setKernel(consumer);
425            Value * const segmentNoPtr = idb->getScalarFieldPtr(LOGICAL_SEGMENT_NO_SCALAR);
426            idb->CreateStore(segmentNoPtr, idb->CreateGEP(consumerSegNoArray, { idb->getInt32(0), idb->getInt32(i) }));
427        }
428        idb->setKernel(this);
429        Value * const consumerCountPtr = idb->CreateGEP(outputConsumers, {idb->getInt32(0), idb->getInt32(0)});
430        idb->CreateStore(idb->getSize(n), consumerCountPtr);
431        Value * const consumerSegNoArrayPtr = idb->CreateGEP(outputConsumers, {idb->getInt32(0), idb->getInt32(1)});
432        idb->CreateStore(idb->CreatePointerCast(consumerSegNoArray, sizePtrPtrTy), consumerSegNoArrayPtr);
433        args.push_back(outputConsumers);
434    }
435    idb->CreateCall(getInitFunction(idb->getModule()), args);
436}
437
438//  The default doSegment method dispatches to the doBlock routine for
439//  each block of the given number of blocksToDo, and then updates counts.
440
441void BlockOrientedKernel::generateDoSegmentMethod(const std::unique_ptr<KernelBuilder> & idb) {
442    BasicBlock * const entryBlock = idb->GetInsertBlock();
443    BasicBlock * const strideLoopCond = idb->CreateBasicBlock(getName() + "_strideLoopCond");
444    mStrideLoopBody = idb->CreateBasicBlock(getName() + "_strideLoopBody");
445    BasicBlock * const stridesDone = idb->CreateBasicBlock(getName() + "_stridesDone");
446    BasicBlock * const doFinalBlock = idb->CreateBasicBlock(getName() + "_doFinalBlock");
447    BasicBlock * const segmentDone = idb->CreateBasicBlock(getName() + "_segmentDone");
448
449    Value * baseTarget = nullptr;
450    if (idb->supportsIndirectBr()) {
451        baseTarget = idb->CreateSelect(mIsFinal, BlockAddress::get(doFinalBlock), BlockAddress::get(segmentDone));
452    }
453
454    ConstantInt * stride = idb->getSize(idb->getStride());
455    Value * availablePos = mAvailableItemCount[0];
456    Value * processed = idb->getProcessedItemCount(mStreamSetInputs[0].name);
457    Value * itemsAvail = idb->CreateSub(availablePos, processed);
458    Value * stridesToDo = idb->CreateUDiv(itemsAvail, stride);
459
460    idb->CreateBr(strideLoopCond);
461
462    idb->SetInsertPoint(strideLoopCond);
463
464    PHINode * branchTarget = nullptr;
465    if (idb->supportsIndirectBr()) {
466        branchTarget = idb->CreatePHI(baseTarget->getType(), 2, "branchTarget");
467        branchTarget->addIncoming(baseTarget, entryBlock);
468    }
469
470    PHINode * const stridesRemaining = idb->CreatePHI(idb->getSizeTy(), 2, "stridesRemaining");
471    stridesRemaining->addIncoming(stridesToDo, entryBlock);
472    // NOTE: stridesRemaining may go to a negative number in the final block if the generateFinalBlockMethod(...)
473    // calls CreateDoBlockMethodCall(). Do *not* replace the comparator with an unsigned one!
474    Value * notDone = idb->CreateICmpSGT(stridesRemaining, idb->getSize(0));
475    idb->CreateLikelyCondBr(notDone, mStrideLoopBody, stridesDone);
476
477    idb->SetInsertPoint(mStrideLoopBody);
478
479    if (idb->supportsIndirectBr()) {
480        mStrideLoopTarget = idb->CreatePHI(baseTarget->getType(), 2, "strideTarget");
481        mStrideLoopTarget->addIncoming(branchTarget, strideLoopCond);
482    }
483
484    /// GENERATE DO BLOCK METHOD
485
486    writeDoBlockMethod(idb);
487
488    /// UPDATE PROCESSED COUNTS
489
490    processed = idb->getProcessedItemCount(mStreamSetInputs[0].name);
491    Value * itemsDone = idb->CreateAdd(processed, stride);
492    idb->setProcessedItemCount(mStreamSetInputs[0].name, itemsDone);
493
494    stridesRemaining->addIncoming(idb->CreateSub(stridesRemaining, idb->getSize(1)), idb->GetInsertBlock());
495
496    BasicBlock * bodyEnd = idb->GetInsertBlock();
497    if (idb->supportsIndirectBr()) {
498        branchTarget->addIncoming(mStrideLoopTarget, bodyEnd);
499    }
500    idb->CreateBr(strideLoopCond);
501
502    stridesDone->moveAfter(bodyEnd);
503
504    idb->SetInsertPoint(stridesDone);
505
506    // Now conditionally perform the final block processing depending on the doFinal parameter.
507    if (idb->supportsIndirectBr()) {
508        mStrideLoopBranch = idb->CreateIndirectBr(branchTarget, 3);
509        mStrideLoopBranch->addDestination(doFinalBlock);
510        mStrideLoopBranch->addDestination(segmentDone);
511    } else {
512        idb->CreateUnlikelyCondBr(mIsFinal, doFinalBlock, segmentDone);
513    }
514
515    doFinalBlock->moveAfter(stridesDone);
516
517    idb->SetInsertPoint(doFinalBlock);
518
519    Value * remainingItems = idb->CreateSub(mAvailableItemCount[0], idb->getProcessedItemCount(mStreamSetInputs[0].name));
520
521    writeFinalBlockMethod(idb, remainingItems);
522
523    itemsDone = mAvailableItemCount[0];
524    idb->setProcessedItemCount(mStreamSetInputs[0].name, itemsDone);
525    idb->setTerminationSignal();
526    idb->CreateBr(segmentDone);
527
528    segmentDone->moveAfter(idb->GetInsertBlock());
529
530    idb->SetInsertPoint(segmentDone);
531
532    // Update the branch prediction metadata to indicate that the likely target will be segmentDone
533    if (idb->supportsIndirectBr()) {
534        MDBuilder mdb(idb->getContext());
535        const auto destinations = mStrideLoopBranch->getNumDestinations();
536        uint32_t weights[destinations];
537        for (unsigned i = 0; i < destinations; ++i) {
538            weights[i] = (mStrideLoopBranch->getDestination(i) == segmentDone) ? 100 : 1;
539        }
540        ArrayRef<uint32_t> bw(weights, destinations);
541        mStrideLoopBranch->setMetadata(LLVMContext::MD_prof, mdb.createBranchWeights(bw));
542    }
543
544}
545
546inline void BlockOrientedKernel::writeDoBlockMethod(const std::unique_ptr<KernelBuilder> & idb) {
547
548    Value * const self = getInstance();
549    Function * const cp = mCurrentMethod;
550    auto ip = idb->saveIP();
551    std::vector<Value *> availableItemCount(0);
552
553    /// Check if the do block method is called and create the function if necessary   
554    if (!idb->supportsIndirectBr()) {
555
556        std::vector<Type *> params;
557        params.reserve(1 + mAvailableItemCount.size());
558        params.push_back(self->getType());
559        for (Value * avail : mAvailableItemCount) {
560            params.push_back(avail->getType());
561        }
562
563        FunctionType * const type = FunctionType::get(idb->getVoidTy(), params, false);
564        mCurrentMethod = Function::Create(type, GlobalValue::InternalLinkage, getName() + DO_BLOCK_SUFFIX, idb->getModule());
565        mCurrentMethod->setCallingConv(CallingConv::C);
566        mCurrentMethod->setDoesNotThrow();
567        mCurrentMethod->setDoesNotCapture(1);
568        auto args = mCurrentMethod->arg_begin();
569        args->setName("self");
570        setInstance(&*args);
571        availableItemCount.reserve(mAvailableItemCount.size());
572        while (++args != mCurrentMethod->arg_end()) {
573            availableItemCount.push_back(&*args);
574        }
575        assert (availableItemCount.size() == mAvailableItemCount.size());
576        mAvailableItemCount.swap(availableItemCount);
577        idb->SetInsertPoint(BasicBlock::Create(idb->getContext(), "entry", mCurrentMethod));
578    }
579
580    std::vector<Value *> priorProduced;
581    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
582        if (isa<CircularCopybackBuffer>(mStreamSetOutputBuffers[i]) || isa<SwizzledCopybackBuffer>(mStreamSetOutputBuffers[i]))  {
583            priorProduced.push_back(idb->getProducedItemCount(mStreamSetOutputs[i].name));
584        }
585    }
586
587    generateDoBlockMethod(idb); // must be implemented by the BlockOrientedKernelBuilder subtype
588
589    unsigned priorIdx = 0;
590    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
591        Value * log2BlockSize = idb->getSize(std::log2(idb->getBitBlockWidth()));
592        if (SwizzledCopybackBuffer * const cb = dyn_cast<SwizzledCopybackBuffer>(mStreamSetOutputBuffers[i]))  {
593            BasicBlock * copyBack = idb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBack");
594            BasicBlock * done = idb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBackDone");
595            Value * newlyProduced = idb->CreateSub(idb->getProducedItemCount(mStreamSetOutputs[i].name), priorProduced[priorIdx]);
596            Value * priorBlock = idb->CreateLShr(priorProduced[priorIdx], log2BlockSize);
597            Value * priorOffset = idb->CreateAnd(priorProduced[priorIdx], idb->getSize(idb->getBitBlockWidth() - 1));
598            Value * instance = idb->getStreamSetBufferPtr(mStreamSetOutputs[i].name);
599            Value * accessibleBlocks = idb->getLinearlyAccessibleBlocks(mStreamSetOutputs[i].name, priorBlock);
600            Value * accessible = idb->CreateSub(idb->CreateShl(accessibleBlocks, log2BlockSize), priorOffset);
601            Value * wraparound = idb->CreateICmpULT(accessible, newlyProduced);
602            idb->CreateCondBr(wraparound, copyBack, done);
603            idb->SetInsertPoint(copyBack);
604            Value * copyItems = idb->CreateSub(newlyProduced, accessible);
605            cb->createCopyBack(idb.get(), instance, copyItems);
606            idb->CreateBr(done);
607            idb->SetInsertPoint(done);
608            priorIdx++;
609        }
610        if (CircularCopybackBuffer * const cb = dyn_cast<CircularCopybackBuffer>(mStreamSetOutputBuffers[i]))  {
611            BasicBlock * copyBack = idb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBack");
612            BasicBlock * done = idb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBackDone");
613            Value * instance = idb->getStreamSetBufferPtr(mStreamSetOutputs[i].name);
614            Value * newlyProduced = idb->CreateSub(idb->getProducedItemCount(mStreamSetOutputs[i].name), priorProduced[priorIdx]);
615            Value * accessible = idb->getLinearlyAccessibleItems(mStreamSetOutputs[i].name, priorProduced[priorIdx]);
616            Value * wraparound = idb->CreateICmpULT(accessible, newlyProduced);
617            idb->CreateCondBr(wraparound, copyBack, done);
618            idb->SetInsertPoint(copyBack);
619            Value * copyItems = idb->CreateSub(newlyProduced, accessible);
620            cb->createCopyBack(idb.get(), instance, copyItems);
621            idb->CreateBr(done);
622            idb->SetInsertPoint(done);
623            priorIdx++;
624        }
625    }
626
627    if (!idb->supportsIndirectBr()) {
628        // Restore the DoSegment function state then call the DoBlock method
629        idb->CreateRetVoid();
630        mDoBlockMethod = mCurrentMethod;
631        idb->restoreIP(ip);
632        setInstance(self);
633        mCurrentMethod = cp;
634        mAvailableItemCount.swap(availableItemCount);
635        CreateDoBlockMethodCall(idb);
636    }
637
638}
639
640inline void BlockOrientedKernel::writeFinalBlockMethod(const std::unique_ptr<KernelBuilder> & idb, Value * remainingItems) {
641
642    Value * const self = getInstance();
643    Function * const cp = mCurrentMethod;
644    Value * const remainingItemCount = remainingItems;
645    auto ip = idb->saveIP();
646    std::vector<Value *> availableItemCount(0);
647
648    if (!idb->supportsIndirectBr()) {
649        std::vector<Type *> params;
650        params.reserve(2 + mAvailableItemCount.size());
651        params.push_back(self->getType());
652        params.push_back(idb->getSizeTy());
653        for (Value * avail : mAvailableItemCount) {
654            params.push_back(avail->getType());
655        }
656        FunctionType * const type = FunctionType::get(idb->getVoidTy(), params, false);
657        mCurrentMethod = Function::Create(type, GlobalValue::InternalLinkage, getName() + FINAL_BLOCK_SUFFIX, idb->getModule());
658        mCurrentMethod->setCallingConv(CallingConv::C);
659        mCurrentMethod->setDoesNotThrow();
660        mCurrentMethod->setDoesNotCapture(1);
661        auto args = mCurrentMethod->arg_begin();
662        args->setName("self");
663        setInstance(&*args);
664        remainingItems = &*(++args);
665        remainingItems->setName("remainingItems");
666        availableItemCount.reserve(mAvailableItemCount.size());
667        while (++args != mCurrentMethod->arg_end()) {
668            availableItemCount.push_back(&*args);
669        }
670        assert (availableItemCount.size() == mAvailableItemCount.size());
671        mAvailableItemCount.swap(availableItemCount);
672        idb->SetInsertPoint(BasicBlock::Create(idb->getContext(), "entry", mCurrentMethod));
673    }
674
675    generateFinalBlockMethod(idb, remainingItems); // may be implemented by the BlockOrientedKernel subtype
676
677    RecursivelyDeleteTriviallyDeadInstructions(remainingItems); // if remainingItems was not used, this will eliminate it.
678
679    if (!idb->supportsIndirectBr()) {
680        idb->CreateRetVoid();
681        idb->restoreIP(ip);
682        setInstance(self);
683        mAvailableItemCount.swap(availableItemCount);
684        // Restore the DoSegment function state then call the DoFinal method
685        std::vector<Value *> args;
686        args.reserve(2 + mAvailableItemCount.size());
687        args.push_back(self);
688        args.push_back(remainingItemCount);
689        for (Value * avail : mAvailableItemCount) {
690            args.push_back(avail);
691        }
692        idb->CreateCall(mCurrentMethod, args);
693        mCurrentMethod = cp;
694    }
695
696}
697
698//  The default finalBlock method simply dispatches to the doBlock routine.
699void BlockOrientedKernel::generateFinalBlockMethod(const std::unique_ptr<KernelBuilder> & idb, Value * /* remainingItems */) {
700    CreateDoBlockMethodCall(idb);
701}
702
703void BlockOrientedKernel::CreateDoBlockMethodCall(const std::unique_ptr<KernelBuilder> & idb) {
704    if (idb->supportsIndirectBr()) {
705        BasicBlock * bb = idb->CreateBasicBlock("resume");
706        mStrideLoopBranch->addDestination(bb);
707        mStrideLoopTarget->addIncoming(BlockAddress::get(bb), idb->GetInsertBlock());
708        idb->CreateBr(mStrideLoopBody);
709        bb->moveAfter(idb->GetInsertBlock());
710        idb->SetInsertPoint(bb);
711    } else {
712        std::vector<Value *> args;
713        args.reserve(1 + mAvailableItemCount.size());
714        args.push_back(getInstance());
715        for (Value * avail : mAvailableItemCount) {
716            args.push_back(avail);
717        }
718        idb->CreateCall(mDoBlockMethod, args);
719    }
720}
721
722void MultiBlockKernel::generateDoSegmentMethod(const std::unique_ptr<KernelBuilder> & kb) {
723
724    // Stream set and buffer analysis.  When near the end of buffers
725    // or for final block processing, data for each streamset may need
726    // to be copied into temporary buffers to ensure linear access.
727    // Data is always copied as a number of whole blocks, dependent
728    // on the stream set processing rate.
729   
730    const unsigned bitBlockWidth = kb->getBitBlockWidth();
731    const unsigned inputSetCount = mStreamSetInputs.size();
732    const unsigned outputSetCount = mStreamSetOutputs.size();
733    const unsigned totalSetCount = inputSetCount + outputSetCount;
734   
735    int maxBlocksToCopy[totalSetCount];
736    for (unsigned i = 0; i < totalSetCount; i++) {
737        if (mIsDerived[i]) {
738            if (mItemsPerStride[i] % bitBlockWidth == 0) {
739                maxBlocksToCopy[i] = mItemsPerStride[i] / bitBlockWidth;
740            }
741            else {
742                // May not be block aligned, can overlap partial blocks at both ends.
743                maxBlocksToCopy[i] = mItemsPerStride[i]/bitBlockWidth + 2;
744            }
745        }
746        else {
747            // For variable input stream sets, we make a single stride of items
748            // available, if possible, but this stride could be nonaligned.
749            maxBlocksToCopy[i] = mStride / bitBlockWidth + 2;
750        }
751    }
752    auto ip = kb->saveIP();
753    Function * const cp = mCurrentMethod;
754    const auto saveInstance = getInstance();
755
756    // First prepare the multi-block method that will be used.
757
758    std::vector<Type *> multiBlockParmTypes;
759    multiBlockParmTypes.push_back(mKernelStateType->getPointerTo());
760    multiBlockParmTypes.push_back(kb->getSizeTy());
761    for (unsigned i = 1; i < mStreamSetInputs.size(); i++) {
762        if (!mIsDerived[i]) multiBlockParmTypes.push_back(kb->getSizeTy());
763    }
764    for (auto buffer : mStreamSetInputBuffers) {
765        multiBlockParmTypes.push_back(buffer->getPointerType());
766    }
767    for (auto buffer : mStreamSetOutputBuffers) {
768        multiBlockParmTypes.push_back(buffer->getPointerType());
769    }
770
771    FunctionType * const type = FunctionType::get(kb->getVoidTy(), multiBlockParmTypes, false);
772    Function * multiBlockFunction = Function::Create(type, GlobalValue::InternalLinkage, getName() + MULTI_BLOCK_SUFFIX, kb->getModule());
773    multiBlockFunction->setCallingConv(CallingConv::C);
774    multiBlockFunction->setDoesNotThrow();
775    mCurrentMethod = multiBlockFunction;
776    kb->SetInsertPoint(BasicBlock::Create(kb->getContext(), "multiBlockEntry", multiBlockFunction, 0));
777
778    auto args = multiBlockFunction->arg_begin();
779    args->setName("self");
780    setInstance(&*args);
781    (++args)->setName("itemsToDo");
782    for (unsigned i = 1; i < mStreamSetInputs.size(); i++) {
783        if (!mIsDerived[i]) (++args)->setName(mStreamSetInputs[i].name + "_availItems");
784    }
785    for (auto binding : mStreamSetInputs) {
786        (++args)->setName(binding.name + "BufPtr");
787    }
788    for (auto binding : mStreamSetOutputs) {
789        (++args)->setName(binding.name + "BufPtr");
790    }
791
792    // Now use the generateMultiBlockLogic method of the MultiBlockKernelBuilder subtype to
793    // provide the required multi-block kernel logic.
794    generateMultiBlockLogic(kb);
795
796    kb->CreateRetVoid();
797
798    kb->restoreIP(ip);
799    mCurrentMethod = cp;
800    setInstance(saveInstance);
801
802    // Now proceed with creation of the doSegment method.
803
804    BasicBlock * const entry = kb->GetInsertBlock();
805    BasicBlock * const doSegmentOuterLoop = kb->CreateBasicBlock(getName() + "_doSegmentOuterLoop");
806    BasicBlock * const doMultiBlockCall = kb->CreateBasicBlock(getName() + "_doMultiBlockCall");
807    BasicBlock * const tempBlockCheck = kb->CreateBasicBlock(getName() + "_tempBlockCheck");
808    BasicBlock * const doTempBufferBlock = kb->CreateBasicBlock(getName() + "_doTempBufferBlock");
809    BasicBlock * const segmentDone = kb->CreateBasicBlock(getName() + "_segmentDone");
810
811    Value * blockBaseMask = kb->CreateNot(kb->getSize(kb->getBitBlockWidth() - 1));
812    //
813    // Define and allocate the temporary buffer area.
814    //
815    Type * tempBuffers[totalSetCount];
816    for (unsigned i = 0; i < totalSetCount; i++) {
817        unsigned blocks = maxBlocksToCopy[i];
818        Type * bufType = i < inputSetCount ? mStreamSetInputBuffers[i]->getStreamSetBlockType() : mStreamSetOutputBuffers[i -inputSetCount]->getStreamSetBlockType();
819        if (blocks > 1) {
820            tempBuffers[i] = ArrayType::get(bufType, blocks);
821        }
822        else {
823            tempBuffers[i] = bufType;
824        }
825    }
826    Type * tempParameterStructType = StructType::create(kb->getContext(), ArrayRef<Type *>(tempBuffers, totalSetCount), "tempBuf");
827    Value * tempParameterArea = kb->CreateCacheAlignedAlloca(tempParameterStructType);
828    ConstantInt * blockSize = kb->getSize(kb->getBitBlockWidth());
829    ConstantInt * strideSize = kb->getSize(mStride);
830   
831    Value * availablePos = mAvailableItemCount[0];
832    Value * itemsAvail = availablePos;
833
834    //  Make sure that corresponding data is available depending on processing rate
835    //  for all derived input stream sets.
836    for (unsigned i = 1; i < mStreamSetInputs.size(); i++) {
837        Value * a = mAvailableItemCount[i];
838        auto & rate = mStreamSetInputs[i].rate;
839        if (mIsDerived[i]) {
840            Value * maxItems = rate.CreateMaxReferenceItemsCalculation(kb.get(), a);
841            itemsAvail = kb->CreateSelect(kb->CreateICmpULT(itemsAvail, maxItems), itemsAvail, maxItems);
842        }
843    }
844
845    Value * processed = kb->getProcessedItemCount(mStreamSetInputs[0].name);
846    Value * itemsToDo = kb->CreateSub(itemsAvail, processed);
847    Value * fullStridesToDo = kb->CreateUDiv(itemsToDo, strideSize);
848    Value * excessItems = kb->CreateURem(itemsToDo, strideSize);
849
850    //  Now we iteratively process these blocks using the doMultiBlock method.
851    //  In each iteration, we process the maximum number of linearly accessible
852    //  blocks on the principal input, reduced to ensure that the corresponding
853    //  data is linearly available at the specified processing rates for the other inputs,
854    //  and that each of the output buffers has sufficient linearly available space
855    //  (using overflow areas, if necessary) for the maximum output that can be
856    //  produced.
857
858    kb->CreateBr(doSegmentOuterLoop);
859    kb->SetInsertPoint(doSegmentOuterLoop);
860    PHINode * const stridesRemaining = kb->CreatePHI(kb->getSizeTy(), 2, "stridesRemaining");
861    stridesRemaining->addIncoming(fullStridesToDo, entry);
862
863    // For each input buffer, determine the processedItemCount, the block pointer for the
864    // buffer block containing the next item, and the number of linearly available items.
865
866    std::vector<Value *> processedItemCount;
867    std::vector<Value *> inputBlockPtr;
868    std::vector<Value *> producedItemCount;
869    std::vector<Value *> outputBlockPtr;
870
871    //  Now determine the linearly available blocks, based on blocks remaining reduced
872    //  by limitations of linearly available input buffer space.
873
874    Value * linearlyAvailStrides = stridesRemaining;
875    for (unsigned i = 0; i < inputSetCount; i++) {
876        Value * p = kb->getProcessedItemCount(mStreamSetInputs[i].name);
877        Value * blkNo = kb->CreateUDiv(p, blockSize);
878        Value * b = kb->getInputStreamBlockPtr(mStreamSetInputs[i].name, kb->getInt32(0));
879        processedItemCount.push_back(p);
880        inputBlockPtr.push_back(b);
881        auto & rate = mStreamSetInputs[i].rate;
882        if (rate.isUnknownRate()) continue;  // No calculation possible for unknown rates.
883        Value * maxReferenceItems = nullptr;
884        if ((rate.isFixedRatio()) && (rate.getRatioNumerator() == rate.getRatioDenominator())) {
885            maxReferenceItems = kb->CreateMul(kb->getLinearlyAccessibleBlocks(mStreamSetInputs[i].name, blkNo), blockSize);
886
887        } else {
888            Value * linearlyAvailItems = kb->getLinearlyAccessibleItems(mStreamSetInputs[i].name, p);
889            maxReferenceItems = rate.CreateMaxReferenceItemsCalculation(kb.get(), linearlyAvailItems);
890        }
891        Value * maxStrides = kb->CreateUDiv(maxReferenceItems, strideSize);
892        linearlyAvailStrides = kb->CreateSelect(kb->CreateICmpULT(maxStrides, linearlyAvailStrides), maxStrides, linearlyAvailStrides);
893    }
894    //  Now determine the linearly writeable blocks, based on available blocks reduced
895    //  by limitations of output buffer space.
896    Value * linearlyWritableStrides = linearlyAvailStrides;
897    for (unsigned i = 0; i < outputSetCount; i++) {
898        Value * p = kb->getProducedItemCount(mStreamSetOutputs[i].name);
899        Value * blkNo = kb->CreateUDiv(p, blockSize);
900        Value * b = kb->getOutputStreamBlockPtr(mStreamSetOutputs[i].name, kb->getInt32(0));
901        producedItemCount.push_back(p);
902        outputBlockPtr.push_back(b);
903        auto & rate = mStreamSetOutputs[i].rate;
904        if (rate.isUnknownRate()) continue;  // No calculation possible for unknown rates.
905        Value * maxReferenceItems = nullptr;
906        if ((rate.isFixedRatio()) && (rate.getRatioNumerator() == rate.getRatioDenominator())) {
907            maxReferenceItems = kb->CreateMul(kb->getLinearlyWritableBlocks(mStreamSetOutputs[i].name, blkNo), blockSize);
908        } else {
909            Value * writableItems = kb->getLinearlyWritableItems(mStreamSetOutputs[i].name, p);
910            maxReferenceItems = rate.CreateMaxReferenceItemsCalculation(kb.get(), writableItems);
911        }
912        Value * maxStrides = kb->CreateUDiv(maxReferenceItems, strideSize);
913        linearlyWritableStrides = kb->CreateSelect(kb->CreateICmpULT(maxStrides, linearlyWritableStrides), maxStrides, linearlyWritableStrides);
914    }
915    Value * haveStrides = kb->CreateICmpUGT(linearlyWritableStrides, kb->getSize(0));
916    kb->CreateCondBr(haveStrides, doMultiBlockCall, tempBlockCheck);
917
918    //  At this point we have verified the availability of one or more blocks of input data and output buffer space for all stream sets.
919    //  Now prepare the doMultiBlock call.
920    kb->SetInsertPoint(doMultiBlockCall);
921
922    Value * linearlyAvailItems = kb->CreateMul(linearlyWritableStrides, strideSize);
923
924    std::vector<Value *> doMultiBlockArgs;
925    doMultiBlockArgs.push_back(getInstance());
926    doMultiBlockArgs.push_back(linearlyAvailItems);
927    for (unsigned i = 1; i < mStreamSetInputs.size(); i++) {
928        if (!mIsDerived[i]) {
929            Value * avail = kb->CreateSub(mAvailableItemCount[i], processedItemCount[i]);
930            Value * linearlyAvail = kb->getLinearlyAccessibleItems(mStreamSetInputs[i].name, processedItemCount[i]);
931            doMultiBlockArgs.push_back(kb->CreateSelect(kb->CreateICmpULT(avail, linearlyAvail), avail, linearlyAvail));
932        }
933    }
934    for (unsigned i = 0; i < mStreamSetInputs.size(); i++) {
935        Value * bufPtr = kb->CreatePointerCast(inputBlockPtr[i], mStreamSetInputBuffers[i]->getPointerType());
936        doMultiBlockArgs.push_back(bufPtr);
937    }
938    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
939        Value * bufPtr = kb->CreatePointerCast(outputBlockPtr[i], mStreamSetOutputBuffers[i]->getPointerType());
940        doMultiBlockArgs.push_back(bufPtr);
941    }
942
943    kb->CreateCall(multiBlockFunction, doMultiBlockArgs);
944    // Do copybacks if necessary.
945    unsigned priorIdx = 0;
946    for (unsigned i = 0; i < mStreamSetOutputs.size(); i++) {
947        Value * log2BlockSize = kb->getSize(std::log2(kb->getBitBlockWidth()));
948        if (auto cb = dyn_cast<SwizzledCopybackBuffer>(mStreamSetOutputBuffers[i]))  {
949            BasicBlock * copyBack = kb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBack");
950            BasicBlock * done = kb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBackDone");
951            Value * newlyProduced = kb->CreateSub(kb->getProducedItemCount(mStreamSetOutputs[i].name), producedItemCount[i]);
952            Value * priorBlock = kb->CreateLShr(producedItemCount[i], log2BlockSize);
953            Value * priorOffset = kb->CreateAnd(producedItemCount[i], kb->getSize(kb->getBitBlockWidth() - 1));
954            Value * instance = kb->getStreamSetBufferPtr(mStreamSetOutputs[i].name);
955            Value * accessibleBlocks = kb->getLinearlyAccessibleBlocks(mStreamSetOutputs[i].name, priorBlock);
956            Value * accessible = kb->CreateSub(kb->CreateShl(accessibleBlocks, log2BlockSize), priorOffset);
957            Value * wraparound = kb->CreateICmpULT(accessible, newlyProduced);
958            kb->CreateCondBr(wraparound, copyBack, done);
959            kb->SetInsertPoint(copyBack);
960            Value * copyItems = kb->CreateSub(newlyProduced, accessible);
961            cb->createCopyBack(kb.get(), instance, copyItems);
962            kb->CreateBr(done);
963            kb->SetInsertPoint(done);
964            priorIdx++;
965        }
966        if (auto cb = dyn_cast<CircularCopybackBuffer>(mStreamSetOutputBuffers[i]))  {
967            BasicBlock * copyBack = kb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBack");
968            BasicBlock * done = kb->CreateBasicBlock(mStreamSetOutputs[i].name + "_copyBackDone");
969            Value * instance = kb->getStreamSetBufferPtr(mStreamSetOutputs[i].name);
970            Value * newlyProduced = kb->CreateSub(kb->getProducedItemCount(mStreamSetOutputs[i].name), producedItemCount[i]);
971            Value * accessible = kb->getLinearlyAccessibleItems(mStreamSetOutputs[i].name, producedItemCount[i]);
972            Value * wraparound = kb->CreateICmpULT(accessible, newlyProduced);
973            kb->CreateCondBr(wraparound, copyBack, done);
974            kb->SetInsertPoint(copyBack);
975            Value * copyItems = kb->CreateSub(newlyProduced, accessible);
976            cb->createCopyBack(kb.get(), instance, copyItems);
977            kb->CreateBr(done);
978            kb->SetInsertPoint(done);
979            priorIdx++;
980        }
981    }
982
983    Value * nowProcessed = kb->CreateAdd(processedItemCount[0], linearlyAvailItems);
984    kb->setProcessedItemCount(mStreamSetInputs[0].name, nowProcessed);
985    Value * reducedStridesToDo = kb->CreateSub(stridesRemaining, linearlyWritableStrides);
986    BasicBlock * multiBlockFinal = kb->GetInsertBlock();
987    stridesRemaining->addIncoming(reducedStridesToDo, multiBlockFinal);
988    kb->CreateBr(doSegmentOuterLoop);
989    //
990    // We use temporary buffers in 3 different cases that preclude full block processing.
991    // (a) One or more input buffers does not have a sufficient number of input items linearly available.
992    // (b) One or more output buffers does not have sufficient linearly available buffer space.
993    // (c) We have processed all the full blocks of input and only the excessItems remain.
994    // In each case we set up temporary buffers for input and output and then
995    // call the Multiblock routine.
996    //
997
998    kb->SetInsertPoint(tempBlockCheck);
999    haveStrides = kb->CreateICmpUGT(stridesRemaining, kb->getSize(0));
1000    kb->CreateCondBr(kb->CreateOr(mIsFinal, haveStrides), doTempBufferBlock, segmentDone);
1001
1002    kb->SetInsertPoint(doTempBufferBlock);
1003    Value * tempBlockItems = kb->CreateSelect(haveStrides, strideSize, excessItems);
1004    Value * doFinal = kb->CreateNot(haveStrides);
1005
1006    // Begin constructing the doMultiBlock args.
1007    std::vector<Value *> tempArgs;
1008    tempArgs.push_back(getInstance());
1009    tempArgs.push_back(tempBlockItems);
1010    // For non-derived inputs, add the available items.
1011    for (unsigned i = 1; i < mStreamSetInputs.size(); i++) {
1012        if (!mIsDerived[i]) {
1013            Value * avail = kb->CreateSub(mAvailableItemCount[i], processedItemCount[i]);
1014            tempArgs.push_back(kb->CreateSelect(kb->CreateICmpULT(avail, strideSize), avail, strideSize));
1015        }
1016    }
1017    // Prepare the temporary buffer area.
1018    //
1019    // First zero it out.
1020    Constant * const tempAreaSize = ConstantExpr::getIntegerCast(ConstantExpr::getSizeOf(tempParameterStructType), kb->getSizeTy(), false);
1021    kb->CreateMemZero(tempParameterArea, tempAreaSize);
1022    // For each input and output buffer, copy over necessary data starting from the last
1023    // block boundary.
1024    Value * finalItemCountNeeded[inputSetCount];
1025    finalItemCountNeeded[0] = kb->CreateAdd(processedItemCount[0], tempBlockItems);
1026
1027    for (unsigned i = 0; i < mStreamSetInputBuffers.size(); i++) {
1028        Type * bufPtrType = mStreamSetInputBuffers[i]->getPointerType();
1029        if (mItemsPerStride[i] != 0) {
1030            Value * tempBufPtr = kb->CreateGEP(tempParameterArea, {kb->getInt32(0), kb->getInt32(i)});
1031            tempBufPtr = kb->CreatePointerCast(tempBufPtr, bufPtrType);
1032            ConstantInt * strideItems = kb->getSize(mItemsPerStride[i]);
1033            Value * strideBasePos = kb->CreateSub(processedItemCount[i], kb->CreateURem(processedItemCount[i], strideItems));
1034            Value * blockBasePos = (mItemsPerStride[i] % bitBlockWidth == 0) ? strideBasePos : kb->CreateAnd(strideBasePos, blockBaseMask);
1035
1036            // The number of items to copy is determined by the processing rate requirements.
1037            if (i >= 1) {
1038                auto & rate = mStreamSetInputs[i].rate;
1039                std::string refSet = mStreamSetInputs[i].rate.referenceStreamSet();
1040                Port port; unsigned ssIdx;
1041                std::tie(port, ssIdx) = getStreamPort(refSet);
1042                finalItemCountNeeded[i] = rate.CreateRatioCalculation(kb.get(), finalItemCountNeeded[ssIdx], doFinal);
1043            }
1044           
1045            Value * inputPtr = kb->CreatePointerCast(kb->getRawInputPointer(mStreamSetInputs[i].name, kb->getInt32(0), blockBasePos), bufPtrType);
1046           
1047            if (maxBlocksToCopy[i] == 1) {
1048                // copy one block
1049                mStreamSetInputBuffers[i]->createBlockCopy(kb.get(), tempBufPtr, inputPtr, kb->getSize(1));
1050            }
1051            else {
1052                Value * neededItems = kb->CreateSub(finalItemCountNeeded[i], blockBasePos);
1053                Value * availFromBase = kb->getLinearlyAccessibleItems(mStreamSetInputs[i].name, blockBasePos);
1054                Value * allAvail = kb->CreateICmpULE(neededItems, availFromBase);
1055                Value * copyItems1 = kb->CreateSelect(allAvail, neededItems, availFromBase);
1056                //mStreamSetInputBuffers[i]->createBlockAlignedCopy(kb.get(), tempBufPtr, inputPtr, copyItems1);
1057                Value * copyBlocks1 = kb->CreateUDivCeil(copyItems1, blockSize);
1058                mStreamSetInputBuffers[i]->createBlockCopy(kb.get(), tempBufPtr, inputPtr, copyBlocks1);
1059                BasicBlock * copyRemaining = kb->CreateBasicBlock("copyRemaining");
1060                BasicBlock * copyDone = kb->CreateBasicBlock("copyDone");
1061                kb->CreateCondBr(allAvail, copyDone, copyRemaining);
1062                kb->SetInsertPoint(copyRemaining);
1063                Value * copyItems2 = kb->CreateSub(neededItems, copyItems1);
1064                Value * copyBlocks2 = kb->CreateUDivCeil(copyItems2, blockSize);
1065                //Value * nextBasePos = kb->CreateAdd(blockBasePos, copyItems1);
1066                Value * nextBasePos = kb->CreateAdd(blockBasePos, kb->CreateMul(copyBlocks2, blockSize));
1067                Value * nextInputPtr = kb->CreatePointerCast(kb->getRawInputPointer(mStreamSetInputs[i].name, kb->getInt32(0), nextBasePos), bufPtrType);
1068                Value * nextBufPtr = kb->CreateGEP(tempBufPtr, kb->CreateUDiv(copyItems1, blockSize));
1069                //mStreamSetInputBuffers[i]->createBlockAlignedCopy(kb.get(), nextBufPtr, nextInputPtr, copyItems2);
1070                mStreamSetInputBuffers[i]->createBlockCopy(kb.get(), nextBufPtr, nextInputPtr, copyBlocks2);
1071                kb->CreateBr(copyDone);
1072                kb->SetInsertPoint(copyDone);
1073            }
1074            tempArgs.push_back(tempBufPtr);
1075        }
1076        else {
1077            Value * bufPtr = kb->getInputStreamBlockPtr(mStreamSetInputs[i].name, kb->getInt32(0));
1078            bufPtr = kb->CreatePointerCast(bufPtr, mStreamSetInputBuffers[i]->getPointerType());
1079            tempArgs.push_back(bufPtr);           
1080        }
1081    }
1082    Value * outputBasePos[outputSetCount];
1083    for (unsigned i = 0; i < mStreamSetOutputBuffers.size(); i++) {
1084        Value * tempBufPtr = kb->CreateGEP(tempParameterArea,  {kb->getInt32(0), kb->getInt32(mStreamSetInputs.size() + i)});
1085        Type * bufPtrType = mStreamSetOutputBuffers[i]->getPointerType();
1086        tempBufPtr = kb->CreatePointerCast(tempBufPtr, bufPtrType);
1087        producedItemCount[i] = kb->getProducedItemCount(mStreamSetOutputs[i].name);
1088        outputBasePos[i] = kb->CreateAnd(producedItemCount[i], blockBaseMask);
1089        //mStreamSetOutputBuffers[i]->createBlockAlignedCopy(kb.get(), tempBufPtr, outputBlockPtr[i], kb->CreateSub(producedItemCount[i], outputBasePos[i]));
1090        Value * copyBlocks = kb->CreateUDivCeil(kb->CreateSub(producedItemCount[i], outputBasePos[i]), blockSize);
1091        mStreamSetOutputBuffers[i]->createBlockCopy(kb.get(), tempBufPtr, outputBlockPtr[i], copyBlocks);
1092        tempArgs.push_back(tempBufPtr);
1093    }
1094
1095    kb->CreateCall(multiBlockFunction, tempArgs);
1096   
1097    //  The items have been processed and output generated to the temporary areas.
1098    //  Update the processed item count (and hence all the counts derived automatically
1099    //  therefrom).
1100    kb->setProcessedItemCount(mStreamSetInputs[0].name, finalItemCountNeeded[0]);
1101   
1102    // Copy back data to the actual output buffers.
1103    for (unsigned i = 0; i < mStreamSetOutputBuffers.size(); i++) {
1104        Value * tempBufPtr = kb->CreateGEP(tempParameterArea,  {kb->getInt32(0), kb->getInt32(mStreamSetInputs.size() + i)});
1105        tempBufPtr = kb->CreatePointerCast(tempBufPtr, mStreamSetOutputBuffers[i]->getPointerType());
1106        Value * finalOutputItems = kb->getProducedItemCount(mStreamSetOutputs[i].name);
1107        Value * copyItems = kb->CreateSub(finalOutputItems, outputBasePos[i]);
1108        // Round up to exact multiple of block size.
1109        //copyItems = kb->CreateAnd(kb->CreateAdd(copyItems, kb->getSize(bitBlockWidth - 1)), blockBaseMask);
1110        Value * writableFromBase = kb->getLinearlyWritableItems(mStreamSetOutputs[i].name, outputBasePos[i]); // must be a whole number of blocks.
1111        Value * allWritable = kb->CreateICmpULE(copyItems, writableFromBase);
1112        Value * copyItems1 = kb->CreateSelect(allWritable, copyItems, writableFromBase);
1113        //mStreamSetOutputBuffers[i]->createBlockAlignedCopy(kb.get(), outputBlockPtr[i], tempBufPtr, copyItems1);
1114        Value * copyBlocks1 = kb->CreateUDivCeil(copyItems1, blockSize);
1115        mStreamSetOutputBuffers[i]->createBlockCopy(kb.get(), outputBlockPtr[i], tempBufPtr, copyBlocks1);
1116        BasicBlock * copyBackRemaining = kb->CreateBasicBlock("copyBackRemaining");
1117        BasicBlock * copyBackDone = kb->CreateBasicBlock("copyBackDone");
1118        kb->CreateCondBr(allWritable, copyBackDone, copyBackRemaining);
1119        kb->SetInsertPoint(copyBackRemaining);
1120        Value * copyItems2 = kb->CreateSub(copyItems, copyItems1);
1121        Value * nextBasePos = kb->CreateAdd(outputBasePos[i], copyItems1);
1122        Type * bufPtrType = mStreamSetOutputBuffers[i]->getPointerType();
1123        Value * nextOutputPtr = kb->CreatePointerCast(kb->getRawOutputPointer(mStreamSetOutputs[i].name, kb->getInt32(0), nextBasePos), bufPtrType);
1124        tempBufPtr = kb->CreateGEP(tempBufPtr, kb->CreateUDiv(copyItems1, blockSize));
1125        //mStreamSetOutputBuffers[i]->createBlockAlignedCopy(kb.get(), nextOutputPtr, tempBufPtr, copyItems2);
1126        Value * copyBlocks2 = kb->CreateUDivCeil(copyItems2, blockSize);
1127        mStreamSetOutputBuffers[i]->createBlockCopy(kb.get(), nextOutputPtr, tempBufPtr, copyBlocks2);
1128        kb->CreateBr(copyBackDone);
1129        kb->SetInsertPoint(copyBackDone);
1130    }
1131
1132
1133    //  We've dealt with the partial block processing and copied information back into the
1134    //  actual buffers.  If this isn't the final block, loop back for more multiblock processing.
1135    //
1136    stridesRemaining->addIncoming(kb->CreateSub(stridesRemaining, kb->CreateZExt(haveStrides, kb->getSizeTy())), kb->GetInsertBlock());
1137    BasicBlock * setTermination = kb->CreateBasicBlock("mBsetTermination");
1138    kb->CreateCondBr(haveStrides, doSegmentOuterLoop, setTermination);
1139    kb->SetInsertPoint(setTermination);
1140    kb->setTerminationSignal();
1141    kb->CreateBr(segmentDone);
1142    kb->SetInsertPoint(segmentDone);
1143}
1144
1145void Kernel::finalizeInstance(const std::unique_ptr<KernelBuilder> & idb) {
1146    assert ("KernelBuilder does not have a valid IDISA Builder" && idb);
1147    mOutputScalarResult = idb->CreateCall(getTerminateFunction(idb->getModule()), { getInstance() });
1148}
1149
1150Kernel::StreamPort Kernel::getStreamPort(const std::string & name) const {
1151    const auto f = mStreamMap.find(name);
1152    if (LLVM_UNLIKELY(f == mStreamMap.end())) {
1153        report_fatal_error(getName() + " does not contain stream set " + name);
1154    }
1155    return f->second;
1156}
1157
1158static inline std::string annotateKernelNameWithDebugFlags(std::string && name) {
1159    if (codegen::EnableAsserts) {
1160        name += "_EA";
1161    }
1162    return name;
1163}
1164
1165// CONSTRUCTOR
1166Kernel::Kernel(std::string && kernelName,
1167               std::vector<Binding> && stream_inputs,
1168               std::vector<Binding> && stream_outputs,
1169               std::vector<Binding> && scalar_parameters,
1170               std::vector<Binding> && scalar_outputs,
1171               std::vector<Binding> && internal_scalars)
1172: KernelInterface(annotateKernelNameWithDebugFlags(std::move(kernelName))
1173                  , std::move(stream_inputs), std::move(stream_outputs)
1174                  , std::move(scalar_parameters), std::move(scalar_outputs)
1175                  , std::move(internal_scalars))
1176, mCurrentMethod(nullptr)
1177, mNoTerminateAttribute(false)
1178, mIsGenerated(false)
1179, mIsFinal(nullptr)
1180, mOutputScalarResult(nullptr)
1181, mStride(0) {
1182
1183}
1184
1185Kernel::~Kernel() {
1186
1187}
1188
1189// CONSTRUCTOR
1190BlockOrientedKernel::BlockOrientedKernel(std::string && kernelName,
1191                                         std::vector<Binding> && stream_inputs,
1192                                         std::vector<Binding> && stream_outputs,
1193                                         std::vector<Binding> && scalar_parameters,
1194                                         std::vector<Binding> && scalar_outputs,
1195                                         std::vector<Binding> && internal_scalars)
1196: Kernel(std::move(kernelName), std::move(stream_inputs), std::move(stream_outputs), std::move(scalar_parameters), std::move(scalar_outputs), std::move(internal_scalars))
1197, mDoBlockMethod(nullptr)
1198, mStrideLoopBody(nullptr)
1199, mStrideLoopBranch(nullptr)
1200, mStrideLoopTarget(nullptr) {
1201
1202}
1203
1204// CONSTRUCTOR
1205MultiBlockKernel::MultiBlockKernel(std::string && kernelName,
1206                                   std::vector<Binding> && stream_inputs,
1207                                   std::vector<Binding> && stream_outputs,
1208                                   std::vector<Binding> && scalar_parameters,
1209                                   std::vector<Binding> && scalar_outputs,
1210                                   std::vector<Binding> && internal_scalars)
1211: Kernel(std::move(kernelName), std::move(stream_inputs), std::move(stream_outputs), std::move(scalar_parameters), std::move(scalar_outputs), std::move(internal_scalars)) {
1212}
1213
1214// CONSTRUCTOR
1215SegmentOrientedKernel::SegmentOrientedKernel(std::string && kernelName,
1216                                             std::vector<Binding> && stream_inputs,
1217                                             std::vector<Binding> && stream_outputs,
1218                                             std::vector<Binding> && scalar_parameters,
1219                                             std::vector<Binding> && scalar_outputs,
1220                                             std::vector<Binding> && internal_scalars)
1221: Kernel(std::move(kernelName), std::move(stream_inputs), std::move(stream_outputs), std::move(scalar_parameters), std::move(scalar_outputs), std::move(internal_scalars)) {
1222   
1223}
1224   
1225}
Note: See TracBrowser for help on using the repository browser.