-# IEEE Floating Point Divider Pipeline
+"""IEEE754 Floating Point Divider Pipeline
-from nmigen import Module
-from nmigen.cli import main, verilog
+Copyright (C) 2019 Luke Kenneth Casson Leighton <lkcl@lkcl.net>
+Copyright (C) 2019 Jacob Lifshay
+
+Relevant bugreports:
+* http://bugs.libre-riscv.org/show_bug.cgi?id=99
+* http://bugs.libre-riscv.org/show_bug.cgi?id=43
+* http://bugs.libre-riscv.org/show_bug.cgi?id=44
+
+Stack looks like this:
+
+scnorm - FPDIVSpecialCasesDeNorm ispec FPBaseData
+------ ospec FPSCData
+
+ StageChain: FPDIVSpecialCasesMod,
+ FPAddDeNormMod
+
+pipediv0 - FPDivStagesSetup ispec FPSCData
+-------- ospec DivPipeInterstageData
+
+ StageChain: FPDivPreFPAdjust,
+ DivPipeSetupStage,
+ DivPipeCalculateStage,
+ ...
+ DivPipeCalculateStage
+
+pipediv1 - FPDivStagesIntermediate ispec DivPipeInterstageData
+-------- ospec DivPipeInterstageData
+
+ StageChain: DivPipeCalculateStage,
+ ...
+ DivPipeCalculateStage
+...
+...
+
+pipediv5 - FPDivStageFinal ispec FPDivStage0Data
+-------- ospec FPPostCalcData
+
+ StageChain: DivPipeCalculateStage,
+ ...
+ DivPipeCalculateStage,
+ DivPipeFinalStage,
+ FPDivPostToFPFormat
+
+normpack - FPNormToPack ispec FPPostCalcData
+-------- ospec FPPackData
+
+ StageChain: Norm1ModSingle,
+ RoundMod,
+ CorrectionsMod,
+ PackMod
+
+the number of combinatorial StageChains (n_comb_stages) in
+FPDivStages is an argument arranged to get the length of the whole
+pipeline down to sane numbers. it specifies the number of "blocks"
+that will be combinatorially chained together.
+
+the reason for keeping the number of stages down is that for every
+pipeline clock delay, a corresponding ReservationStation is needed.
+if there are 24 pipeline stages, we need a whopping TWENTY FOUR
+RS's. that's far too many. 6 is just about an acceptable number.
+even 8 is starting to get alarmingly high.
+"""
from nmutil.singlepipe import ControlBase
from nmutil.concurrentunit import ReservationStations, num_bits
-from ieee754.fpcommon.getop import FPADDBaseData
-from ieee754.fpcommon.denorm import FPSCData
-from ieee754.fpcommon.pack import FPPackData
+from ieee754.fpcommon.fpbase import FPFormat
from ieee754.fpcommon.normtopack import FPNormToPack
-from .specialcases import FPDivSpecialCasesDeNorm
-from .divstages import FPDivStages
-
+from ieee754.fpdiv.specialcases import FPDIVSpecialCasesDeNorm
+from ieee754.fpdiv.divstages import (FPDivStagesSetup,
+ FPDivStagesIntermediate,
+ FPDivStagesFinal)
+from ieee754.pipeline import PipelineSpec
+from ieee754.div_rem_sqrt_rsqrt.core import DivPipeCoreConfig
+from nmutil.dynamicpipe import MaskCancellableRedir
class FPDIVBasePipe(ControlBase):
- def __init__(self, width, id_wid):
- ControlBase.__init__(self)
- self.pipe1 = FPDivSpecialCasesDeNorm(width, id_wid)
- self.pipe2 = FPDivStages(width, id_wid)
- self.pipe3 = FPNormToPack(width, id_wid)
+ def __init__(self, pspec):
+ self.pspec = pspec
+ ControlBase.__init__(self, maskwid=pspec.maskwid)
+
+ pipechain = []
+ # to which the answer: "as few as possible"
+ # is required. too many ReservationStations
+ # means "big problems".
+
+ # get number of stages, set up loop.
+ n_stages = pspec.core_config.n_stages
+ max_n_comb_stages = self.pspec.n_comb_stages
+ print("n_stages", n_stages)
+ stage_idx = 0
- self._eqs = self.connect([self.pipe1, self.pipe2, self.pipe3])
+ end = False
+ while not end:
+
+ n_comb_stages = max_n_comb_stages
+ # needs to convert input from pipestart ospec
+ if stage_idx == 0:
+ n_comb_stages -= 1
+ kls = FPDivStagesSetup # does n_comb_stages-1 calcs as well
+
+ # needs to convert output to pipeend ispec
+ elif stage_idx + n_comb_stages >= n_stages:
+ kls = FPDivStagesFinal # does n_comb_stages-1 calcs as well
+ end = True
+ n_comb_stages = n_stages - stage_idx
+
+ # intermediary stage
+ else:
+ kls = FPDivStagesIntermediate # does n_comb_stages calcs
+
+ # create (in each pipe) a StageChain n_comb_stages in length
+ pipechain.append(kls(self.pspec, n_comb_stages, stage_idx))
+ stage_idx += n_comb_stages # increment so that each CalcStage
+ # gets a (correct) unique index
+
+ self.pipechain = pipechain
+
+ # start and end: unpack/specialcases then normalisation/packing
+ self.pipestart = pipestart = FPDIVSpecialCasesDeNorm(self.pspec)
+ self.pipeend = pipeend = FPNormToPack(self.pspec)
+
+ self._eqs = self.connect([pipestart] + pipechain + [pipeend])
def elaborate(self, platform):
m = ControlBase.elaborate(self, platform)
- m.submodules.scnorm = self.pipe1
- m.submodules.divstages = self.pipe2
- m.submodules.normpack = self.pipe3
+
+ # add submodules
+ m.submodules.scnorm = self.pipestart
+ for i, p in enumerate(self.pipechain):
+ setattr(m.submodules, "pipediv%d" % i, p)
+ m.submodules.normpack = self.pipeend
+
+ # ControlBase.connect creates the "eqs" needed to connect each pipe
m.d.comb += self._eqs
+
return m
+def roundup(x, mod):
+ return x if x % mod == 0 else x + mod - x % mod
+
+
class FPDIVMuxInOut(ReservationStations):
""" Reservation-Station version of FPDIV pipeline.
- * fan-in on inputs (an array of FPADDBaseData: a,b,mid)
+ * fan-in on inputs (an array of FPBaseData: a,b,mid)
* N-stage divider pipeline
* fan-out on outputs (an array of FPPackData: z,mid)
Fan-in and Fan-out are combinatorial.
+
+ :op_wid: - set this to the width of an operator which can
+ then be used to change the behaviour of the pipeline.
"""
- def __init__(self, width, num_rows):
- self.width = width
- self.id_wid = num_bits(width)
- self.alu = FPDIVBasePipe(width, self.id_wid)
- ReservationStations.__init__(self, num_rows)
- def i_specfn(self):
- return FPADDBaseData(self.width, self.id_wid)
+ def __init__(self, width, num_rows, op_wid=2):
+ self.id_wid = num_bits(num_rows)
+ self.pspec = PipelineSpec(width, self.id_wid, op_wid)
+
+ # get the standard mantissa width, store in the pspec
+ fmt = FPFormat.standard(width)
+ log2_radix = 3 # tested options so far: 1, 2 and 3.
+ n_comb_stages = 2 # 2 compute stages per pipeline stage
+ maskwid = 1 # SIMD width effectively
+
+ # extra bits needed: guard + round (sticky comes from remainer.bool())
+ fraction_width = fmt.fraction_width
+ fraction_width += 2
+
+ # rounding width to a multiple of log2_radix is not needed,
+ # DivPipeCoreCalculateStage just internally reduces log2_radix on
+ # the last stage
+ cfg = DivPipeCoreConfig(fmt.width, fraction_width, log2_radix)
+
+ self.pspec.pipekls = MaskCancellableRedir
+ self.pspec.maskwid = maskwid * num_rows # RS gets just maskwid
+ self.pspec.fpformat = fmt
+ self.pspec.n_comb_stages = n_comb_stages
+ self.pspec.core_config = cfg
- def o_specfn(self):
- return FPPackData(self.width, self.id_wid)
+ # XXX TODO - a class (or function?) that takes the pspec (right here)
+ # and creates... "something". that "something" MUST have an eq function
+ # new_pspec = deepcopy(self.pspec)
+ # new_pspec.opkls = DivPipeCoreOperation
+ # self.alu = FPDIVBasePipe(new_pspec)
+ self.alu = FPDIVBasePipe(self.pspec)
+ ReservationStations.__init__(self, num_rows, maskwid=maskwid)