# on ne remplit les 3 tableaux précédents que pour les axes qui sont découpés
listOfSubdivisionDirection['0'] = 1
listOfValuesForSubdivision['0'] = [150, 10000]
listOfNumberOfPointsForSubdivision['0'] = [
3, 9
] # 3 points entre [debut, 150], 9 points entre [150, 10000], le reste pour [10000, fin]. A priori c'est inutile et redondant
listOfSubdivisionDirection['2'] = 1
listOfValuesForSubdivision['2'] = [0.71692097187]
listOfNumberOfPointsForSubdivision['2'] = [5]
###===================================================
### Pré-définir les discrétisations pour "d" (d = dimension) grilles de Tucker
###===================================================
D = discret.discretizationChoice(listOfDomainBorders, listOfNbPointsOnEachIntervalForFinalDiscretization,\
listOfSubdivisionDirection, listOfValuesForSubdivision, listOfNumberOfPointsForSubdivision)
D.getDiscretization(
) ### Update the listOfDomainBorders with the listOfSubdivisionDirection
listOfIntegralMethods = D.listOfIntegralMethods
###===================================================
### Toutes les dicrétisations fines sur "d" axes sont calculées et mettre dans "listOfTuckerGridNodes"
###===================================================
File_infor.write("===== GENERAL INFORMATION =====\n")
File_infor.write("\n")
File_infor.write("Test case: %s\n" % (testCase))
File_infor.write("\n")
File_infor.write("listOfCrossSectionNamesWithArgs %s\n" %
(listOfCrossSectionNamesWithArgs))
File_infor.write("\n")
### CHANGE HERE FOR EACH CHOICE OF LIST OF DKLIBS
###===================================================
listOfCouplageDirection = {}
listOfValuesForCouplage = {}
listOfNumberOfPointsForCouplage = {}
###===================================================
### CHANGE HERE FOR EACH DIRECTION, E.X :
###===================================================
#listOfCouplageDirection['0'] = 1
#listOfValuesForCouplage['0'] = [0.25, 0.5]
#listOfNumberOfPointsForCouplage['0'] = [9,9]
###===================================================
###===================================================
D = discret.discretizationChoice(listOfDomainBorders, listOfNbPointsOnEachIntervalForFinalDiscretization,\
listOfCouplageDirection, listOfValuesForCouplage, listOfNumberOfPointsForCouplage)
D.getDiscretization(
) ### Update the listOfDomainBorders with the listOfCouplageDirection
listOfMethodsIntegrals = D.listOfMethodsIntegrals
#print "====================="
#print "listOfMethodsIntegrals", listOfMethodsIntegrals
#print "====================="
#raw_input()
listOfTuckerGridNodes = {}
#print listOfCrossSectionNamesWithArgs
csName = listOfCrossSectionNamesWithArgs[0]
listOfKLDecompositionToDefineListOfTuckerGridNodes = {}
for Axis_k in range(dimension):