Skip to content

Commit

Permalink
Fred's Office Computer Updates to mcqmc 2022 figures
Browse files Browse the repository at this point in the history
  • Loading branch information
@fjhickernell
fjhickernell committed Mar 17, 2023
1 parent 48ae050 commit c75e3ce
Showing 1 changed file with 18 additions and 7 deletions.
25 changes: 18 additions & 7 deletions demos/talk_paper_demos/MCQMC2022_Article_Figures.ipynb
View file
Original file line number Diff line number Diff line change
Expand Up @@ -24,10 +24,21 @@
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"QMCPy Version 1.3.2\n"
"ename": "AttributeError",
"evalue": "dlsym(0x7ff9188541d0, gen_digitalnetb2): symbol not found",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mAttributeError\u001b[0m Traceback (most recent call last)",
"Cell \u001b[0;32mIn[1], line 3\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;21;01mmatplotlib\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mpyplot\u001b[39;00m \u001b[38;5;28;01mas\u001b[39;00m \u001b[38;5;21;01mplt\u001b[39;00m\n\u001b[1;32m 2\u001b[0m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;21;01mnumpy\u001b[39;00m \u001b[38;5;28;01mas\u001b[39;00m \u001b[38;5;21;01mnp\u001b[39;00m\n\u001b[0;32m----> 3\u001b[0m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;21;01mqmcpy\u001b[39;00m \u001b[38;5;28;01mas\u001b[39;00m \u001b[38;5;21;01mqp\u001b[39;00m\n\u001b[1;32m 4\u001b[0m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;21;01mtime\u001b[39;00m \u001b[38;5;66;03m#timing routines\u001b[39;00m\n\u001b[1;32m 5\u001b[0m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;21;01mwarnings\u001b[39;00m \u001b[38;5;66;03m#to suppress warnings when needed\u001b[39;00m\n",
"File \u001b[0;32m~/SoftwareRepositories/QMCSoftware/qmcpy/__init__.py:1\u001b[0m\n\u001b[0;32m----> 1\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mdiscrete_distribution\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;241m*\u001b[39m\n\u001b[1;32m 2\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mtrue_measure\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;241m*\u001b[39m\n\u001b[1;32m 3\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mintegrand\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;241m*\u001b[39m\n",
"File \u001b[0;32m~/SoftwareRepositories/QMCSoftware/qmcpy/discrete_distribution/__init__.py:4\u001b[0m\n\u001b[1;32m 2\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01miid_std_uniform\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m IIDStdUniform\n\u001b[1;32m 3\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mlattice\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m Lattice\n\u001b[0;32m----> 4\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mdigital_net_b2\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m DigitalNetB2, Sobol, SobolSciPy\n\u001b[1;32m 5\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mhalton\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m Halton\n",
"File \u001b[0;32m~/SoftwareRepositories/QMCSoftware/qmcpy/discrete_distribution/digital_net_b2/__init__.py:1\u001b[0m\n\u001b[0;32m----> 1\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mdigital_net_b2\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m DigitalNetB2,Sobol\n\u001b[1;32m 2\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01msobol_scipy\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m SobolSciPy\n",
"File \u001b[0;32m~/SoftwareRepositories/QMCSoftware/qmcpy/discrete_distribution/digital_net_b2/digital_net_b2.py:10\u001b[0m\n\u001b[1;32m 6\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01mnumpy\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;241m*\u001b[39m\n\u001b[1;32m 7\u001b[0m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;21;01mwarnings\u001b[39;00m\n\u001b[0;32m---> 10\u001b[0m \u001b[38;5;28;01mclass\u001b[39;00m \u001b[38;5;21;01mDigitalNetB2\u001b[39;00m(LD):\n\u001b[1;32m 11\u001b[0m \u001b[38;5;250m \u001b[39m\u001b[38;5;124;03m\"\"\"\u001b[39;00m\n\u001b[1;32m 12\u001b[0m \u001b[38;5;124;03m Quasi-Random digital nets in base 2.\u001b[39;00m\n\u001b[1;32m 13\u001b[0m \u001b[38;5;124;03m \u001b[39;00m\n\u001b[0;32m (...)\u001b[0m\n\u001b[1;32m 88\u001b[0m \u001b[38;5;124;03m DOI:https://doi.org/10.1145/42288.214372\u001b[39;00m\n\u001b[1;32m 89\u001b[0m \u001b[38;5;124;03m \"\"\"\u001b[39;00m\n\u001b[1;32m 91\u001b[0m dnb2_cf \u001b[38;5;241m=\u001b[39m c_lib\u001b[38;5;241m.\u001b[39mgen_digitalnetb2\n",
"File \u001b[0;32m~/SoftwareRepositories/QMCSoftware/qmcpy/discrete_distribution/digital_net_b2/digital_net_b2.py:91\u001b[0m, in \u001b[0;36mDigitalNetB2\u001b[0;34m()\u001b[0m\n\u001b[1;32m 10\u001b[0m \u001b[38;5;28;01mclass\u001b[39;00m \u001b[38;5;21;01mDigitalNetB2\u001b[39;00m(LD):\n\u001b[1;32m 11\u001b[0m \u001b[38;5;250m \u001b[39m\u001b[38;5;124;03m\"\"\"\u001b[39;00m\n\u001b[1;32m 12\u001b[0m \u001b[38;5;124;03m Quasi-Random digital nets in base 2.\u001b[39;00m\n\u001b[1;32m 13\u001b[0m \u001b[38;5;124;03m \u001b[39;00m\n\u001b[0;32m (...)\u001b[0m\n\u001b[1;32m 88\u001b[0m \u001b[38;5;124;03m DOI:https://doi.org/10.1145/42288.214372\u001b[39;00m\n\u001b[1;32m 89\u001b[0m \u001b[38;5;124;03m \"\"\"\u001b[39;00m\n\u001b[0;32m---> 91\u001b[0m dnb2_cf \u001b[38;5;241m=\u001b[39m \u001b[43mc_lib\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mgen_digitalnetb2\u001b[49m\n\u001b[1;32m 92\u001b[0m dnb2_cf\u001b[38;5;241m.\u001b[39margtypes \u001b[38;5;241m=\u001b[39m [\n\u001b[1;32m 93\u001b[0m ctypes\u001b[38;5;241m.\u001b[39mc_ulong, \u001b[38;5;66;03m# n\u001b[39;00m\n\u001b[1;32m 94\u001b[0m ctypes\u001b[38;5;241m.\u001b[39mc_ulong, \u001b[38;5;66;03m# n0\u001b[39;00m\n\u001b[0;32m (...)\u001b[0m\n\u001b[1;32m 102\u001b[0m ctypeslib\u001b[38;5;241m.\u001b[39mndpointer(ctypes\u001b[38;5;241m.\u001b[39mc_double, flags\u001b[38;5;241m=\u001b[39m\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mC_CONTIGUOUS\u001b[39m\u001b[38;5;124m'\u001b[39m), \u001b[38;5;66;03m# x\u001b[39;00m\n\u001b[1;32m 103\u001b[0m ctypeslib\u001b[38;5;241m.\u001b[39mndpointer(ctypes\u001b[38;5;241m.\u001b[39mc_double, flags\u001b[38;5;241m=\u001b[39m\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mC_CONTIGUOUS\u001b[39m\u001b[38;5;124m'\u001b[39m)] \u001b[38;5;66;03m# xr\u001b[39;00m\n\u001b[1;32m 104\u001b[0m dnb2_cf\u001b[38;5;241m.\u001b[39mrestype \u001b[38;5;241m=\u001b[39m ctypes\u001b[38;5;241m.\u001b[39mc_uint32\n",
"File \u001b[0;32m~/miniconda3/envs/qmcpy/lib/python3.9/ctypes/__init__.py:387\u001b[0m, in \u001b[0;36mCDLL.__getattr__\u001b[0;34m(self, name)\u001b[0m\n\u001b[1;32m 385\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m name\u001b[38;5;241m.\u001b[39mstartswith(\u001b[38;5;124m'\u001b[39m\u001b[38;5;124m__\u001b[39m\u001b[38;5;124m'\u001b[39m) \u001b[38;5;129;01mand\u001b[39;00m name\u001b[38;5;241m.\u001b[39mendswith(\u001b[38;5;124m'\u001b[39m\u001b[38;5;124m__\u001b[39m\u001b[38;5;124m'\u001b[39m):\n\u001b[1;32m 386\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mAttributeError\u001b[39;00m(name)\n\u001b[0;32m--> 387\u001b[0m func \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[38;5;21;43m__getitem__\u001b[39;49m\u001b[43m(\u001b[49m\u001b[43mname\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 388\u001b[0m \u001b[38;5;28msetattr\u001b[39m(\u001b[38;5;28mself\u001b[39m, name, func)\n\u001b[1;32m 389\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m func\n",
"File \u001b[0;32m~/miniconda3/envs/qmcpy/lib/python3.9/ctypes/__init__.py:392\u001b[0m, in \u001b[0;36mCDLL.__getitem__\u001b[0;34m(self, name_or_ordinal)\u001b[0m\n\u001b[1;32m 391\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21m__getitem__\u001b[39m(\u001b[38;5;28mself\u001b[39m, name_or_ordinal):\n\u001b[0;32m--> 392\u001b[0m func \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43m_FuncPtr\u001b[49m\u001b[43m(\u001b[49m\u001b[43m(\u001b[49m\u001b[43mname_or_ordinal\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m)\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 393\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28misinstance\u001b[39m(name_or_ordinal, \u001b[38;5;28mint\u001b[39m):\n\u001b[1;32m 394\u001b[0m func\u001b[38;5;241m.\u001b[39m\u001b[38;5;18m__name__\u001b[39m \u001b[38;5;241m=\u001b[39m name_or_ordinal\n",
"\u001b[0;31mAttributeError\u001b[0m: dlsym(0x7ff9188541d0, gen_digitalnetb2): symbol not found"
]
}
],
Expand Down Expand Up @@ -122,7 +133,7 @@
"n = 16 #number of points\n",
"ld = qp.Lattice(d) #define the generator\n",
"fig = plot_successive_points(ld,'Lattice',first_n=n,n_cols=4)\n",
"fig.savefig(figpath+'latticepts.eps',format='eps')"
"fig.savefig(figpath+'latticepts.eps',format='eps',bbox_inches='tight')"
]
},
{
Expand Down Expand Up @@ -154,7 +165,7 @@
"ld = qp.Sobol(d) #define the generator\n",
"xpts_Sobol = ld.gen_samples(n) #generate points\n",
"fig = plot_successive_points(ld,'Sobol\\'',first_n=n,n_cols=4)\n",
"fig.savefig(figpath+'sobolpts.eps',format='eps')"
"fig.savefig(figpath+'sobolpts.eps',format='eps',bbox_inches='tight')"
]
},
{
Expand Down Expand Up @@ -187,7 +198,7 @@
"source": [
"iid = qp.IIDStdUniform(d) #define the generator\n",
"fig = plot_successive_points(iid,'IID',first_n=n,n_cols=4)\n",
"fig.savefig(figpath+'iidpts.eps',format='eps')"
"fig.savefig(figpath+'iidpts.eps',format='eps',bbox_inches='tight')"
]
},
{
Expand Down

0 comments on commit c75e3ce

Please sign in to comment.