import re import math import random PREAMBLE = """ # WARNING: This file is auto-generated. Do NOT modify it manually, but rather # modify the generating script file. Otherwise changes will be lost! """[1:] class CaseGroup(object): def __init__(self, name, description, children): self.name = name self.description = description self.children = children class ShaderCase(object): def __init__(self): pass g_processedCases = {} def indentTextBlock(text, indent): indentStr = indent * "\t" lines = text.split("\n") lines = [indentStr + line for line in lines] lines = [ ["", line][line.strip() != ""] for line in lines] return "\n".join(lines) def writeCase(f, case, indent, prefix): print " %s" % (prefix + case.name) if isinstance(case, CaseGroup): f.write(indentTextBlock('group %s "%s"\n\n' % (case.name, case.description), indent)) for child in case.children: writeCase(f, child, indent + 1, prefix + case.name + ".") f.write(indentTextBlock("\nend # %s\n" % case.name, indent)) else: # \todo [petri] Fix hack. fullPath = prefix + case.name assert (fullPath not in g_processedCases) g_processedCases[fullPath] = None f.write(indentTextBlock(str(case) + "\n", indent)) def writeAllCases(fileName, caseList): # Write all cases to file. print " %s.." % fileName f = file(fileName, "wb") f.write(PREAMBLE + "\n") for case in caseList: writeCase(f, case, 0, "") f.close() print "done! (%d cases written)" % len(g_processedCases) # Template operations. def genValues(inputs, outputs): res = [] for (name, values) in inputs: res.append("input %s = [ %s ];" % (name, " | ".join([str(v) for v in values]).lower())) for (name, values) in outputs: res.append("output %s = [ %s ];" % (name, " | ".join([str(v) for v in values]).lower())) return ("\n".join(res)) def fillTemplate(template, params): s = template for (key, value) in params.items(): m = re.search(r"^(\s*)\$\{\{%s\}\}$" % key, s, re.M) if m is not None: start = m.start(0) end = m.end(0) ws = m.group(1) if value is not None: repl = "\n".join(["%s%s" % (ws, line) for line in value.split("\n")]) s = s[:start] + repl + s[end:] else: s = s[:start] + s[end+1:] # drop the whole line else: s = s.replace("${{%s}}" % key, value) return s # Return shuffled version of list def shuffled(lst): tmp = lst[:] random.shuffle(tmp) return tmp def repeatToLength(lst, toLength): return (toLength / len(lst)) * lst + lst[: toLength % len(lst)] # Helpers to convert a list of Scalar/Vec values into another type. def toFloat(lst): return [Scalar(float(v.x)) for v in lst] def toInt(lst): return [Scalar(int(v.x)) for v in lst] def toUint(lst): return [Uint(int(v.x)) for v in lst] def toBool(lst): return [Scalar(bool(v.x)) for v in lst] def toVec4(lst): return [v.toFloat().toVec4() for v in lst] def toVec3(lst): return [v.toFloat().toVec3() for v in lst] def toVec2(lst): return [v.toFloat().toVec2() for v in lst] def toIVec4(lst): return [v.toInt().toVec4() for v in lst] def toIVec3(lst): return [v.toInt().toVec3() for v in lst] def toIVec2(lst): return [v.toInt().toVec2() for v in lst] def toBVec4(lst): return [v.toBool().toVec4() for v in lst] def toBVec3(lst): return [v.toBool().toVec3() for v in lst] def toBVec2(lst): return [v.toBool().toVec2() for v in lst] def toUVec4(lst): return [v.toUint().toUVec4() for v in lst] def toUVec3(lst): return [v.toUint().toUVec3() for v in lst] def toUVec2(lst): return [v.toUint().toUVec2() for v in lst] def toMat2(lst): return [v.toMat2() for v in lst] def toMat2x3(lst): return [v.toMat2x3() for v in lst] def toMat2x4(lst): return [v.toMat2x4() for v in lst] def toMat3x2(lst): return [v.toMat3x2() for v in lst] def toMat3(lst): return [v.toMat3() for v in lst] def toMat3x4(lst): return [v.toMat3x4() for v in lst] def toMat4x2(lst): return [v.toMat4x2() for v in lst] def toMat4x3(lst): return [v.toMat4x3() for v in lst] def toMat4(lst): return [v.toMat4() for v in lst] # Random value generation. class GenRandom(object): def __init__(self): pass def uniformVec4(self, count, mn, mx): ret = [Vec4(random.uniform(mn, mx), random.uniform(mn, mx), random.uniform(mn, mx), random.uniform(mn, mx)) for x in xrange(count)] ret[0].x = mn ret[1].x = mx ret[2].x = (mn + mx) * 0.5 return ret def uniformBVec4(self, count): ret = [Vec4(random.random() >= 0.5, random.random() >= 0.5, random.random() >= 0.5, random.random() >= 0.5) for x in xrange(count)] ret[0].x = True ret[1].x = False return ret # def uniform(self, # Math operating on Scalar/Vector types. def glslSign(a): return 0.0 if (a == 0) else +1.0 if (a > 0.0) else -1.0 def glslMod(x, y): return x - y*math.floor(x/y) def glslClamp(x, mn, mx): return mn if (x < mn) else mx if (x > mx) else x class GenMath(object): @staticmethod def unary(func): return lambda val: val.applyUnary(func) @staticmethod def binary(func): return lambda a, b: (b.expandVec(a)).applyBinary(func, a.expandVec(b)) @staticmethod def frac(val): return val.applyUnary(lambda x: x - math.floor(x)) @staticmethod def exp2(val): return val.applyUnary(lambda x: math.pow(2.0, x)) @staticmethod def log2(val): return val.applyUnary(lambda x: math.log(x, 2.0)) @staticmethod def rsq(val): return val.applyUnary(lambda x: 1.0 / math.sqrt(x)) @staticmethod def sign(val): return val.applyUnary(glslSign) @staticmethod def isEqual(a, b): return Scalar(a.isEqual(b)) @staticmethod def isNotEqual(a, b): return Scalar(not a.isEqual(b)) @staticmethod def step(a, b): return (b.expandVec(a)).applyBinary(lambda edge, x: [1.0, 0.0][x < edge], a.expandVec(b)) @staticmethod def length(a): return a.length() @staticmethod def distance(a, b): return a.distance(b) @staticmethod def dot(a, b): return a.dot(b) @staticmethod def cross(a, b): return a.cross(b) @staticmethod def normalize(a): return a.normalize() @staticmethod def boolAny(a): return a.boolAny() @staticmethod def boolAll(a): return a.boolAll() @staticmethod def boolNot(a): return a.boolNot() @staticmethod def abs(a): return a.abs() # .. class Scalar(object): def __init__(self, x): self.x = x def applyUnary(self, func): return Scalar(func(self.x)) def applyBinary(self, func, other): return Scalar(func(self.x, other.x)) def isEqual(self, other): assert isinstance(other, Scalar); return (self.x == other.x) def expandVec(self, val): return val def toScalar(self): return Scalar(self.x) def toVec2(self): return Vec2(self.x, self.x) def toVec3(self): return Vec3(self.x, self.x, self.x) def toVec4(self): return Vec4(self.x, self.x, self.x, self.x) def toUVec2(self): return UVec2(self.x, self.x) def toUVec3(self): return UVec3(self.x, self.x, self.x) def toUVec4(self): return UVec4(self.x, self.x, self.x, self.x) def toMat2(self): return Mat.fromScalar(2, 2, float(self.x)) def toMat2x3(self): return Mat.fromScalar(2, 3, float(self.x)) def toMat2x4(self): return Mat.fromScalar(2, 4, float(self.x)) def toMat3x2(self): return Mat.fromScalar(3, 2, float(self.x)) def toMat3(self): return Mat.fromScalar(3, 3, float(self.x)) def toMat3x4(self): return Mat.fromScalar(3, 4, float(self.x)) def toMat4x2(self): return Mat.fromScalar(4, 2, float(self.x)) def toMat4x3(self): return Mat.fromScalar(4, 3, float(self.x)) def toMat4(self): return Mat.fromScalar(4, 4, float(self.x)) def toFloat(self): return Scalar(float(self.x)) def toInt(self): return Scalar(int(self.x)) def toUint(self): return Uint(int(self.x)) def toBool(self): return Scalar(bool(self.x)) def getNumScalars(self): return 1 def getScalars(self): return [self.x] def typeString(self): if isinstance(self.x, bool): return "bool" elif isinstance(self.x, int): return "int" elif isinstance(self.x, float): return "float" else: assert False def vec4Swizzle(self): return "" def __str__(self): return str(self.x).lower() def __float__(self): return float(self.x) def length(self): return Scalar(abs(self.x)) def distance(self, v): assert isinstance(v, Scalar) return Scalar(abs(self.x - v.x)) def dot(self, v): assert isinstance(v, Scalar) return Scalar(self.x * v.x) def normalize(self): return Scalar(glslSign(self.x)) def abs(self): if isinstance(self.x, bool): return Scalar(self.x) else: return Scalar(abs(self.x)) def __neg__(self): return Scalar(-self.x) def __add__(self, val): if not isinstance(val, Scalar): print val assert isinstance(val, Scalar) return Scalar(self.x + val.x) def __sub__(self, val): return self + (-val) def __mul__(self, val): if isinstance(val, Scalar): return Scalar(self.x * val.x) elif isinstance(val, Vec2): return Vec2(self.x * val.x, self.x * val.y) elif isinstance(val, Vec3): return Vec3(self.x * val.x, self.x * val.y, self.x * val.z) elif isinstance(val, Vec4): return Vec4(self.x * val.x, self.x * val.y, self.x * val.z, self.x * val.w) else: assert False def __div__(self, val): if isinstance(val, Scalar): return Scalar(self.x / val.x) elif isinstance(val, Vec2): return Vec2(self.x / val.x, self.x / val.y) elif isinstance(val, Vec3): return Vec3(self.x / val.x, self.x / val.y, self.x / val.z) elif isinstance(val, Vec4): return Vec4(self.x / val.x, self.x / val.y, self.x / val.z, self.x / val.w) else: assert False class Uint(Scalar): def __init__(self, x): assert x >= 0 self.x = x def typeString(self): return "uint" def abs(self): return Scalar.abs(self).toUint() def __neg__(self): return Scalar.__neg__(self).toUint() def __add__(self, val): return Scalar.__add__(self, val).toUint() def __sub__(self, val): return self + (-val) def __mul__(self, val): return Scalar.__mul__(self, val).toUint() def __div__(self, val): return Scalar.__div__(self, val).toUint() class Vec(object): @staticmethod def fromScalarList(lst): assert (len(lst) >= 1 and len(lst) <= 4) if (len(lst) == 1): return Scalar(lst[0]) elif (len(lst) == 2): return Vec2(lst[0], lst[1]) elif (len(lst) == 3): return Vec3(lst[0], lst[1], lst[2]) else: return Vec4(lst[0], lst[1], lst[2], lst[3]) def isEqual(self, other): assert isinstance(other, Vec); return (self.getScalars() == other.getScalars()) def length(self): return Scalar(math.sqrt(self.dot(self).x)) def normalize(self): return self * Scalar(1.0 / self.length().x) def swizzle(self, indexList): inScalars = self.getScalars() outScalars = map(lambda ndx: inScalars[ndx], indexList) return Vec.fromScalarList(outScalars) def __init__(self): pass def __eq__(self, other): return self.isEqual(other) def __ne__(self, other): return not self.isEqual(other) class Vec2(Vec): def __init__(self, x, y): assert(x.__class__ == y.__class__) self.x = x self.y = y def applyUnary(self, func): return Vec2(func(self.x), func(self.y)) def applyBinary(self, func, other): return Vec2(func(self.x, other.x), func(self.y, other.y)) def expandVec(self, val): return val.toVec2() def toScalar(self): return Scalar(self.x) def toVec2(self): return Vec2(self.x, self.y) def toVec3(self): return Vec3(self.x, self.y, 0.0) def toVec4(self): return Vec4(self.x, self.y, 0.0, 0.0) def toUVec2(self): return UVec2(self.x, self.y) def toUVec3(self): return UVec3(self.x, self.y, 0.0) def toUVec4(self): return UVec4(self.x, self.y, 0.0, 0.0) def toMat2(self): return Mat2(float(self.x), 0.0, 0.0, float(self.y)); def toFloat(self): return Vec2(float(self.x), float(self.y)) def toInt(self): return Vec2(int(self.x), int(self.y)) def toUint(self): return UVec2(int(self.x), int(self.y)) def toBool(self): return Vec2(bool(self.x), bool(self.y)) def getNumScalars(self): return 2 def getScalars(self): return [self.x, self.y] def typeString(self): if isinstance(self.x, bool): return "bvec2" elif isinstance(self.x, int): return "ivec2" elif isinstance(self.x, float): return "vec2" else: assert False def vec4Swizzle(self): return ".xyxy" def __str__(self): if isinstance(self.x, bool): return "bvec2(%s, %s)" % (str(self.x).lower(), str(self.y).lower()) elif isinstance(self.x, int): return "ivec2(%i, %i)" % (self.x, self.y) elif isinstance(self.x, float): return "vec2(%s, %s)" % (self.x, self.y) else: assert False def distance(self, v): assert isinstance(v, Vec2) return (self - v).length() def dot(self, v): assert isinstance(v, Vec2) return Scalar(self.x*v.x + self.y*v.y) def abs(self): if isinstance(self.x, bool): return Vec2(self.x, self.y) else: return Vec2(abs(self.x), abs(self.y)) def __neg__(self): return Vec2(-self.x, -self.y) def __add__(self, val): if isinstance(val, Scalar): return Vec2(self.x + val, self.y + val) elif isinstance(val, Vec2): return Vec2(self.x + val.x, self.y + val.y) else: assert False def __sub__(self, val): return self + (-val) def __mul__(self, val): if isinstance(val, Scalar): val = val.toVec2() assert isinstance(val, Vec2) return Vec2(self.x * val.x, self.y * val.y) def __div__(self, val): if isinstance(val, Scalar): return Vec2(self.x / val.x, self.y / val.x) else: assert isinstance(val, Vec2) return Vec2(self.x / val.x, self.y / val.y) def boolAny(self): return Scalar(self.x or self.y) def boolAll(self): return Scalar(self.x and self.y) def boolNot(self): return Vec2(not self.x, not self.y) class UVec2(Vec2): def __init__(self, x, y): assert isinstance(x, int) and isinstance(y, int) assert x >= 0 and y >= 0 Vec2.__init__(self, x, y) def typeString(self): return "uvec2" def __str__(self): return "uvec2(%i, %i)" % (self.x, self.y) def abs(self): return Vec2.abs(self).toUint() class Vec3(Vec): def __init__(self, x, y, z): assert((x.__class__ == y.__class__) and (x.__class__ == z.__class__)) self.x = x self.y = y self.z = z def applyUnary(self, func): return Vec3(func(self.x), func(self.y), func(self.z)) def applyBinary(self, func, other): return Vec3(func(self.x, other.x), func(self.y, other.y), func(self.z, other.z)) def expandVec(self, val): return val.toVec3() def toScalar(self): return Scalar(self.x) def toVec2(self): return Vec2(self.x, self.y) def toVec3(self): return Vec3(self.x, self.y, self.z) def toVec4(self): return Vec4(self.x, self.y, self.z, 0.0) def toUVec2(self): return UVec2(self.x, self.y) def toUVec3(self): return UVec3(self.x, self.y, self.z) def toUVec4(self): return UVec4(self.x, self.y, self.z, 0.0) def toMat3(self): return Mat3(float(self.x), 0.0, 0.0, 0.0, float(self.y), 0.0, 0.0, 0.0, float(self.z)); def toFloat(self): return Vec3(float(self.x), float(self.y), float(self.z)) def toInt(self): return Vec3(int(self.x), int(self.y), int(self.z)) def toUint(self): return UVec3(int(self.x), int(self.y), int(self.z)) def toBool(self): return Vec3(bool(self.x), bool(self.y), bool(self.z)) def getNumScalars(self): return 3 def getScalars(self): return [self.x, self.y, self.z] def typeString(self): if isinstance(self.x, bool): return "bvec3" elif isinstance(self.x, int): return "ivec3" elif isinstance(self.x, float): return "vec3" else: assert False def vec4Swizzle(self): return ".xyzx" def __str__(self): if isinstance(self.x, bool): return "bvec3(%s, %s, %s)" % (str(self.x).lower(), str(self.y).lower(), str(self.z).lower()) elif isinstance(self.x, int): return "ivec3(%i, %i, %i)" % (self.x, self.y, self.z) elif isinstance(self.x, float): return "vec3(%s, %s, %s)" % (self.x, self.y, self.z) else: assert False def distance(self, v): assert isinstance(v, Vec3) return (self - v).length() def dot(self, v): assert isinstance(v, Vec3) return Scalar(self.x*v.x + self.y*v.y + self.z*v.z) def cross(self, v): assert isinstance(v, Vec3) return Vec3(self.y*v.z - v.y*self.z, self.z*v.x - v.z*self.x, self.x*v.y - v.x*self.y) def abs(self): if isinstance(self.x, bool): return Vec3(self.x, self.y, self.z) else: return Vec3(abs(self.x), abs(self.y), abs(self.z)) def __neg__(self): return Vec3(-self.x, -self.y, -self.z) def __add__(self, val): if isinstance(val, Scalar): return Vec3(self.x + val, self.y + val) elif isinstance(val, Vec3): return Vec3(self.x + val.x, self.y + val.y, self.z + val.z) else: assert False def __sub__(self, val): return self + (-val) def __mul__(self, val): if isinstance(val, Scalar): val = val.toVec3() assert isinstance(val, Vec3) return Vec3(self.x * val.x, self.y * val.y, self.z * val.z) def __div__(self, val): if isinstance(val, Scalar): return Vec3(self.x / val.x, self.y / val.x, self.z / val.x) elif isinstance(val, Vec3): return Vec3(self.x / val.x, self.y / val.y, self.z / val.z) else: assert False def boolAny(self): return Scalar(self.x or self.y or self.z) def boolAll(self): return Scalar(self.x and self.y and self.z) def boolNot(self): return Vec3(not self.x, not self.y, not self.z) class UVec3(Vec3): def __init__(self, x, y, z): assert isinstance(x, int) and isinstance(y, int) and isinstance(z, int) assert x >= 0 and y >= 0 and z >= 0 Vec3.__init__(self, x, y, z) def typeString(self): return "uvec3" def __str__(self): return "uvec3(%i, %i, %i)" % (self.x, self.y, self.z) def abs(self): return Vec3.abs(self).toUint() class Vec4(Vec): def __init__(self, x, y, z, w): assert((x.__class__ == y.__class__) and (x.__class__ == z.__class__) and (x.__class__ == w.__class__)) self.x = x self.y = y self.z = z self.w = w def applyUnary(self, func): return Vec4(func(self.x), func(self.y), func(self.z), func(self.w)) def applyBinary(self, func, other): return Vec4(func(self.x, other.x), func(self.y, other.y), func(self.z, other.z), func(self.w, other.w)) def expandVec(self, val): return val.toVec4() def toScalar(self): return Scalar(self.x) def toVec2(self): return Vec2(self.x, self.y) def toVec3(self): return Vec3(self.x, self.y, self.z) def toVec4(self): return Vec4(self.x, self.y, self.z, self.w) def toUVec2(self): return UVec2(self.x, self.y) def toUVec3(self): return UVec3(self.x, self.y, self.z) def toUVec4(self): return UVec4(self.x, self.y, self.z, self.w) def toMat2(self): return Mat2(float(self.x), float(self.y), float(self.z), float(self.w)) def toMat4(self): return Mat4(float(self.x), 0.0, 0.0, 0.0, 0.0, float(self.y), 0.0, 0.0, 0.0, 0.0, float(self.z), 0.0, 0.0, 0.0, 0.0, float(self.w)); def toFloat(self): return Vec4(float(self.x), float(self.y), float(self.z), float(self.w)) def toInt(self): return Vec4(int(self.x), int(self.y), int(self.z), int(self.w)) def toUint(self): return UVec4(int(self.x), int(self.y), int(self.z), int(self.w)) def toBool(self): return Vec4(bool(self.x), bool(self.y), bool(self.z), bool(self.w)) def getNumScalars(self): return 4 def getScalars(self): return [self.x, self.y, self.z, self.w] def typeString(self): if isinstance(self.x, bool): return "bvec4" elif isinstance(self.x, int): return "ivec4" elif isinstance(self.x, float): return "vec4" else: assert False def vec4Swizzle(self): return "" def __str__(self): if isinstance(self.x, bool): return "bvec4(%s, %s, %s, %s)" % (str(self.x).lower(), str(self.y).lower(), str(self.z).lower(), str(self.w).lower()) elif isinstance(self.x, int): return "ivec4(%i, %i, %i, %i)" % (self.x, self.y, self.z, self.w) elif isinstance(self.x, float): return "vec4(%s, %s, %s, %s)" % (self.x, self.y, self.z, self.w) else: assert False def distance(self, v): assert isinstance(v, Vec4) return (self - v).length() def dot(self, v): assert isinstance(v, Vec4) return Scalar(self.x*v.x + self.y*v.y + self.z*v.z + self.w*v.w) def abs(self): if isinstance(self.x, bool): return Vec4(self.x, self.y, self.z, self.w) else: return Vec4(abs(self.x), abs(self.y), abs(self.z), abs(self.w)) def __neg__(self): return Vec4(-self.x, -self.y, -self.z, -self.w) def __add__(self, val): if isinstance(val, Scalar): return Vec3(self.x + val, self.y + val) elif isinstance(val, Vec4): return Vec4(self.x + val.x, self.y + val.y, self.z + val.z, self.w + val.w) else: assert False def __sub__(self, val): return self + (-val) def __mul__(self, val): if isinstance(val, Scalar): val = val.toVec4() assert isinstance(val, Vec4) return Vec4(self.x * val.x, self.y * val.y, self.z * val.z, self.w * val.w) def __div__(self, val): if isinstance(val, Scalar): return Vec4(self.x / val.x, self.y / val.x, self.z / val.x, self.w / val.x) elif isinstance(val, Vec4): return Vec4(self.x / val.x, self.y / val.y, self.z / val.z, self.w / val.w) else: assert False def boolAny(self): return Scalar(self.x or self.y or self.z or self.w) def boolAll(self): return Scalar(self.x and self.y and self.z and self.w) def boolNot(self): return Vec4(not self.x, not self.y, not self.z, not self.w) class UVec4(Vec4): def __init__(self, x, y, z, w): assert isinstance(x, int) and isinstance(y, int) and isinstance(z, int) and isinstance(w, int) assert x >= 0 and y >= 0 and z >= 0 and w >= 0 Vec4.__init__(self, x, y, z, w) def typeString(self): return "uvec4" def __str__(self): return "uvec4(%i, %i, %i, %i)" % (self.x, self.y, self.z, self.w) def abs(self): return Vec4.abs(self).toUint() # \note Column-major storage. class Mat(object): def __init__ (self, numCols, numRows, scalars): assert len(scalars) == numRows*numCols self.numCols = numCols self.numRows = numRows self.scalars = scalars @staticmethod def fromScalar (numCols, numRows, scalar): scalars = [] for col in range(0, numCols): for row in range(0, numRows): scalars.append(scalar if col == row else 0.0) return Mat(numCols, numRows, scalars) @staticmethod def identity (numCols, numRows): return Mat.fromScalar(numCols, numRows, 1.0) def get (self, colNdx, rowNdx): assert 0 <= colNdx and colNdx < self.numCols assert 0 <= rowNdx and rowNdx < self.numRows return self.scalars[colNdx*self.numRows + rowNdx] def set (self, colNdx, rowNdx, scalar): assert 0 <= colNdx and colNdx < self.numCols assert 0 <= rowNdx and rowNdx < self.numRows self.scalars[colNdx*self.numRows + rowNdx] = scalar def toMatrix (self, numCols, numRows): res = Mat.identity(numCols, numRows) for col in range(0, min(self.numCols, numCols)): for row in range(0, min(self.numRows, numRows)): res.set(col, row, self.get(col, row)) return res def toMat2 (self): return self.toMatrix(2, 2) def toMat2x3 (self): return self.toMatrix(2, 3) def toMat2x4 (self): return self.toMatrix(2, 4) def toMat3x2 (self): return self.toMatrix(3, 2) def toMat3 (self): return self.toMatrix(3, 3) def toMat3x4 (self): return self.toMatrix(3, 4) def toMat4x2 (self): return self.toMatrix(4, 2) def toMat4x3 (self): return self.toMatrix(4, 3) def toMat4 (self): return self.toMatrix(4, 4) def typeString(self): if self.numRows == self.numCols: return "mat%d" % self.numRows else: return "mat%dx%d" % (self.numCols, self.numRows) def __str__(self): return "%s(%s)" % (self.typeString(), ", ".join(["%s" % s for s in self.scalars])) def isTypeEqual (self, other): return isinstance(other, Mat) and self.numRows == other.numRows and self.numCols == other.numCols def isEqual(self, other): assert self.isTypeEqual(other) return (self.scalars == other.scalars) def compMul(self, val): assert self.isTypeEqual(val) return Mat(self.numRows, self.numCols, [self.scalars(i) * val.scalars(i) for i in range(self.numRows*self.numCols)]) class Mat2(Mat): def __init__(self, m00, m01, m10, m11): Mat.__init__(self, 2, 2, [m00, m10, m01, m11]) class Mat3(Mat): def __init__(self, m00, m01, m02, m10, m11, m12, m20, m21, m22): Mat.__init__(self, 3, 3, [m00, m10, m20, m01, m11, m21, m02, m12, m22]) class Mat4(Mat): def __init__(self, m00, m01, m02, m03, m10, m11, m12, m13, m20, m21, m22, m23, m30, m31, m32, m33): Mat.__init__(self, 4, 4, [m00, m10, m20, m30, m01, m11, m21, m31, m02, m12, m22, m32, m03, m13, m23, m33])