Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] View of /branches/ein16/synth/d2/obj_coeff.py
ViewVC logotype

View of /branches/ein16/synth/d2/obj_coeff.py

Parent Directory Parent Directory | Revision Log Revision Log

Revision 4456 - (download) (as text) (annotate)
Thu Aug 25 04:03:27 2016 UTC (2 years, 11 months ago) by cchiw
File size: 3466 byte(s)
added missing cubic
import sys
import re
import os
from obj_apply import *
from obj_ty import *
from obj_operator import *
from obj_ex import *
from set import *
import random
import sympy
from sympy import *
x,y,z =symbols('x y z')

#-------------------- Synthetic field --------------------
# coefficients for synthetic field
# a, b*x + c*y
# d*x*x, e*x*y,f*y*y
# g*y*x*x, h*x*y*y, i*x*x*y*y
zero = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
def mk_coeffs(c):

    def get_int():
        return random.randint(lcoeff, ucoeff)
    for i in range(c):
        coeffs.append(random.randint(lcoeff, ucoeff))
    #zero = [random.randint(lcoeff, ucoeff)]
    #for i in range(c-1):
    #    zero.append(0)
    #return zero
    return coeffs

def mk_CoeffSet():
    if (coeff_style == coeff_cubic):
        return mk_coeffs(16)
    elif(coeff_style == coeff_quadratic):
        return mk_coeffs(9)+[0, 0, 0, 0, 0, 0, 0]

def mk_multiple_coeffs(n):
    for i in range(n):
    return coeffs

#get coeffs depending on dimension and coeff_order
def get_coeffs(dim):
    if (dim==0): #tensor type
        return [random.randint(lcoeff, ucoeff),0,0,0,0,0,0,0,0]
    if (dim==1):

        # 1+ x+ x^2 +x^3
        if (coeff_style == coeff_cubic):
            return mk_coeffs(4)
        elif(coeff_style == coeff_quadratic):
            # forces cubic index to be 0
            return mk_coeffs(3)+[0]
            raise Exception("coeff style ", coeff_style, "unsupported")
    if (dim==2):
        return mk_CoeffSet()
    if (dim==3):
        for i in range(dim):
        return coeffs
        #return [coeffs[0],zero+[0],zero+[0]]
        raise "dimension is not supported"

# unravel pack of coeffs
def toExp_cubic(coeff):
    # digest coefficients
    #print "inside to exp coeff",coeff
    #print "len",len(coeff)
    [a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p] = coeff
    # transform to expression
    tA = a + b*y + c*x*y+ d*x;
    tB = ((e+f*x+g*(x*x))*y*y) + (h*(x*x)*y);
    tC = i*(x*x) + (j+k*x+l*(x*x))*y*y*y;
    tD = (x*x*x)*((m*y*y*y)+(n*y*y)+(o*y)+p);
    return  tA+tB+tC+tD;

# unravel pack of coeffs
def toExp_quad(coeff, zvalue):
    # digest coefficients
    #print "inside to exp coeff",coeff
    [a,b,c,d,e,f,g,h,i] = coeff
    # transform to expression
    t0 = b*x + c*y
    t1 = d*x*x + g*y*x*x
    t2 = f*y*y + h*x*y*y
    t3 = a + e*x*y + i*(y*y)*(x*x)
    exp = zvalue*(t0+t1+t2+t3)
    return exp

def coeffToExp(coeff, dim):
    #print "coeff",coeff
    #print "dim",dim
    if (dim==0):
        zvalue = 1 # no z is use
        return toExp_quad(coeff, zvalue)
    elif (dim==1):
        [a, b, c, d] = coeff
        return a+b*x+c*x*x+d*x*x*x
    elif (dim==2):
        zvalue = 1 # no z is used
        return toExp_cubic(coeff)
        #digest coefficients
        [z0,z1,z2] = coeff
        e0 = toExp_cubic(z0) # no z is used
        e1 = toExp_cubic(z1) * z # multiply with one z
        e2 = toExp_cubic(z2) * z*z #z^2
        exp = e0+e1+e2
        #print "\n ******* exp\n e0:",e0,"\n e1:", e1,"\n e2:",e2
        return exp
    else :
        raise "unsupported field dimension"

def mk_exp(dim):
    # print "\n\nmk-exp--------------"
    coeff1= get_coeffs(dim)
    exp1 = coeffToExp(coeff1, dim)
    # print "mk-exp-tree",coeff1,exp1
    # print "exp1",exp1
    return (coeff1,exp1)

ViewVC Help
Powered by ViewVC 1.0.0