Skip to search.

Breaking News Visit Yahoo! News for the latest.

×Close this window

simpson-simmol · The SIMPSON - SIMMOL Discussion Forum

The Yahoo! Groups Product Blog

Check it out!

Group Information

  • Members: 363
  • Category: Chemistry
  • Founded: Dec 10, 2001
  • Language: English
? Already a member? Sign in to Yahoo!

Yahoo! Groups Tips

Did you know...
Message search is now enhanced, find messages faster. Take it for a spin.

Messages

Advanced
Messages Help
Messages 919 - 948 of 1016   Oldest  |  < Older  |  Newer >  |  Newest
Messages: Show Message Summaries Sort by Date ^  
#919 From: Bingwen Hu <bw.hu@...>
Date: Sat Mar 12, 2011 2:23 am
Subject: simpson 3.0 for windows doesn't support cluster mode
bw.hu
Send Email Send Email
 
Dear developper,

     It looks that simpson 3.0 for windows doesn't support cluster mode any more. Why? SIMPSON 1.1 supports.

Best
Bingwen


#920 From: Thomas Vosegaard <tv@...>
Date: Sat Mar 12, 2011 5:15 am
Subject: Re: simpson 3.0 for windows doesn't support cluster mode
vosegaard
Send Email Send Email
 
There was a significant time cost in the cluster implementation, so in version 3.0 we use MPI for cluster purposes instead of the old custom cluster functionality.

Thomas

On Mar 12, 2011, at 03:23 , Bingwen Hu wrote:

 

Dear developper,

     It looks that simpson 3.0 for windows doesn't support cluster mode any more. Why? SIMPSON 1.1 supports.

Best
Bingwen




#921 From: "enrico_ravera" <eric.ravera@...>
Date: Thu Mar 17, 2011 1:56 pm
Subject: setting spectral width in ROCSA simulation
enrico_ravera
Send Email Send Email
 
Hi everyone,

I'm sorry for the very stupid question but I'm stuck and I need some helping
hand...
  I'm simulating ROCSA on carbonyls (basing on the tar file downloaded from this
group) and I'm getting folded profiles. How can I increase the spectral width
and get rid of this problem?
(i tried increasing the sw parameter but it seems not to affect anything)

regards

E.

#922 From: Julien TREBOSC <simpson.trebosc@...>
Date: Mon Mar 21, 2011 6:55 pm
Subject: Re: Re: Syntax for manual propagation
jtrebosc
Send Email Send Email
 
Another possibility is to apply a filter that select all elements. This forces density matrix propagation and does not need the actual recording of a point.
for example :

matrix set 1 notelements {}
# select all elements in density
filter 1
matrix get density

Julien TREBOSC

Le 12/01/2011 21:17, kjharris_ssnmr a écrit :
 

Thanks for the quick reply Thomas, and that clears that issue up.

Kris

--- In simpson-simmol@yahoogroups.com, Thomas Vosegaard <tv@...> wrote:
>
> Hi,
>
> The issue you're observing is due to the fact that when propagating a pulse sequence, only the propagator is calculated. To save matrix multiplications, the propagator is only applied to evolve the density operator when the "acq" command is invoked. Thus, if you are using "matrix get density" for debugging, you should always put in an "acq" command before.
>
> By the way, the sequence
>
> > pulseid 8 62500 x
> > store 5
> > prop 5
> >
>
> corresponds to two 180 deg. pulses, one by the pulseid command, the second from "prop 5".
>
> Thomas
>
> On Jan 12, 2011, at 20:21 , kjharris_ssnmr wrote:
>
> > Hello all,
> >
> > I would like to follow changes in the density matrix at various points using the "putmatrix" command, but I'm having trouble. I can print out the current density matrix, but I can't control when it is propagated forward in time by the pulses/delays. I'm trying to use the "prop" command to do this, but can't figure the syntax.
> >
> > e.g. a 180 degree pulse on a single 1H nucleus:
> > -------------------------------------
> > spinsys {
> > channels 1H
> > nuclei 1H
> > }
> >
> > par {
> > method direct
> > spin_rate 0
> > proton_frequency 400.0e6
> > crystal_file alpha0beta0
> > start_operator I1z
> > detect_operator I1p
> > }
> >
> > proc pulseq {} {
> > reset
> >
> > putmatrix [matrix get density] #Just print starting p, = Iz
> >
> > pulseid 8 62500 x
> > store 5
> > prop 5
> > putmatrix [matrix get density] #Should print p after evolving with pulse
> >
> > acq
> > putmatrix [matrix get density] #Print p after acq command to see the difference
> > }
> >
> > proc main {} {
> > set f [fsimpson]
> > funload $f
> > }
> > -------------------------------------
> >
> > The result is:
> > -------------------------------------
> > ( 0.5, 0) ( 0, 0)
> > ( 0, 0) ( -0.5, 0)
> >
> > ( 0.5, 0) ( 0, 0)
> > ( 0, 0) ( -0.5, 0)
> >
> > ( -0.5, 0) ( 0, 2.83e-16)
> > ( 0,-2.83e-16) ( 0.5, 0)
> > -------------------------------------
> >
> > So the propagator doesn't alter the current density matrix after the "prop" command, only
> > after the "acq" command is called
> >
> > thanks for any tips,
> > Kris
> >
> >
>



#923 From: Julien TREBOSC <simpson.trebosc@...>
Date: Mon Mar 21, 2011 7:12 pm
Subject: Re: Quadrupolar Echo Simulation
jtrebosc
Send Email Send Email
 
There is a long standing bug in simpson since they introduced the ability of changing the MAS axis during the sequence (maybe not related to the bug). All versions 1.2 and later (2, 3) have this bug.

Check this message to make a test:
http://tech.groups.yahoo.com/group/simpson-simmol/message/815?gi=0

Basically the time origin is wrong when using quadrupole to the second order. This makes impossible to simulate more than one pulse.

Personnally I use simpson 1.1.1 whenever I simulate quadrupolar nuclei. However there is a lot of work to clean the code for compilation with modern API.

I have an update version. If you need it I can send the archive. However many things are disabled (minuit, simmol) and I tested it only under linux.

Best,

Julien TREBOSC

Le 27/01/2011 18:24, kjharris_ssnmr a écrit :
 

Hello SIMPSON users,

Hopefully you guys can spot what I'm doing wrong here, I'm a little confused. When simulating an echo in a 14N spectrum, it only works when the quadrupolar interaction is set to 1st order, not 2nd. I get a weird lineshape that can't be phased in the second case.

With direct acquisition (no echo) the exact same lineshape results at 1st/2nd order (as expected, the 2nd-order contribution for spin-1 nucleus isn't much against the large 1st-order). So... if the lineshape is the same, I don't think the 2nd order perturbation should affect the formation of an echo.

Anybody know what might be going on here?

Here's the simulation code I used, where I (Start as Ix)-(delay 50 &#956;s)-(ideal &#960;/2 pulse)-(delay 50 &#956;s)-observe I+.
This version works fine, but fails if I swap in this line:
quadrupole 1 2 1.14e6 0.3 0.0 0.0 0.0

-----------------
spinsys {
channels 14N
nuclei 14N
shift 1 0.0 0.0 0.0 0.0 0.0 0.0
quadrupole 1 1 1.14e6 0.3 0.0 0.0 0.0
}

par {
method direct
spin_rate 0
sw 12.0e6
np 512
proton_frequency 400.0e6
crystal_file zcw28656.cry
start_operator I1x
detect_operator I1p
verbose 1111
}

proc pulseq {} {
global par
set dw [expr 1.0e6/$par(sw)]

set d1 50.0

reset
delay $dw
store 1

reset
delay $d1
pulseid 4 62500 x
delay $d1

acq $par(np) 1

}

proc main {} {
global par
set f [fsimpson]
fsave $f $par(name).fid
fft $f
fsave $f $par(name).spe
funload $f
}



#924 From: "vosegaard" <tv@...>
Date: Fri Mar 25, 2011 3:08 pm
Subject: SIMPSON 3.1.0 released
vosegaard
Send Email Send Email
 
Dear all

It is a pleasure for me to announce the release of SIMPSON version 3.1.0. This is a bug fix release that fixes a few bugs in the OC routines and for quadrupolar nuclei with second-order effects. Furthermore, the Makefile has been shined up a little for those compiling SIMPSON on their own. Otherwise, just download the installers.

http://www.bionmr.chem.au.dk/simpson
 
Thomas

#925 From: "Lipton, Andrew S" <as_lipton@...>
Date: Fri Mar 25, 2011 6:35 pm
Subject: Re: SIMPSON 3.1.0 released
aslipton
Send Email Send Email
 
Thomas,

on the mac version I still can't get it to pass the frms call. It just returns
"Abort" and stops.

Andy

On Mar 25, 2011, at 8:08 AM, vosegaard wrote:



Dear all

It is a pleasure for me to announce the release of SIMPSON version 3.1.0. This
is a bug fix release that fixes a few bugs in the OC routines and for
quadrupolar nuclei with second-order effects. Furthermore, the Makefile has been
shined up a little for those compiling SIMPSON on their own. Otherwise, just
download the installers.

<http://www.bionmr.chem.au.dk/simpson>
http://www.bionmr.chem.au.dk/simpson

Thomas

#926 From: Bingwen Hu <bw.hu@...>
Date: Mon Apr 11, 2011 4:18 am
Subject: Compiler Error with INTEL MKL 11.1
bw.hu
Send Email Send Email
 
Dear ALL,
  Does anyone have met this problem when compiling SIMPSON 3.1 with INTEL MKL 11.1:

/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_in_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_critical'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_ordered'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_atomic_fixed4_add'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_num_threads'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_fini_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_fini'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `ompc_set_nested'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_serialized_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_atomic_fixed8_add'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_reduce_nowait'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_ordered'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_fini_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_nested'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_thread_num'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_next_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_init_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_num_procs'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_flush'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_barrier'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_reduce_nowait'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_serialized_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_init_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_single'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_critical'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_init_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_fork_call'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_init_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_next_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_push_num_threads'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_single'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_global_thread_num'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_ok_to_fork'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_max_threads'
make: *** [simpson] Error 1


Conditions:
INCLUDES = -I/usr/include/ -I/data/soft/compiler/intel/Compiler/11.1/056/mkl/include/ -I/data/soft/compiler/mpi/openmpi/1.3.3/gcc.pgf90/include/

LIBRARIES = -L/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64 -L/data/soft/compiler/mpi/openmpi/1.3.3/gcc.pgf90/lib -lmkl_intel_lp64 -lmkl_core -lmkl_lapack -fopenmp -ltcl8.4 -lmkl_intel_thread -lm

CC =gcc
EXTRA_FLAGS = -DMKL_INTEL -DMPI

  Looking forward to your reply.

Best
Bingwen




#927 From: Zdenek Tosner <tosner@...>
Date: Mon Apr 11, 2011 6:22 am
Subject: Re: Compiler Error with INTEL MKL 11.1
ztosner
Send Email Send Email
 
Hello,
this looks like you are missing some MKL library. Try to add -lguide and let us know.
Zdenek

Dne 11.4.2011 6:18, Bingwen Hu napsal(a):
 
Dear ALL,
  Does anyone have met this problem when compiling SIMPSON 3.1 with INTEL MKL 11.1:

/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_in_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_critical'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_ordered'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_atomic_fixed4_add'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_num_threads'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_fini_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_fini'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `ompc_set_nested'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_serialized_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_atomic_fixed8_add'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_reduce_nowait'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_ordered'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_fini_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_nested'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_thread_num'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_next_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_init_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_num_procs'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_flush'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_barrier'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_reduce_nowait'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_serialized_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_init_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_single'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_critical'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_init_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_fork_call'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_init_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_next_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_push_num_threads'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_single'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_global_thread_num'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_ok_to_fork'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_max_threads'
make: *** [simpson] Error 1


Conditions:
INCLUDES = -I/usr/include/ -I/data/soft/compiler/intel/Compiler/11.1/056/mkl/include/ -I/data/soft/compiler/mpi/openmpi/1.3.3/gcc.pgf90/include/

LIBRARIES = -L/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64 -L/data/soft/compiler/mpi/openmpi/1.3.3/gcc.pgf90/lib -lmkl_intel_lp64 -lmkl_core -lmkl_lapack -fopenmp -ltcl8.4 -lmkl_intel_thread -lm

CC =gcc
EXTRA_FLAGS = -DMKL_INTEL -DMPI

  Looking forward to your reply.

Best
Bingwen





#928 From: Bingwen Hu <bw.hu@...>
Date: Mon Apr 11, 2011 7:54 am
Subject: Re: Compiler Error with INTEL MKL 11.1
bw.hu
Send Email Send Email
 

Hello, Zdenek,

  Strange, when I add " -lguide", the compile of SIMPSON successes and the SIMPSON works well now. I don't know why.
  Thank you very much.

Best
Bingwen


From: Zdenek Tosner <tosner@...>
To: simpson-simmol@yahoogroups.com
Sent: Mon, April 11, 2011 2:22:16 PM
Subject: Re: [simpson-simmol] Compiler Error with INTEL MKL 11.1

 

Hello,
this looks like you are missing some MKL library. Try to add -lguide and let us know.
Zdenek

Dne 11.4.2011 6:18, Bingwen Hu napsal(a):

 
Dear ALL,
  Does anyone have met this problem when compiling SIMPSON 3.1 with INTEL MKL 11.1:

/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_in_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_critical'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_ordered'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_atomic_fixed4_add'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_num_threads'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_fini_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_fini'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `ompc_set_nested'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_serialized_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_atomic_fixed8_add'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_reduce_nowait'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_ordered'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_fini_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_nested'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_thread_num'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_next_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_init_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_num_procs'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_flush'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_barrier'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_reduce_nowait'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_serialized_parallel'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_init_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_single'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_critical'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_init_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_fork_call'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_for_static_init_8'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_dispatch_next_4'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_push_num_threads'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_end_single'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_global_thread_num'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `__kmpc_ok_to_fork'
/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64t/libmkl_intel_thread.so: undefined reference to `omp_get_max_threads'
make: *** [simpson] Error 1


Conditions:
INCLUDES = -I/usr/include/ -I/data/soft/compiler/intel/Compiler/11.1/056/mkl/include/ -I/data/soft/compiler/mpi/openmpi/1.3.3/gcc.pgf90/include/

LIBRARIES = -L/data/soft/compiler/intel/Compiler/11.1/056/mkl/lib/em64 -L/data/soft/compiler/mpi/openmpi/1.3.3/gcc.pgf90/lib -lmkl_intel_lp64 -lmkl_core -lmkl_lapack -fopenmp -ltcl8.4 -lmkl_intel_thread -lm

CC =gcc
EXTRA_FLAGS = -DMKL_INTEL -DMPI

  Looking forward to your reply.

Best
Bingwen






#929 From: "Kris Harris" <kjharris_ssnmr@...>
Date: Wed Apr 20, 2011 6:54 pm
Subject: Re: Quadrupolar Echo Simulation
kjharris_ssnmr
Send Email Send Email
 
Thanks very much for the info, I wasn't aware of that bug, but I'll keep it in
mind now.

thanks,
Kris


--- In simpson-simmol@yahoogroups.com, Julien TREBOSC <simpson.trebosc@...>
wrote:
>
> There is a long standing bug in simpson since they introduced the
> ability of changing the MAS axis during the sequence (maybe not related
> to the bug). All versions 1.2 and later (2, 3) have this bug.
>
> Check this message to make a test:
> http://tech.groups.yahoo.com/group/simpson-simmol/message/815?gi=0
>
> Basically the time origin is wrong when using quadrupole to the second
> order. This makes impossible to simulate more than one pulse.
>
> Personnally I use simpson 1.1.1 whenever I simulate quadrupolar nuclei.
> However there is a lot of work to clean the code for compilation with
> modern API.
>
> I have an update version. If you need it I can send the archive. However
> many things are disabled (minuit, simmol) and I tested it only under linux.
>
> Best,
>
> Julien TREBOSC
>
> Le 27/01/2011 18:24, kjharris_ssnmr a écrit :
> >
> > Hello SIMPSON users,
> >
> > Hopefully you guys can spot what I'm doing wrong here, I'm a little
> > confused. When simulating an echo in a 14N spectrum, it only works
> > when the quadrupolar interaction is set to 1st order, not 2nd. I get a
> > weird lineshape that can't be phased in the second case.
> >
> > With direct acquisition (no echo) the exact same lineshape results at
> > 1st/2nd order (as expected, the 2nd-order contribution for spin-1
> > nucleus isn't much against the large 1st-order). So... if the
> > lineshape is the same, I don't think the 2nd order perturbation should
> > affect the formation of an echo.
> >
> > Anybody know what might be going on here?
> >
> > Here's the simulation code I used, where I (Start as Ix)-(delay 50
> > μs)-(ideal π/2 pulse)-(delay 50 μs)-observe I+.
> > This version works fine, but fails if I swap in this line:
> > quadrupole 1 2 1.14e6 0.3 0.0 0.0 0.0
> >
> > -----------------
> > spinsys {
> > channels 14N
> > nuclei 14N
> > shift 1 0.0 0.0 0.0 0.0 0.0 0.0
> > quadrupole 1 1 1.14e6 0.3 0.0 0.0 0.0
> > }
> >
> > par {
> > method direct
> > spin_rate 0
> > sw 12.0e6
> > np 512
> > proton_frequency 400.0e6
> > crystal_file zcw28656.cry
> > start_operator I1x
> > detect_operator I1p
> > verbose 1111
> > }
> >
> > proc pulseq {} {
> > global par
> > set dw [expr 1.0e6/$par(sw)]
> >
> > set d1 50.0
> >
> > reset
> > delay $dw
> > store 1
> >
> > reset
> > delay $d1
> > pulseid 4 62500 x
> > delay $d1
> >
> > acq $par(np) 1
> >
> > }
> >
> > proc main {} {
> > global par
> > set f [fsimpson]
> > fsave $f $par(name).fid
> > fft $f
> > fsave $f $par(name).spe
> > funload $f
> > }
> >
> >
>

#930 From: Thomas Vosegaard <tv@...>
Date: Wed Apr 20, 2011 7:33 pm
Subject: Re: Re: Quadrupolar Echo Simulation
vosegaard
Send Email Send Email
 
Could you try to run it with the latest release 3.1.0 and see if the problem persists? In this release, we have fixed a bug related to quadrupolar nuclei, so hopefully your problem is also solved.

Thomas


On Apr 20, 2011, at 20:54 , Kris Harris wrote:

 

Thanks very much for the info, I wasn't aware of that bug, but I'll keep it in mind now.

thanks,
Kris

--- In simpson-simmol@yahoogroups.com, Julien TREBOSC <simpson.trebosc@...> wrote:
>
> There is a long standing bug in simpson since they introduced the
> ability of changing the MAS axis during the sequence (maybe not related
> to the bug). All versions 1.2 and later (2, 3) have this bug.
>
> Check this message to make a test:
> http://tech.groups.yahoo.com/group/simpson-simmol/message/815?gi=0
>
> Basically the time origin is wrong when using quadrupole to the second
> order. This makes impossible to simulate more than one pulse.
>
> Personnally I use simpson 1.1.1 whenever I simulate quadrupolar nuclei.
> However there is a lot of work to clean the code for compilation with
> modern API.
>
> I have an update version. If you need it I can send the archive. However
> many things are disabled (minuit, simmol) and I tested it only under linux.
>
> Best,
>
> Julien TREBOSC
>
> Le 27/01/2011 18:24, kjharris_ssnmr a écrit :
> >
> > Hello SIMPSON users,
> >
> > Hopefully you guys can spot what I'm doing wrong here, I'm a little
> > confused. When simulating an echo in a 14N spectrum, it only works
> > when the quadrupolar interaction is set to 1st order, not 2nd. I get a
> > weird lineshape that can't be phased in the second case.
> >
> > With direct acquisition (no echo) the exact same lineshape results at
> > 1st/2nd order (as expected, the 2nd-order contribution for spin-1
> > nucleus isn't much against the large 1st-order). So... if the
> > lineshape is the same, I don't think the 2nd order perturbation should
> > affect the formation of an echo.
> >
> > Anybody know what might be going on here?
> >
> > Here's the simulation code I used, where I (Start as Ix)-(delay 50
> > &#956;s)-(ideal &#960;/2 pulse)-(delay 50 &#956;s)-observe I+.
> > This version works fine, but fails if I swap in this line:
> > quadrupole 1 2 1.14e6 0.3 0.0 0.0 0.0
> >
> > -----------------
> > spinsys {
> > channels 14N
> > nuclei 14N
> > shift 1 0.0 0.0 0.0 0.0 0.0 0.0
> > quadrupole 1 1 1.14e6 0.3 0.0 0.0 0.0
> > }
> >
> > par {
> > method direct
> > spin_rate 0
> > sw 12.0e6
> > np 512
> > proton_frequency 400.0e6
> > crystal_file zcw28656.cry
> > start_operator I1x
> > detect_operator I1p
> > verbose 1111
> > }
> >
> > proc pulseq {} {
> > global par
> > set dw [expr 1.0e6/$par(sw)]
> >
> > set d1 50.0
> >
> > reset
> > delay $dw
> > store 1
> >
> > reset
> > delay $d1
> > pulseid 4 62500 x
> > delay $d1
> >
> > acq $par(np) 1
> >
> > }
> >
> > proc main {} {
> > global par
> > set f [fsimpson]
> > fsave $f $par(name).fid
> > fft $f
> > fsave $f $par(name).spe
> > funload $f
> > }
> >
> >
>



#931 From: "Kris Harris" <kjharris_ssnmr@...>
Date: Wed Apr 27, 2011 8:45 pm
Subject: Re: Quadrupolar Echo Simulation
kjharris_ssnmr
Send Email Send Email
 
Hi Thomas,

I just downloaded and tried the new version (3.1.0 for mac), but I get the same
result, so I guess that wasn't the issue.  Here's a link of the results of
straight observation versus echo for 1st and 2nd order quadrupolar interactions:
http://img684.imageshack.us/i/picture10gm.png/

cheers,
Kris


--- In simpson-simmol@yahoogroups.com, Thomas Vosegaard <tv@...> wrote:
>
> Could you try to run it with the latest release 3.1.0 and see if the problem
persists? In this release, we have fixed a bug related to quadrupolar nuclei, so
hopefully your problem is also solved.
>
> Thomas
>
>
> On Apr 20, 2011, at 20:54 , Kris Harris wrote:
>
> > Thanks very much for the info, I wasn't aware of that bug, but I'll keep it
in mind now.
> >
> > thanks,
> > Kris
> >
> > --- In simpson-simmol@yahoogroups.com, Julien TREBOSC <simpson.trebosc@>
wrote:
> > >
> > > There is a long standing bug in simpson since they introduced the
> > > ability of changing the MAS axis during the sequence (maybe not related
> > > to the bug). All versions 1.2 and later (2, 3) have this bug.
> > >
> > > Check this message to make a test:
> > > http://tech.groups.yahoo.com/group/simpson-simmol/message/815?gi=0
> > >
> > > Basically the time origin is wrong when using quadrupole to the second
> > > order. This makes impossible to simulate more than one pulse.
> > >
> > > Personnally I use simpson 1.1.1 whenever I simulate quadrupolar nuclei.
> > > However there is a lot of work to clean the code for compilation with
> > > modern API.
> > >
> > > I have an update version. If you need it I can send the archive. However
> > > many things are disabled (minuit, simmol) and I tested it only under
linux.
> > >
> > > Best,
> > >
> > > Julien TREBOSC
> > >
> > > Le 27/01/2011 18:24, kjharris_ssnmr a écrit :
> > > >
> > > > Hello SIMPSON users,
> > > >
> > > > Hopefully you guys can spot what I'm doing wrong here, I'm a little
> > > > confused. When simulating an echo in a 14N spectrum, it only works
> > > > when the quadrupolar interaction is set to 1st order, not 2nd. I get a
> > > > weird lineshape that can't be phased in the second case.
> > > >
> > > > With direct acquisition (no echo) the exact same lineshape results at
> > > > 1st/2nd order (as expected, the 2nd-order contribution for spin-1
> > > > nucleus isn't much against the large 1st-order). So... if the
> > > > lineshape is the same, I don't think the 2nd order perturbation should
> > > > affect the formation of an echo.
> > > >
> > > > Anybody know what might be going on here?
> > > >
> > > > Here's the simulation code I used, where I (Start as Ix)-(delay 50
> > > > μs)-(ideal π/2 pulse)-(delay 50 μs)-observe I+.
> > > > This version works fine, but fails if I swap in this line:
> > > > quadrupole 1 2 1.14e6 0.3 0.0 0.0 0.0
> > > >
> > > > -----------------
> > > > spinsys {
> > > > channels 14N
> > > > nuclei 14N
> > > > shift 1 0.0 0.0 0.0 0.0 0.0 0.0
> > > > quadrupole 1 1 1.14e6 0.3 0.0 0.0 0.0
> > > > }
> > > >
> > > > par {
> > > > method direct
> > > > spin_rate 0
> > > > sw 12.0e6
> > > > np 512
> > > > proton_frequency 400.0e6
> > > > crystal_file zcw28656.cry
> > > > start_operator I1x
> > > > detect_operator I1p
> > > > verbose 1111
> > > > }
> > > >
> > > > proc pulseq {} {
> > > > global par
> > > > set dw [expr 1.0e6/$par(sw)]
> > > >
> > > > set d1 50.0
> > > >
> > > > reset
> > > > delay $dw
> > > > store 1
> > > >
> > > > reset
> > > > delay $d1
> > > > pulseid 4 62500 x
> > > > delay $d1
> > > >
> > > > acq $par(np) 1
> > > >
> > > > }
> > > >
> > > > proc main {} {
> > > > global par
> > > > set f [fsimpson]
> > > > fsave $f $par(name).fid
> > > > fft $f
> > > > fsave $f $par(name).spe
> > > > funload $f
> > > > }
> > > >
> > > >
> > >
> >
> >
>

#932 From: Julien TREBOSC <simpson.trebosc@...>
Date: Thu May 19, 2011 1:38 pm
Subject: Re: Re: Quadrupolar Echo Simulation
jtrebosc
Send Email Send Email
 
Hi,

I don't think there is problem with simpson. It's most probably a problem of phase cycling the sequence. Evolution during delay is different under 1st and 2nd order interaction. The echo refocuses the interaction for +1 -> -1 pathway. However even a perfect 90 pulse do not transfer 100% of magnetization from +1 -> -1 and more important do not remove all magnetization from -1 coherence order. So if you do not appy phase cycling you still observe evolution for -1 -> -1 pathway which is different under HQ(2) and HQ(1).

If you apply a proper filter then both simulations are equal.

Julien TREBOSC

Le 27/04/2011 22:45, Kris Harris a écrit :
 

Hi Thomas,

I just downloaded and tried the new version (3.1.0 for mac), but I get the same result, so I guess that wasn't the issue. Here's a link of the results of straight observation versus echo for 1st and 2nd order quadrupolar interactions:
http://img684.imageshack.us/i/picture10gm.png/

cheers,
Kris

--- In simpson-simmol@yahoogroups.com, Thomas Vosegaard <tv@...> wrote:
>
> Could you try to run it with the latest release 3.1.0 and see if the problem persists? In this release, we have fixed a bug related to quadrupolar nuclei, so hopefully your problem is also solved.
>
> Thomas
>
>
> On Apr 20, 2011, at 20:54 , Kris Harris wrote:
>
> > Thanks very much for the info, I wasn't aware of that bug, but I'll keep it in mind now.
> >
> > thanks,
> > Kris
> >
> > --- In simpson-simmol@yahoogroups.com, Julien TREBOSC <simpson.trebosc@> wrote:
> > >
> > > There is a long standing bug in simpson since they introduced the
> > > ability of changing the MAS axis during the sequence (maybe not related
> > > to the bug). All versions 1.2 and later (2, 3) have this bug.
> > >
> > > Check this message to make a test:
> > > http://tech.groups.yahoo.com/group/simpson-simmol/message/815?gi=0
> > >
> > > Basically the time origin is wrong when using quadrupole to the second
> > > order. This makes impossible to simulate more than one pulse.
> > >
> > > Personnally I use simpson 1.1.1 whenever I simulate quadrupolar nuclei.
> > > However there is a lot of work to clean the code for compilation with
> > > modern API.
> > >
> > > I have an update version. If you need it I can send the archive. However
> > > many things are disabled (minuit, simmol) and I tested it only under linux.
> > >
> > > Best,
> > >
> > > Julien TREBOSC
> > >
> > > Le 27/01/2011 18:24, kjharris_ssnmr a écrit :
> > > >
> > > > Hello SIMPSON users,
> > > >
> > > > Hopefully you guys can spot what I'm doing wrong here, I'm a little
> > > > confused. When simulating an echo in a 14N spectrum, it only works
> > > > when the quadrupolar interaction is set to 1st order, not 2nd. I get a
> > > > weird lineshape that can't be phased in the second case.
> > > >
> > > > With direct acquisition (no echo) the exact same lineshape results at
> > > > 1st/2nd order (as expected, the 2nd-order contribution for spin-1
> > > > nucleus isn't much against the large 1st-order). So... if the
> > > > lineshape is the same, I don't think the 2nd order perturbation should
> > > > affect the formation of an echo.
> > > >
> > > > Anybody know what might be going on here?
> > > >
> > > > Here's the simulation code I used, where I (Start as Ix)-(delay 50
> > > > &#956;s)-(ideal &#960;/2 pulse)-(delay 50 &#956;s)-observe I+.
> > > > This version works fine, but fails if I swap in this line:
> > > > quadrupole 1 2 1.14e6 0.3 0.0 0.0 0.0
> > > >
> > > > -----------------
> > > > spinsys {
> > > > channels 14N
> > > > nuclei 14N
> > > > shift 1 0.0 0.0 0.0 0.0 0.0 0.0
> > > > quadrupole 1 1 1.14e6 0.3 0.0 0.0 0.0
> > > > }
> > > >
> > > > par {
> > > > method direct
> > > > spin_rate 0
> > > > sw 12.0e6
> > > > np 512
> > > > proton_frequency 400.0e6
> > > > crystal_file zcw28656.cry
> > > > start_operator I1x
> > > > detect_operator I1p
> > > > verbose 1111
> > > > }
> > > >
> > > > proc pulseq {} {
> > > > global par
> > > > set dw [expr 1.0e6/$par(sw)]
> > > >
> > > > set d1 50.0
> > > >
> > > > reset
> > > > delay $dw
> > > > store 1
> > > >
> > > > reset
> > > > delay $d1
> > > > pulseid 4 62500 x
> > > > delay $d1
> > > >
> > > > acq $par(np) 1
> > > >
> > > > }
> > > >
> > > > proc main {} {
> > > > global par
> > > > set f [fsimpson]
> > > > fsave $f $par(name).fid
> > > > fft $f
> > > > fsave $f $par(name).spe
> > > > funload $f
> > > > }
> > > >
> > > >
> > >
> >
> >
>



#933 From: "vosegaard" <tv@...>
Date: Mon May 23, 2011 7:44 am
Subject: Optimization/Minimization package
vosegaard
Send Email Send Email
 
Hi all

As you may have noticed, we removed the minuit functions in version 3.0 because they caused instabilities and errors.

I have just released an optimization package that may be included in SIMPSON and that allows to optimize a function (e.g., rms deviation between experimental and simulated spectrum), scan parameters and calculate 95% confidence intervals.

The optimization package can be downloaded here


An example of fitting an experimental spectrum can be found here


I hope you find this package useful.

Thomas

#934 From: "rs_comp_chem" <roel.sanchez1@...>
Date: Fri Aug 19, 2011 12:59 pm
Subject: NMR spectrum with large quadrupolar coupling constants > 15 MHz
rs_comp_chem
Send Email Send Email
 
Dear all,

I'm wondering if any one might have a SIMPSON input file to simulate the
spectrum of a system with a very large quadrupolar coupling constant > 15 MHz;
something like:

-----------
quadrupole       1 2 18e6 0.9 0 0 0
-----------

Also, if one ends up having large quadrupolar coupling constants from a
Density-Functional-Theory simulation of the NMR parameters of a molecular solid,
what NMR technique is the most appropriate to get significant NMR information
with the SIMPSON code?

Every time that I change the quadrupolar coupling constant to something like 16
or 18 MHz, my NMR gets too broad. So, I'm curious to know if people in the
community have done similar simulations and will be willing to share their
SIMPSON input files and/or provide some feedback as to how to simulate/predict
the NMR spectrum of system with very large quadrupolar coupling constants.

Thanks,
rs

#935 From: evgeny.markhasin
Date: Mon Aug 29, 2011 5:01 pm
Subject: SIMPSON for Windows Filters Backslashes from the Path
evgeny.markh...
 
Hello,

SIMPSON for windows filters backslashes from command line parameters, which is
improper behavior on windows. Please fix it.

#936 From: Zdenek Tosner <tosner@...>
Date: Tue Aug 30, 2011 7:35 am
Subject: Re: SIMPSON for Windows Filters Backslashes from the Path
ztosner
Send Email Send Email
 
Can you please give me an example? (both correct behavior and improper behavior) Command line parameters are not receiving special treatment inside simpson, it is operating system and Tcl job...

Zdenek

On 08/29/2011 07:01 PM, evgeny.markhasin wrote:
 

Hello,

SIMPSON for windows filters backslashes from command line parameters, which is improper behavior on windows. Please fix it.



#937 From: evgeny.markhasin
Date: Wed Aug 31, 2011 2:35 pm
Subject: Re: SIMPSON for Windows Filters Backslashes from the Path
evgeny.markh...
 
Say, I am trying to start a simulation like this:
    ?>simpson.exe C:\DATA\SEDOR\SEDOR.in
Instead of starting the simulation, it issues an error:
    error: unable to open input file 'C:DATASEDORSEDOR.in',
where all backslashes have been filtered out. This is what one would expect on
Linux, but should never been done on Windows.



--- In simpson-simmol@yahoogroups.com, Zdenek Tosner <tosner@...> wrote:
>
> Can you please give me an example? (both correct behavior and improper
> behavior) Command line parameters are not receiving special treatment
> inside simpson, it is operating system and Tcl job...
>
> Zdenek
>
> On 08/29/2011 07:01 PM, evgeny.markhasin wrote:
> >
> > Hello,
> >
> > SIMPSON for windows filters backslashes from command line parameters,
> > which is improper behavior on windows. Please fix it.
> >
> >
>

#938 From: Zdenek Tosner <tosner@...>
Date: Wed Aug 31, 2011 7:12 pm
Subject: Re: Re: SIMPSON for Windows Filters Backslashes from the Path
ztosner
Send Email Send Email
 
I see, this somehow managed to find its way through our tests... Will be fixed in new version coming out soon.
Zdenek

Dne 31.8.2011 16:35, evgeny.markhasin napsal(a):
 

Say, I am trying to start a simulation like this:
?>simpson.exe C:\DATA\SEDOR\SEDOR.in
Instead of starting the simulation, it issues an error:
error: unable to open input file 'C:DATASEDORSEDOR.in',
where all backslashes have been filtered out. This is what one would expect on Linux, but should never been done on Windows.

--- In simpson-simmol@yahoogroups.com, Zdenek Tosner <tosner@...> wrote:
>
> Can you please give me an example? (both correct behavior and improper
> behavior) Command line parameters are not receiving special treatment
> inside simpson, it is operating system and Tcl job...
>
> Zdenek
>
> On 08/29/2011 07:01 PM, evgeny.markhasin wrote:
> >
> > Hello,
> >
> > SIMPSON for windows filters backslashes from command line parameters,
> > which is improper behavior on windows. Please fix it.
> >
> >
>



#939 From: evgeny.markhasin
Date: Mon Sep 5, 2011 3:14 pm
Subject: Updated Simplex.c file
evgeny.markh...
 
Hello,

I have modified the "Simplex.c" file and emailed the updated version to
Professor Vosegaard. I was wondering if my message made it through.

Thanks,
Evgeny Markhasin

#940 From: evgeny.markhasin
Date: Wed Sep 7, 2011 3:50 am
Subject: Decetion of multiple individual operators
evgeny.markh...
 
If I have, say, {13C 13C 1H} spin system. Is it possible to detect
simultaneously multiple operators, e.g. I1p, I2p, I3x, rather than just their
sum?

#941 From: Thomas Vosegaard <tv@...>
Date: Wed Sep 7, 2011 3:58 am
Subject: Re: Decetion of multiple individual operators
vosegaard
Send Email Send Email
 
Hi,

I have used the following to achieve simultaneous detection of several operators:

set oper {I1x I1y I1z}
set noper [llength $oper]

set npeff [expr $par(np)/$noper]

for {set i 1} {$i <= $npeff} {incr i} {
  # Do the pulse program
  ...
  foreach o $oper {
    matrix set detect operator $o
    acq
  }
}

then you just need to specify np three times as long and separate out the fid points afterwards. Try to see the attached file.

Thomas


On Sep 7, 2011, at 05:50 , evgeny.markhasin wrote:

 

If I have, say, {13C 13C 1H} spin system. Is it possible to detect simultaneously multiple operators, e.g. I1p, I2p, I3x, rather than just their sum?



1 of 1 File(s)


#942 From: Thomas Vosegaard <tv@...>
Date: Wed Sep 7, 2011 4:01 am
Subject: Re: Decetion of multiple individual operators [1 Attachment]
vosegaard
Send Email Send Email
 
To run the attached file you may replace the line

source ccacb.spinsys

by this spin system

spinsys {
#      1      2      3 
#     9C    8CA   22CB 
#
 channels 13C 
 nuclei   13C 13C 13C 
 shift 1 170p -76p 0.9 -71.727 78.37 -103.63
 shift 2 50p -20p 0.43 39.857 18.335 127.98
 shift 3 30p 10p 0 -98.452 122.65 84.753
 dipole 2 1 -2142.4 0 133.91 177.8
 dipole 2 3 -2159.4 0 57.347 -95.247
 dipole 1 3 -489.19 0 41.169 -53.333
 jcoupling 2 1 55 0 0 0 0 0
 jcoupling 2 3 35 0 0 0 0 0
}

Thomas

On Sep 7, 2011, at 05:58 , Thomas Vosegaard wrote:

Hi,

I have used the following to achieve simultaneous detection of several operators:

set oper {I1x I1y I1z}
set noper [llength $oper]

set npeff [expr $par(np)/$noper]

for {set i 1} {$i <= $npeff} {incr i} {
  # Do the pulse program
  ...
  foreach o $oper {
    matrix set detect operator $o
    acq
  }
}

then you just need to specify np three times as long and separate out the fid points afterwards. Try to see the attached file.

Thomas

<c7-2.in>

On Sep 7, 2011, at 05:50 , evgeny.markhasin wrote:

 

If I have, say, {13C 13C 1H} spin system. Is it possible to detect simultaneously multiple operators, e.g. I1p, I2p, I3x, rather than just their sum?




#943 From: evgeny.markhasin
Date: Sun Sep 18, 2011 8:11 pm
Subject: Is Cosine Transform available?
evgeny.markh...
 
Can SIMPSON do cosine transform of real data?

#944 From: Thomas Vosegaard <tv@...>
Date: Mon Sep 19, 2011 4:16 am
Subject: Re: Is Cosine Transform available?
vosegaard
Send Email Send Email
 
No, it does not have a cosine transformation implemented, but if you take real data and zero the imaginary part, the builtin FT routine should give you the cosine transformation as the real part of the transformed spectrum.

Thomas

On Sep 18, 2011, at 22:11 , evgeny.markhasin wrote:

 

Can SIMPSON do cosine transform of real data?



#945 From: evgeny.markhasin
Date: Tue Oct 11, 2011 12:45 pm
Subject: gcompute and "delay 9999"
evgeny.markh...
 
Could you please comment on implicit acquisition employed by gcompute?

Specifically, how does acquisition occur when the pulse sequence only consists
of something like "delay 9999" Does it start at t=0? When I have other pulses
and delays, when does acquisition start?

Also if I do not have "delay 9999" and leave the pulse sequence empty, all fid
points will be the same and equal to the first point of fid detected with "delay
9999" in place. Why does it happen?

What is significance of "9999" in "delay 9999"? Can I pick a different number? I
need to understand where this number comes from, what the logic is behind it.

Thank you very much.

#946 From: evgeny.markhasin
Date: Tue Oct 11, 2011 1:12 pm
Subject: gcompute and "delay 9999" more
evgeny.markh...
 
Hello, I come up with a couple more questions.

If I put "delay 20" in the pulse sequence, the result is the same as with "delay
9999". If I put "delay 2" or two delays "delay 10", the results will be
different. Could you please explain this behavior and what the different results
mean? (My dwell time is 10 usec.)

Also, if I want to have a specific propagator during acquisition, e.g. I want to
apply decoupling or disable some terms of the internal Hamiltonian, can I do it
with gcompute and how?

#947 From: Zdenek Tosner <tosner@...>
Date: Tue Oct 11, 2011 1:40 pm
Subject: Re: gcompute and "delay 9999" more
ztosner
Send Email Send Email
 
I agree, this is very confusing. My suggestion is to forget about delay 9999 and use the current version of simpson with igcompute method described on the web

http://www.bionmr.chem.au.dk/bionmr/software/simpson.php

We are preparing new version with somewhat enhanced fool-proof checks and different inner code design,  but will be grateful for your experience with the acq_block command! Here you can put your decoupling or Hamiltonian manipulation - please give me feedback...

Zdenek T.
 
Dne 11.10.2011 15:12, evgeny.markhasin napsal(a):
 

Hello, I come up with a couple more questions.

If I put "delay 20" in the pulse sequence, the result is the same as with "delay 9999". If I put "delay 2" or two delays "delay 10", the results will be different. Could you please explain this behavior and what the different results mean? (My dwell time is 10 usec.)

Also, if I want to have a specific propagator during acquisition, e.g. I want to apply decoupling or disable some terms of the internal Hamiltonian, can I do it with gcompute and how?



#948 From: evgeny.markhasin
Date: Mon Oct 17, 2011 2:32 am
Subject: Re: gcompute and "delay 9999" more
evgeny.markh...
 
Thank you for your suggestion. I have actually tried the new igcompute a little
bit. This new version is indeed very intuitive and and allows perform
simulations, which, as far as I understand, are not possible using the old
gcompute. So far, I have not experienced any issues. I will give more feedback
as I run more simulations.

I understand that the new modes are in the development. I am missing the ability
to run a 2D experiment, thus, having acq_block{} in a cycle. Also, I wished I
could use multiple acq_block{}, with each acquiring a specified number of points
like "acq n". This feature would be convenient for CPMG experiment, though, I
understand that this application may not be sufficiently important.

--- In simpson-simmol@yahoogroups.com, Zdenek Tosner <tosner@...> wrote:
>
> I agree, this is very confusing. My suggestion is to forget about delay
> 9999 and use the current version of simpson with igcompute method
> described on the web
>
> http://www.bionmr.chem.au.dk/bionmr/software/simpson.php
>
> We are preparing new version with somewhat enhanced fool-proof checks
> and different inner code design,  but will be grateful for your
> experience with the acq_block command! Here you can put your decoupling
> or Hamiltonian manipulation - please give me feedback...
>
> Zdenek T.
>
> Dne 11.10.2011 15:12, evgeny.markhasin napsal(a):
> >
> > Hello, I come up with a couple more questions.
> >
> > If I put "delay 20" in the pulse sequence, the result is the same as
> > with "delay 9999". If I put "delay 2" or two delays "delay 10", the
> > results will be different. Could you please explain this behavior and
> > what the different results mean? (My dwell time is 10 usec.)
> >
> > Also, if I want to have a specific propagator during acquisition, e.g.
> > I want to apply decoupling or disable some terms of the internal
> > Hamiltonian, can I do it with gcompute and how?
> >
> >
>

Messages 919 - 948 of 1016   Oldest  |  < Older  |  Newer >  |  Newest
Add to My Yahoo!      XML What's This?

Copyright © 2010 Yahoo! Inc. All rights reserved.
Privacy Policy - Terms of Service - Guidelines NEW - Help