Skip to main content

initialization - Why can't I import a file in the init.m file?


I would like to have some data read in every time the kernel starts, but it doesn't seem to work. Say I have a CSV file I want to have read in, for example I use:


Export["/home/jason/temp.csv", Table[{n, 2 n, 3 n}, {n, 5}]];
Quit[];

Then I open up the init.m file in a text editor and add the line


temp123=Import["/home/jason/temp.csv"];


I save the file, then go back to a notebook and reopen the kernel and ask for the definition of temp123, it comes back as Null


?temp123

(*
 Global`temp123

 temp123=Null
*)


What is the best way to do this? I've tried putting the call to Import the CSV file in another file which is called by init.m, but the result is the same.


Thanks



Answer



I think that the reason is that -- depending on the format -- Import does either use external programs, java libraries via JLink or the frontend for imports of most "nontrivial" formats. Only some -- mainly simple ascii or plain binary -- formats and of course mathematicas own syntax are directly "imported" with pure kernel functionality. As you correctly suggested init.m is read early in the setup when the kernel not yet can make full use of all these features and thus it just isn't the right place to do that. Your own findings indicate that the setup of Import itself isn't even done before init.m is loaded. You should also note that Import is by definition a relatively slow operation (it often needs to start extra processes, load libraries, interact with the OS to read files and it furthermore tries to be smart about analyzing the file contents). Thus I would not suggest to make an import part of an init.m even if it would work, except you really enjoy the meditative tranquility when waiting for the kernel to start up :-).


What would be a relative simple way to achieve something that would provide functionality similar to what you want is to just predefine how to import the data you need in init.m, but defer the actual import to when the data is first used. This would typically be done at a later time when the kernel setup has finished and the import should work. You would just have to change the code in your init.m to:


temp123:=temp123=Import["/home/jason/temp.csv"];

init.m would then not be slowing down the startup process and the data would be imported on first use (and not at all if it isn't used in a particular session). Of course that'd fail if the data is already needed in the furhter code in init.m itself, then your current solution proposed in your own answer seems to be the best you can do. Another possibility would be to not use import but read the file as ascii text and parse the content with builtin string manipulation features which might be simple enough depending on the variability and complexity of your csv files...


Comments

Post a Comment

Popular posts from this blog

mathematical optimization - Minimizing using indices, error: Part::pkspec1: The expression cannot be used as a part specification

I want to use Minimize where the variables to minimize are indices pointing into an array. Here a MWE that hopefully shows what my problem is. vars = u@# & /@ Range[3]; cons = Flatten@ { Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; Minimize[{Total@((vec1[[#]] - vec2[[u[#]]])^2 & /@ Range[1, 3]), cons}, vars, Integers] The error I get: Part::pkspec1: The expression u[1] cannot be used as a part specification. >> Answer Ok, it seems that one can get around Mathematica trying to evaluate vec2[[u[1]]] too early by using the function Indexed[vec2,u[1]] . The working MWE would then look like the following: vars = u@# & /@ Range[3]; cons = Flatten@{ Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; NMinimize[ {Total@((vec1[[#]] - Indexed[vec2, u[#]])^2 & /@ R...
Post a Comment
Read more

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more

What is and isn't a valid variable specification for Manipulate?

I have an expression whose terms have arguments (representing subscripts), like this: myExpr = A[0] + V[1,T] I would like to put it inside a Manipulate to see its value as I move around the parameters. (The goal is eventually to plot it wrt one of the variables inside.) However, Mathematica complains when I set V[1,T] as a manipulated variable: Manipulate[Evaluate[myExpr], {A[0], 0, 1}, {V[1, T], 0, 1}] (*Manipulate::vsform: Manipulate argument {V[1,T],0,1} does not have the correct form for a variable specification. >> *) As a workaround, if I get rid of the symbol T inside the argument, it works fine: Manipulate[ Evaluate[myExpr /. T -> 15], {A[0], 0, 1}, {V[1, 15], 0, 1}] Why this behavior? Can anyone point me to the documentation that says what counts as a valid variable? And is there a way to get Manpiulate to accept an expression with a symbolic argument as a variable? Investigations I've done so far: I tried using variableQ from this answer , but it says V[1...
Post a Comment
Read more