jeffrey-hokanson / PSDR / 277

pending completion

Build # 277

Build Type

push

travis-ci-com

Committed by

jeffrey-hokanson

Commit Message

Improving keyword support

Run Details

857 of 1169 branches covered (73.31%)

Branch coverage included in aggregate %.

15 of 16 new or added lines in 1 file covered. (93.75%)

18 existing lines in 2 files now uncovered.

3593 of 4038 relevant lines covered (88.98%)

0.89 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

16.95

/psdr/demos/naca0012.py

from __future__ import print_function
import numpy as np

from psdr import Function, BoxDomain

__all__ = ['NACA0012']

class NACA0012(Function):
        r""" A test problem from OpenAeroStruct

        A test problem using OpenAeroStruct similar to that described in [JHM18]_.

        References
        ----------
        .. [JHM18] Open-source coupled aerostructural optimization using Python.
                John P. Jasa, John T. Hwang, and Joaquim R. R. A. Martins.
                Structural and Multidisciplinary Optimization (2018) 57:1815-1827
                DOI: 10.1007/s00158-018-1912-8
        
        """
        def __init__(self, n_lower = 10, n_upper = 10, fraction = 0.1):
                domain = build_hicks_henne_domain(n_lower, n_upper, fraction = fraction)
                Function.__init__(self, naca0012_func, domain, vectorized = False, kwargs = {'n_lower':n_lower, 'n_upper':n_upper})


def build_hicks_henne_domain(n_lower = 10, n_upper = 10, fraction = 0.1):
        dom = BoxDomain(-fraction*np.ones(n_lower), fraction*np.ones(n_lower), names = 'lower bump') * \
                BoxDomain(-fraction*np.ones(n_upper), fraction*np.ones(n_upper), names = 'upper bump')

        return dom


def naca0012_func(x, version = 'v1', workdir = None, verbose = False, keep_data = False, n_lower = 10, n_upper = 10):
        r"""



        """
        # If we use this inside the RedisPool, we need to load the modules
        # internal to this file
        import shutil, subprocess, os, tempfile, shlex, platform
        import numpy as np
        from subprocess import Popen, PIPE, STDOUT

        if workdir is None:
                # Docker cannot access /var by default, so we move the temporary file to
                # /tmp on MacOS
                if platform.system() == 'Darwin':
                        workdir = tempfile.mkdtemp(dir = '/tmp')
                else:
                        workdir = tempfile.mkdtemp()
        else:
                workdir = os.path.abspath(workdir)
                os.makedirs(workdir)

        # Copy the inputs to a file
        np.savetxt(workdir + '/my.input', x, fmt = '%.15e')
        
        call = "docker run -t --rm --mount  type=bind,source='%s',target='/workdir' jeffreyhokanson/naca0012:%s /workdir/my.input" % (workdir, version)
        call += " --nlower %d --nupper %d" % (n_lower, n_upper)
        args = shlex.split(call)
        with open(workdir + '/output.log', 'a') as log:
                p = Popen(args, stdout = PIPE, stderr = STDOUT)
                while True:
                        # Read output from pipe
                        # TODO: this should buffer to end of line rather than fixed size
                        output = p.stdout.readline()
                        log.write(output)

                        if verbose:
                                print(output, end ='')

                        # Check for termination
                        if p.poll() is not None:
                                break
                if p.returncode != 0:
                        print("exited with error code %d" % p.returncode)

        Y = np.loadtxt(workdir + '/my.output')

        if not keep_data:
                shutil.rmtree(workdir) 
        return Y        
        

if __name__ == '__main__':
        naca = NACA0012()
        X = naca.domain.sample(10)
        print(naca.domain_app.names)
        #Y = naca(X, verbose = True)        
        #print(Y)

1	from __future__ import print_function	4✔
2	import numpy as np	4✔
3
4	from psdr import Function, BoxDomain	4✔
5
6	__all__ = ['NACA0012']	4✔
7
8	class NACA0012(Function):	4✔
9	r""" A test problem from OpenAeroStruct
10
11	A test problem using OpenAeroStruct similar to that described in [JHM18]_.
12
13	References
14	----------
15	.. [JHM18] Open-source coupled aerostructural optimization using Python.
16	John P. Jasa, John T. Hwang, and Joaquim R. R. A. Martins.
17	Structural and Multidisciplinary Optimization (2018) 57:1815-1827
18	DOI: 10.1007/s00158-018-1912-8
19
20	"""
21	def __init__(self, n_lower = 10, n_upper = 10, fraction = 0.1):	4✔
UNCOV 22	domain = build_hicks_henne_domain(n_lower, n_upper, fraction = fraction)	×
23	Function.__init__(self, naca0012_func, domain, vectorized = False, kwargs = {'n_lower':n_lower, 'n_upper':n_upper})	×
24
25
26	def build_hicks_henne_domain(n_lower = 10, n_upper = 10, fraction = 0.1):	4✔
27	dom = BoxDomain(-fractionnp.ones(n_lower), fractionnp.ones(n_lower), names = 'lower bump') * \	×
28	BoxDomain(-fractionnp.ones(n_upper), fractionnp.ones(n_upper), names = 'upper bump')
29
UNCOV 30	return dom	×
31
32
33	def naca0012_func(x, version = 'v1', workdir = None, verbose = False, keep_data = False, n_lower = 10, n_upper = 10):	4✔
34	r"""
35
36
37
38	"""
39	# If we use this inside the RedisPool, we need to load the modules
40	# internal to this file
41	import shutil, subprocess, os, tempfile, shlex, platform	×
UNCOV 42	import numpy as np	×
UNCOV 43	from subprocess import Popen, PIPE, STDOUT	×
44
45	if workdir is None:	×
46	# Docker cannot access /var by default, so we move the temporary file to
47	# /tmp on MacOS
UNCOV 48	if platform.system() == 'Darwin':	×
49	workdir = tempfile.mkdtemp(dir = '/tmp')	×
50	else:
UNCOV 51	workdir = tempfile.mkdtemp()	×
52	else:
53	workdir = os.path.abspath(workdir)	×
UNCOV 54	os.makedirs(workdir)	×
55
56	# Copy the inputs to a file
57	np.savetxt(workdir + '/my.input', x, fmt = '%.15e')	×
58
59	call = "docker run -t --rm --mount type=bind,source='%s',target='/workdir' jeffreyhokanson/naca0012:%s /workdir/my.input" % (workdir, version)	×
UNCOV 60	call += " --nlower %d --nupper %d" % (n_lower, n_upper)	×
UNCOV 61	args = shlex.split(call)	×
62	with open(workdir + '/output.log', 'a') as log:	×
63	p = Popen(args, stdout = PIPE, stderr = STDOUT)	×
UNCOV 64	while True:	×
65	# Read output from pipe
66	# TODO: this should buffer to end of line rather than fixed size
UNCOV 67	output = p.stdout.readline()	×
UNCOV 68	log.write(output)	×
69
70	if verbose:	×
71	print(output, end ='')	×
72
73	# Check for termination
74	if p.poll() is not None:	×
UNCOV 75	break	×
76	if p.returncode != 0:	×
77	print("exited with error code %d" % p.returncode)	×
78
UNCOV 79	Y = np.loadtxt(workdir + '/my.output')	×
80
81	if not keep_data:	×
82	shutil.rmtree(workdir)	×
83	return Y	×
84
85
86	if __name__ == '__main__':	4!
UNCOV 87	naca = NACA0012()	×
UNCOV 88	X = naca.domain.sample(10)	×
UNCOV 89	print(naca.domain_app.names)	×
90	#Y = naca(X, verbose = True)
91	#print(Y)

jeffrey-hokanson / PSDR / 277

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous