Skip to content
Snippets Groups Projects
Commit 2f856c75 authored by Paul McCarthy's avatar Paul McCarthy :mountain_bicyclist:
Browse files

TEST: Ported runme_gpu scripts to feedsRun equivalents

parent 0b2d8ad6
No related branches found
No related tags found
1 merge request!32Mnt/chris rorden gpu test
Pipeline #11307 passed
#!/usr/bin/env bash
set -e
outdir=$1
indir=$(pwd)
pushd $outdir
bedpostx_gpu $indir/bedpost
\ No newline at end of file
#!/usr/bin/env bash
set -e
outdir=$1
indir=$(pwd)
pushd ${outdir}
dti=${indir}/DTIap
dti_bvec=${indir}/DTIap.bvec
dti_bval=${indir}/DTIap.bval
dti_b=${outdir}/DTIapb #brain-extracted dti
dti_b0=${outdir}/DTIapb0
dti_u=${outdir}/DTIapu #undistorted dti
dti_b=${outdir}/DTIapb_mask #masked brain-extracted dti
dti_txt=${outdir}DTIap_acq_param.txt
dti_txt2=${outdir}/DTIap_index.txt
#create acq_param: dummies as we will not run TOPUP
printf "0 1 0 0.03388\n0 -1 0 0.03388\n" > $dti_txt
#create index files: all = 1 as we will not run TOPUP
nvol=$(fslnvols $dti)
indx=""
for ((i=1; i<=nvol; i+=1)); do indx="$indx 1"; done
echo $indx > $dti_txt2
fslroi $dti $dti_b0 0 1
bet $dti_b0 $dti_b -f 0.2 -R -n -m
eddy \
--imain=$dti \
--mask=$dti_b --acqp=$dti_txt \
--index=$dti_txt2 \
--bvecs=$dti_bvec \
--bvals=$dti_bval \
--repol \
--out=$dti_u
#!/usr/bin/env bash
set -e
outdir=$1
indir=$(pwd)
pushd $outdir
for i in {1..172}; do
probtrackx2_gpu -x $indir/masks/$i.nii.gz \
--dir=$indir/probtrackx/$i --forcedir -P 200 \
-s $indir/bedpost.bedpostX/merged \
-m $indir/bedpost.bedpostX/nodif_brain_mask \
--opd --pd -l -c 0.2 --distthresh=0"
done
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment