loadODF_VPSC

[creuzige]

I've started to use the VPSC functions you've included in mtex. Overall they work pretty well (Thanks!), but I encountered a few issues to share.

• The loadODF_VPSC function seems to only be looking for the string "TEXTURE AT STRAIN" on the first line to determine format. Many input files (weights or .wts) to VPSC which are the same format may not have this name. This string tends to be written out only when running VPSC. It may be better to remove this format check from this program.

• The function export_VPSC does not include the line "TEXTURE AT STRAIN", so you can't read files back in. I added: fprintf(fid,'TEXTURE AT STRAIN\n\n\nB %d\n',length(ori)); to line 35 of export_VPSC (orientations) as a work around on my system.

• One of the options in VPSC is to output the textures (as weights files) at intermediate deformation points (IWRITE variable in VPSC). This results in a file that has several lines "TEXTURE AT STRAIN" distributed throughout the file. When I've tried to read these in, the loadODF_VPSC program fails. Is it possible to scan through the file, then pick which step to use? Or create ODFs for all of them?

• The 4th line indicates the angular convention. The options used in VPSC are Bunge (B), Roe (R), and Kocks (K). The number that follows are the number of grains. This could be used to truncate readin if you wanted to import more deformation steps, or to set the angular convention at the time of readin. see http://public.lanl.gov/lebenso/VPSC7c_manual.pdf, page 51-52

• For output files, other columns may be included (grain Von Mises strain, Von Mises stress, accumulated plastic work, and Taylor Factor). Not sure if that will crash things or not right now, or if they will just not be included. see http://public.lanl.gov/lebenso/VPSC7c_manual.pdf, page 58-59

• I can see that this part of the code hasn't been as extensively documented. Notes that an ODF will be returned, so you need to set symmetry would be helpful. I'm willing to help out with this, but I'm not sure how to best to do so. cs=crystalSymmetry('mmm'); ss=specimenSymmetry('mmm'); loadODF_VPSC('TEX_PH3-last.txt','cs',cs,'ss',ss)