Dr. Stuart R. Hudson

Hi, The routine "tr00ab" computes the rotational transform on the interfaces. I suggest that this be called "iniota" or something better. Regards, Stuart Dr. Stuart R. Hudson, Acting Head, Theory...

Comment One) What you are calling "slow spec" is actually non-stellarator-symmetric spec, so comparisons between "fast spec", which exploits stellarator symmetry and slow spec are not particularly meaningful. Recall that...

Ultimately, it would be nice to determine a scaling law, e.g., T \sim L^a M^b N^c Nv^d for each subroutine, where L, M and N are the radial, poloidal and...

In http://w3.pppl.gov/~shudson/Spec/ma00aa.pdf I have color-coded the stellarator symmetric terms and the non-stellarator symmetric terms. As of last week, all the terms were being calculated regardless of Istellsym. All I did...

How are you doing the timing? I implemented a timing capability. It might be clumsy, and I don't know if it works. Set Ltiming=1 in the &diagnos namelist and a...

So, do the results from the different timings agree?

Ok, I made with CC=intel_prof, ran xspec, gprofed xspec and got this: % cumulative self self total time seconds seconds calls s/call s/call name 73.54 22.10 22.10 410 0.05 0.06...

I am still curious about the timing. Can the gnu profiling handle handle the calls to the NAG library? Also, is there a way to "profile" the memory requirements? I...

1) I think that the Ltiming option will only be reliable if only one cpu is used. Could this explain your strange result? 2) mp00ac does little more than call...

I have pushed my changes to ma00aa. I will take a break on this topic now, so this would be a good time for you to pull and merge this...