nShape = 72 DVGeo, DVCon = init_func3(nTwist) aero_groups = [] for ipt in xrange(npt): ap = init_func1('fc%02i' % ipt, A_list0[ipt], M_list0[ipt]) sys_group = Assembly('fc%02i' % ipt, input_maps={ 'twist': None, 'shape': None, }, output_maps={}, subsystems=[ IndVar('alpha', value=A_list0[ipt]), SysAeroSolver('sys_aero', ap=ap, nTwist=nTwist, nShape=nShape, DVGeo=DVGeo, init_func=init_func2), ]) aero_groups.append(sys_group) sys_aero_groups = Assembly('sys_aero_groups', subsystems=aero_groups, proc_split=[]) sys_aero_groups.nonlinear_solver = NLNsolverJC(ilimit=1) sys_aero_groups.linear_solver = LINsolverJC(ilimit=1)
init_func = pickle.load(open("../good_int_preopts/funcs_045.pkl", "rb")) DVGeo, DVCon = init_func3(nTwist) aero_groups = [] for ipt in xrange(npt): ap = init_func1('fc%02i' % ipt, A_list0[ipt], M_list0[ipt]) sys_group = Assembly('fc%02i' % ipt, input_maps={ 'twist': None, 'shape': None, }, output_maps={}, subsystems=[ IndVar('alpha', value=A_list0[ipt]), SysAeroSolver('sys_aero', ap=ap, nTwist=nTwist, nShape=nShape, DVGeo=DVGeo, init_func=init_func2), ]) aero_groups.append(sys_group) sys_aero_groups = Assembly('sys_aero_groups', subsystems=aero_groups, proc_split=[]) sys_aero_groups.nonlinear_solver = NLNsolverJC(ilimit=1) sys_aero_groups.linear_solver = LINsolverJC(ilimit=1)
nShape = 72 DVGeo, DVCon = init_func3(nTwist) aero_groups = [] for ipt in xrange(npt): ap = init_func1('fc%02i'%ipt, A_list0[ipt], M_list0[ipt]) sys_group = Assembly('fc%02i'%ipt, input_maps={ 'twist': None, 'shape': None, }, output_maps={}, subsystems=[ IndVar('alpha', value=A_list0[ipt]), SysAeroSolver('sys_aero', ap=ap, nTwist=nTwist, nShape=nShape, DVGeo=DVGeo, init_func=init_func2), ]) aero_groups.append(sys_group) sys_aero_groups = Assembly('sys_aero_groups', subsystems=aero_groups, proc_split=[]) sys_aero_groups.nonlinear_solver = NLNsolverJC(ilimit=1) sys_aero_groups.linear_solver = LINsolverJC(ilimit=1) #problem = 'problem_32rt_3ac_1new.py' #problem = 'problem_4rt_3ac_1new.py'