Plotting 3D level sets in Paraview
Surfaces can be represented as certain levelsets for some 3D functions. Given a set of points with values attached to them a level set associated to a certain number will separate the points into two sets: with values higher or lower than the number considered. It is nice to have good looking plots when working with the levelset method in shape optimization. Paraview is a nice, open source framework, which has the right tools in order to produce high quality plots.
I’ll briefly describe below how to use Paraview to make some nice pictures of levelsets. First of all, you’ll need your data in some format that Paraview can understand. I use vtk file format for which there is a nice automated interface in the software I use for the optimization (FreeFem++). In the vtk file you need to have a set of points and a scalar value attached to them.
If you want to create levelsets associated to certain values, follow the steps below:

 Load your file (containing points with at least a scalar value) and click Apply.
 Next, go to Filters/Alphabetical and select “Cell Data to Point Data” (if you forget to do this, you’ll get a rough surface where you see the discretization). Click Apply.
 Then apply the Contour filter (by clicking the button or going in Filters/Alphabetical). You’ll have to select the field for which the contours will be drawn and then put in the values of the levelsets you want to see. Click Apply. An example can be seen below.
If your levelset cuts the boundary of the domain, Paraview will draw a hole there. If you want to have a closed region instead, you need to use the IsoVolume filter instead of the Contour one. The difference is that you need to specify two values and Paraview will draw the surface enclosing the points corresponding to these values. Many other features are directly available: you can color the level set following another scalar value, you can set the lighting, etc. You can also symmetrize your geometry using the Reflect filter. Below you can see a result built from my work.
You can also create animations in a pretty straightforward way. Just go to View and select the Animation box. Then you’ll see the animation options. Add a Camera object with Orbit field selected. You’ll be presented with multiple options, like the center of rotation, direction of the vertical and initial position. Once everything is set, click the play button to see the animation. Then go to File/Save Animation to save it to a file.
I heard that Paraview could to many things when dealing with visualization aspects, but I hesitated to use it until now since the interface is not straightforward. The use of Filters is not clear in the beginning, but after playing with some examples, everything becomes really easy to use. The next step is to automatize all this using scripts.
Happy New Year!
Putnam 2017 A3 – Solution
Problem A3. Denote . Then is continuous and . We can see that
Now note that on we have so . Furthermore, on we have so multiplying with we get . Therefore
To prove that goes to we can still work with . Note that the negative part cannot get too big:
As for the positive part, taking we have
Next, note that on
We would need that the last term be larger than . This is equivalent to . Since is continuous on , it is bounded above, so some delta small enough exists in order for this to work.
Grouping all of the above we get that
Since we get that goes to .
Putnam 2017 A2 – Solution
Problem A2. We have the following recurrence relation
for , given and . In order to prove that is always a polynomial with integer coefficients we should prove that divides somehow. Recurrence does not seem to work very well. Also, root based arguments might work, but you need to take good care in the computation.
A simpler idea, which is classic in this context, is to try and linearize the recurrence relation. In order to do this, let’s write two consecutive recurrence relations
We add them and we obtain the following relation
which leads straightforward to a telescoping argument. Finally, we are left with a simple linear recurrence with integer coefficient polynomials, and the result follows immediately.
Putnam 2017 – Problem A1
Problem A1. We have and . Therefore .
Next, let’s note what elements cannot be in . Note that taking square roots and squaring cannot change a nonzero remainder modulo into a zero remainder. Therefore, starting from one could never get a multiple of following the allowed operations. Thus we can safely say that multiples of are not in the minimal set .
Furthermore, could only be obtained as a square root of itself with the allowed operations. Starting from , one could never get below by performing square roots or . Therefore, the minimal set does not contain and multiples of .
Now, we show that it contains all the rest. The general idea is as follows: it is enough to find which is the smallest element in a class of remainders modulo to deduce that all larger elements are there (recall the operation ). Now in order to obtain small elements of , one would need to take successive square roots. So if we prove that for some we have for some then we get that .
Now let’s start from the beginning. We have so . Since is in for every , we get that all squares of the form greater than are in . Moreover, so all numbers of the form greater than are in . Since it follows that . Moreover, ends in and is greater than so . Next, we have which ends in and is greater than so it is also in . Therefore .
Finally, we have that , and , . This means that the minimal set is .
Putnam 2017 – Problems
Source: Art of Problem Solving forum
Problem A1. Let be the smallest set of positive integers such that
 a) is in
 b) is in whenever is in and
 c) is in whenever is in
Which positive integers are not in
(The set is “smallest” in the sense that is contained in any other such set.)
Problem A2. Let , and
for all Show that, whenever is a positive integer, is equal to a polynomial with integer coefficients.
Problem A3. Let and be real numbers with and let and be continuous functions from to such that but For every positive integer define
Show that is an increasing sequence with
Problem A4. A class with students took a quiz, on which the possible scores were Each of these scores occurred at least once, and the average score was exactly Show that the class can be divided into two groups of students in such a way that the average score for each group was exactly
Problem A5. Each of the integers from to is written on a separate card, and then the cards are combined into a deck and shuffled. Three players, and take turns in the order choosing one card at random from the deck. (Each card in the deck is equally likely to be chosen.) After a card is chosen, that card and all highernumbered cards are removed from the deck, and the remaining cards are reshuffled before the next turn. Play continues until one of the three players wins the game by drawing the card numbered
Show that for each of the three players, there are arbitrarily large values of for which that player has the highest probability among the three players of winning the game.
Problem A6. The edges of a regular icosahedron are distinguished by labeling them How many different ways are there to paint each edge red, white, or blue such that each of the 20 triangular faces of the icosahedron has two edges of the same color and a third edge of a different color?
Problem B1. Let and be distinct lines in the plane. Prove that and intersect if and only if, for every real number and every point not on or there exist points on and on such that
Problem B2. Suppose that a positive integer can be expressed as the sum of consecutive positive integers
for but for no other values of Considering all positive integers with this property, what is the smallest positive integer that occurs in any of these expressions?
Problem B3. Suppose that
is a power series for which each coefficient is or . Show that if , then must be irrational.
Problem B4. Evaluate the sum
(As usual, denotes the natural logarithm of )
Problem B5. A line in the plane of a triangle is called an equalizer if it divides into two regions having equal area and equal perimeter. Find positive integers with as small as possible, such that there exists a triangle with side lengths that has exactly two distinct equalizers.
Problem B6. Find the number of ordered tuples such that are distinct elements of and
is divisible by
A hint for Project Euler Pb 613
The text for Problem 613 can be found here. The hint is the following picture 🙂