Beispiel #1
0
def getCppSAWLCradii(numPaths = 1,\
				pathLength = 1000.0,\
			  linDensity = 1.0,\
			  persisLength = 1.0,\
			  bumped = True,\
			  locPrecision = 0.0,\
			  linkDiameter = 0.5):
	radii = np.array( PolymerCpp.getSAWLCrgs(numPaths, pathLength, linDensity,\
			    persisLength, int(bumped), locPrecision, linkDiameter) )
	if not bumped:
		return radii
	else: #search for the delimiting value of -1
		i = 0
		while (radii[i]>=0.0):
			i+=1
		return radii[i+1:-1] #return